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CHAPITRE XIV. 

LE PROBLEME DE LL FONCTION PERTURBA.TRICE. 


211. On se rappelle quel est le principe de la methode de la 
variation des constantes que nous avons exposee au Ghapitre IV, 
Tome I. On adopte comrae variables independantes l 7 un des sys- 
temes que nous avons appeles systemes de variables kepleriennes , 
par exemple le systeme des variables L, A, £, (Cf. n os 56, 57, 
78, 79). 

On arrive ainsi aux equations (4 bis) du n° 93 

1 dh dFi di] dFi d% dF x 

\ ~di = — ^ 5 

(0 < 

K J ) dX _ dFp dFj 

[ dt ~~ d\> ^ dL 

Les fonctions F 0 et F* ont ete definies aux n os 36 et suivants; la 
premiere est tres simple, la seconde est ce qu’on appelle la fonc- 
tion perturbatrice . 

Nous avons vu aux n os 83 et 99 que eette fonction peut se d£ve- 
lopper sous la forme 

( i ) jjl F j[ = A cos ( k\ A, -+- ^2^2 — {— h ) OTL , 
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ou k { et A*a sent des entiers. ou A et h dependent seulement des L 
et ou est un monome entier par rapport aux ; et aux r 4 . 

D'apres la methode de Lagrange, il faut aior* : 

r' En premiere approximation, regarder touted les variables 
kepleriennes comme des constantes. sauf les A qui seront des 
fonctions lineaires du temps; on aura ainsi 

U) L^L?, 5,= £?- r kl = r*. 

On aura d'ailleurs 


F 




ill 

2hj 



n t t Li)*=M, 


( Cf . n“> 8 2 et 97). 

2 ° En seconde approximation, rempJacer les variables par leurs 
valeurs (3) dans toutes les derivees parti elles de F<, de sorte que 
les termes des equations (0 ou figurent ces derivees de\iennent 
des fonctions connues du temps. Les equations (i) s’integrent alors 
facilement par quadratures. 

3° En troisieme approximation, remplacer dans Jes derivees de 
F i les variables par leurs valeurs de seconde approximation et inte- 
grer de nouveau les equations (i) par quadratures, et ainsi de 
suite. 


Generalemenl la seconde approximation suffit. Si l’on remplace 
la fonction F, par son developpement (2 ) et qu’on procede a cette 
seconde approximation, on trouve tout de suite 

Li = Lj h— — A cos ( k\ X 1 — h Xo -+- h ) e)R, ; 

*1' = *1? -*-y - sin( X-iXi-t- A-.X.+ A)^, 
jmd v " dr^ 


de sorte que le probleme des perturbations planetaires (au moins 
jusqu’a la seconde approximation inclusivement) serait entiere- 
ment resolu, si Ton connaissait le developpement ( 2 ). 

Or nous avons bien montre que le developpement ( 2 ) est pos- 
sible, mais nous n’avons pas fait voir comment on pouvait effective- 
ment 1’obtenir. Ce dernier probleme, dont on comprend aisdment 
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Fimportance, \ a maintenant nous occuper, et je veux seulement, 
dans ce Chapitre, montrer sous quelles formes di\erses il se pre- 
sente. 


212. Nous avons \u au Chapitre II que la fonction perturba- 
trice pou\ait etre mise sous trois formes diflerentes qui sont celle 
des n us 26 et 43. celle des n L,s 30 et 38. celle du n° 4-4. 

Dans ces trois cas, nous avons decompose la fonction perturba- 
trice en deux parties, la partie principcile et la j tar tie cample - 
meat a ire. 

Au n° 38, nous avons ecrit la partie principale sous la forme 


(4) tfij/724 


WxT+ x'i ~ f-T- ( 37 3 




au n° 43, de meme qu’au n° 44, sous la forme 

( 5 j —m i m k — - - =.. 1 ====== = - • 

V ( ^4 — ,r 1 — 6 — ^3 f 2 

11 est aise de passer d 7 une de ces expressions a Fautre; elles ne 
different en efiet que du premier terme de la parenthese (4)* 

Or ce terme ne depend que des coordonnees x \ , x\ de Fune 
des deux planetes fictives; le developpement peut en etre oblenu 
aisement ainsi que nous le verrons des le Chapitre prochain. 

Passons ala partie complementaire de la fonction perturbatrice ; 
si nous adoptons les valuables du n° 30, nous avons trouve au 
n° 38 Fexpression de cette partie complementaire et nous avons 
expose ensuite, dans le meme n° 38, comment on pouvait passer 
du developpement de la partie principale a celui de la fonction 
perturbatrice complete, et, par consequent, a celui de la partie 
complementaire. Nous reviendrons sur ce point. 

Si nous adoptons au contraire les variables du n° 26, la partie 
complementaire a line expression plus simple et elle s’ecrit comme 
nous Favons vu au n° 43 ; 


< 6 ) 


T,= 


i 

m 1 


(y , iy\-yy , « r '»-+‘ y'*y*)- 


Si enfin on adopte les variables usuelles, c J est-a-dire celles du 
n° 44, la partie complementaire n’est plus la meme pour les deux 
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planetes. Pour Tune, elle s'eerit 

m . , , , , 

<7) ****->- ****)• 

et, pourPautre, 

. . nii ni’ . , . , , . , 

(7 ins ) ^ qF ^ 1 ^ ^3^*6 *• 

Nous verrons plus loin que le developpement des expressions* 
(6), (7), (7 bis) est tres aise. 

Le probleme se trouvera done ramene au developpement de 
[’expression (5) et e’est cette question qui nous retiendra le plus 
longtemps. Nous aurons ensuite a voir comment on peut passer du 
developpement de (5) a celui dela partie complementaire sous ses 
diverses formes et quelles relations il y a entre les developpements 
des diverses expressions ( 5 ), (6), (7), (7 bis). 

213 . Nous avons deiini au n° 59 quatre systemes de variables 
kepleriennes ; ce sont : 

[° Le systeme des elements elliptiques 

a. e , Z, /, g -+- 6, 0 ; 

2 0 Le premier systeme d’elements canoniques 
L, G, 0, Z, g, 6 ; 

3 ° Le deuxieme systeme d’elements canoniques 

L, p l7 A, a),; 

4 ° Le troisieme systeme d’elements canoniques 
L, X, y),. 

Selon que Ton adoptera Pun ou l’autre de ces systemes de va- 
luables, le developpement de la fonction perturbatrice sera ^videm- 
ment different. Heureusement il est ais6 de passer d’un d^velop- 
pement a Pautre. 

Si l’6n adopte Pun des systemes d’el^ments canoniques comme 
nous l’avons fait dans tout le Tome I, on a cet avantage que les. 
Equations du probleme conservent la forme canonique. 
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Cependant les astronomes emploient plus communement les ele- 
ments elliptiques. Les equations, alors, perdent leur forme cano- 
nique. Neanmoins. ainsique nous 1'aAons expose au n° 81 . les equa- 
tions se presentent toujoui'S sous la forme suivante : 

Les derivees par rapport a t des elements elliptiques s’exprime- 
ront lineairement a 1 'aide des derivees partielles de F par rapport 
aux elements elliptiques, les coefficients etant des fonctions con- 
nues des elements elliptiques. 

J’ajouterai que dans ces coefficients, fonctions connues des ele- 
ments elliptiques, figui'eront seulement les elements 

a , e , i, £-- 4 - 6 . 6, 

qui sont des constantes en premiere approximation, tandis que 
Lanomalie moyenne /, qui, en premiere approximation, varie pro- 
portionnellement au temps, ne figurera pas. 

Les equations ainsi obtenues pourront s’appeler les equations de 
Lagrange . 

a o 

Dans le cas des variables canoniques et des equations cano- 
niques, nous n’avons dans le second membre qu' une seule des deri- 
vees partielles de la fonction F et elle est a flee tee seulement du 
coefficient ± i ou — i . Au contrame, dans le cas des variables 
-elliptiques et des equations de Lagrange, nous avons, dans le 
second membre, plusieurs dei'ivees partielles de F et elles sont 
alfectees d un coefficient qui depend lui-meme des elements ellip- 
tiques. Nous ne transcrirons pas ici les equations de Lagrange; 
nous nous bornerons, comme au n° 81 , a renvoyer a LOuvrage de 
Tisserand (t. I, p. 187). 

Tout ce que nous avons etabli dans le Tome I, en nous servant 
des equations canoniques, aurait pu evidemment etre demontre en 
partant des equations de Lagrange, mais les demonstrations au- 
raient ete fort allong^es par la longueur des ecritures, sans que 
rien y soit change d’essentiel. De meme, dans la pratique du 
-calcul, l’emploi des equations canoniques epargnerait bien des 
longueurs, mais la difference ne commencerait k devenir tr£s sen- 
sible qu 5 a la troisieme approximation, tandis qu’on s’arrete d’or- 
dinaire k la seconde. 

Quoi qu’il en soit, il arrive que les elements elliptiques etant 



6 


CHAPITRE XIV. 


plus familiers aux astronomes et plus anciennement employes, la 
plupart des travaux relatifs au developpement de la fonction per- 
turbatrice ont ete executes avec ces elements. 11 n’aurait peut-^tre 
pas ete plus difficile de les faire directement en se servant des 
elements eanoniques, si on les avait adoptes tout d’abord. Mais 
aujourd'hui ce serai t refaire un tra\ail qui a deja ete fait une fois 
et il \aut mieux utiliser les resultats obtenus par pos devanciers. 

Heureusement, comme je Pai dit, il est aise de passer du deve- 
loppement procedant suivant les elements elliptiques a celui qui 
procede suivant Tun ou l’autre des systemes d’elements cano- 
niques. 

Nous ferons peu d'usage du premier systeme d’elements cano- 
niques ; en effet, si les excentricites et les inclinaisons sont tres 
petites, il arrive que quelques-unes des variables du deuxieme et 
du troisieme systemes eanoniques sont egalement tres petites. La 
meme chose n ’arrive pas avec le premier systeme, qui est ainsi peu 
applicable a des methodes d’approximation ou la petitesse des 
excentricites et des inclinaisons joue un role capital. 

C’est ce qui nous engage a le rejeter. 

En revanche, nous serons obliges de considerer un autre systeme 
de variables elliptiques 

a , e, i, u , ^4-0, 6, 

ou l’anomalie moyenne l est remplacee par l’anomalie excen- 
trique a. 

L’emploi de ces variables se prete mal a l’integration, bien que 
Hansen en ait fait quelque usage. En revanche, le developpement 
procedant suivant ces nouvelles variables est beaucoup plus facile 
que celui qui procede suivant les variables elliptiques ordinaires. 
Il est aise de passer de Fun a l’autre et c’est par l’intermediaire du 
premier developpement qu’il est le plus commode de parvenir au 
second. C’est ce qui justifie l’etude detaillee que nous ferons de ce 
premier developpement. 

Nous aurons done a etudier quatre developpements de la fonc- 
tion perturbatrice procedant : 

i" Suivant les variables 


( 8 ) 


a, 6, i, u , g -+- 0, 6 
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a 1 ' Suivant les variables elliptiques 

«, e. i, /. g — 0 , fi, 

on, si Ton aime mieux, puisque 

a = / — g—H, 

suivant les variables elliptiques 

{91 a. e , f. X, g -h 6, 0; 

3 '* Suivant les variables canoniques 
(10 ) L, p z *, X, to/ ; 

4° Suivant les variables canoniques 
(1 * j L, £/. 

et nous aurons principalement a rechercher comment on peut 
passer d 7 un de ces develop pements a un autre. 

! 2 I 4 . La fonction perturbatrice va done se developper sous la 
forme 

(12) B cos ( k\ Xj — f" k %\% ) — 0 sin ( k\ H - A* 2 X2 ), 

si Ton adopte Fun des systernes (9), (10) et (1 1), k\ et k 2 etant des 
entiers, \ { et ies longitudes moyennes. Quant aux coefficients B 
et C, ce sont des fonctions des autres variables, e’est-a-dire des 
elements osculateurs a, e, /, g 6, Q des deux planetes, ou bien 
des L, des p et des co, ou bien encore des L, des q et des y\. 

Si Ton adoptaitles variables (9), le developpement se presente- 
rait sous la forme 

(13) ^ B' cos ( ki Ut *+■ k^u*) C' sin(£i u\ -+* k^u^), 

ou u K et u > sont les anomalies moyennes et ou JB' et G' sont des 
fonctions des 6l6ments 
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Le probleme consiste a determiner ces fonctions B, C, B', O; 
mais on peut 1’envisager de bien des manieres : 

i° On peut chercher a developper ces fonctions suivant les 
puissances des e et des z\ si Ton prend les variables (9); suivant 
celles des p si Ton prend les variables (10); suivant celles des £ et 
des Tj si Ton prend les variables ( 1 1) et envisager separement les 
diilerents termes du developpement. 

On a alors l’avantage d’obtenir des formules analytiques valables 
une fois pour toutes et que Ton peut appliquer ensuite a un cas 
particular quelconque, pourvu que les excentricites et les incli- 
naisons soient assez petites. 

La question qui se pose est alors celle des conditions de con- 
vergence de ces developpements et c’est une de celles que nous 
aurons a examiner. 

2 0 On peut se proposer seulement de determiner, pour deux 
planetes particulieres, les valeurs numeriques de ces fonctions et 
de quelques-unes de leurs derivees. 

Le nombre des termes a calculer se trouve ainsi consid^rable- 
ment diminue, de sorte que le travail est notablement allege, sur- 
tout si Ton ne doit pas pousser trop loin l’approximation. 

A ce point de vue, la question la plus importante sera de deter- 
miner une limite superieure de chaque coefficient, afm de laisser 
de cote les termes dont le coefficient serait certainement trop 
petit. 

3 ° On peut enfin etudier les' proprietes analytiques de ces fonc- 
tions, en vue precisement de faciliter le calcul numdrique dont il 
vient d’etre question. 

Nous verrons que ces fonctions satisfonta des Equations diff£ren- 
tielles lin^aires, dont les coefficients sont des fonctions rationnelles 
des L, £ et 7i, ou encore des 

i g 6 

a , e , tang-, tang—.? tang-- 

Nous verrons ^galement qu’il y a entre ces fonctions de remar- 
quables relations de recurrence. 

215 . J’ai dit que ce qui nous arreterait le plus longtemps, c’est 



LE PROBLEME DE LA FUNCTION PERTURB ATRICE. 


y 


le developpement de la partie de la fonction pert urba trice qui es>t 
de la forme 

75’ 

oil 

D = t x\ — x\ (x % — x’- — (x z — x ?i f. 

Mais nous serons amenes a de\elopper non seulement -4=> mais 

! 1 V 'D 

encore 

i i 

7^ 7 B*’ 

Voici pourqaoi; je suppose que les coordonnees des planetes et, 
par consequent, D, dependent d une quantile tres petite a (par 
exemple l’excentricite i, et que nous voulions de\elopper suivant 
les puissances de cette quantite a. 

Soil D 0 ce que devient D pour a = o, et posons 

D = Du-h s. 

Nous allons d’abord chercher a developper suivant les puissances 
de e et il sera aise ensuite de passer au de\ eloppement suivant les 
puissances de a; on trouve ainsi 

1 I £ I } £- I 

On voit que dans le second membre figure nt non seulement ? 

V Do 

i i 

mats — = 5 -= 

i/D* v/Dl 

Ce n’est pas tout. Dans les equations canoniques (4 bis) et 
(5 bis) du n° 93, par exemple, figurent, non pas la fonction pertur- 
batrice elle-m£me, mais les deidvees partielles de cette fonction. 

Or, si la fonction perturbatrice contient un terme en -d=? la diffe- 
rentiation introduira dans ses derives un terme en —=* 

/D* 

Supposons maintenant qu’on veuille pousser au dela de la 
deuxieme approximation; les seconds membres de nos equations 
prendront, comme nous Tavons vu au n° 100 (t. I, p. 122 ), la 
forme 


^B31V. 
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oil ies B seront des derivees partielles de la fonction perturbatrice, 
qui cette fois ne seront plus du premier ordre, mais d'ordre quel- 


conque, ce qui amenera, non plus seulement un terme ~^==’ mais 


des termes en - 7 =? • * * * 

v/Di 

Voila done deux raisons differentes qui nous obligent a etudier 
le developpement de — !nr> en meme temps que celui de — L* 

™ v/D 3 /D 3 v/d 


216. Reprenons les formules ( 12 ) et considerons le developpe- 
ment de B ou de G >uivant les puissances des excentricites et des 
inclinaisons. Soit m le degre de Tun quelconque des termes de ce 
developpement. 

Nous avons vu au Tome 1 que m est au moins egal a | k\ — k 2 | 
et de meme parite cjue | k A — k>> j. 

Or, si les excentricites et les inclinaisons sont petites, un terme 
est d’autant plus petit que son degre est plus eleve. On aura done 
une bonne approximation de B ou de C en se bornant aux termes 
dont le degre est precisement j k { — ko\. L’ensemble des termes 
constituera ce que nous pourrons appeler la partie prmcipale 
du coefficient B ou C. II sera interessant de chercher ce que devient 
la fonction perturbatrice quand on reduit chacun des coefficients 
de la formule ( 12 ) a sa partie principale. 

Ce que nous venons de dire s’applique encore si, au lieu de 
developper suivant les puissances des excentricites et des inclinai- 
sons, on developpe suivant celles des £ et des r \ . 


217. II arrive quelquefois qu’on est oblige, dans la fonction per- 
turbatrice, de calculer un terme de degre elev£, sans avoir a cal- 
culer les termes de degr^ inf&rieur. En g^n^ral, les coefficients des 
divers termes vont en ddcroissant tres rapidement a mesure que le 
degr£ s’dleve; mais il peut arriver qu’un terme dont le coefficient 
est tres petit ait n^anmoins une grande importance, parce que Tin- 
tegration introduit des petits diviseurs grace auxquels un petit 
terme produit parfois une perturbation considerable. 

Supposons que nous ayons, dans la fonction perturbatrice, un 
terme 

B (303(^1 -f- ki A2), 
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et que les entiers k { et k 2 soient tres grands, mais de telle fagon 
que le rapport — ^ soit voisin du rapport des moyens mouve- 
ments-^* Alors la difference j k { — k 2 j et, par consequent, le 

degre du terme, seront tres eleves; le coefficient B sera, par con- 
sequent, tres petit, mais, en revanche, le petit diviseur k { n K -j- k 2 n 2 
sera egalement tres petit, de sorte que, finalement, la perturbation 
pourra etre tres notable. 

Le calcul exact du coefficient B serait alors tres long, parce qu’il 
exigerait le calcul prealable des termes precedents qui ne sent pas 
directement utiles. On peut Peviter, en se servant des formules 
approchees applicabies seulement aux termes de degre eleve et 
fondees sur les proprietes des fonctions de tres grands nombres, 
soit que ces formules donnent d’emblde une approximation suffi- 
sante, soit qu ? elles permettent seulement de reconnaitre si le terme 
en question est sensible, et si, par consequent, on doit prendre la 
peine de le calculer exactement. 

218. Nous aurons a examiner specialement differents cas parti- 
culiers, dontles principaux sont : 

i° Celui ou les excentricites et les inclinaisons sont nulles; 

2 ° Celui ou les excentricites seulement sont nulles ; dans ces 
deux cas, la distinction entre Panomalie moyenne et l’anomalie 
excentrique et, par consequent, entre les deux ddveloppements ( 12 ) 
et (i3), n’a plus de raison d’etre; 

3° Celui ou les inclinaisons sont nulles. 

Un autre cas particulier qui meritera notre attention est celui de 
la Lune; le rapport des grands axes etant alors tres petit, on a avan- 
tage a developper suivant les puissances de ce rapport, et il en 
resulte une forme particuliere du developpement. 

J’ajoute que les theories de la Lune etant tres nombreuses, les 
formes de developpement qui ont ete proposees dans le cas de cet 
astre le sont Egalement, et nous aurons sans doute Poccasion d J en 
dire quelques mots. 


219. Nousavons dejaeu Poccasion, outre le developpement ( 12 ) 
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qui procede suivant les anomalies moyennes, d’emisager le deve- 
loppement (i3) qui procede suivant les anomalies excentriques. 

Ce dernier developpement a ete employe par Hansen, qui s’est 
aussi servi de developpements procedant suivant Fanomalie 
excentrique de Tune des planetes et Fanomalie moyenne de 
F autre. 

D 5 autre part, Gylden a employe des developpements procedant 
sui\ant les anomalies vraies. 

Quand onprend, comme Font fail ces astronomes, pour variable 
independante Fanomalie, soit vraie, soit excentrique, de Fune des 
planetes, on est oblige a certaines transformations souvent assez 
longues et compliquees. 

En effet, il y a entre les anomalies vraie et excentrique d’une 
planete et Fanomalie moyenne de cette meme planete des relations 
assez simples et bien connues. D’autre part, les anomalies 
moyennes des deux planetes sont liees par une relation lineaire. 

Au contraire, la relation qui lie les deux anomalies excentriques 
(ou bien encore celle qui lie les deux anomalies vraies) est com- 
parativement compliquee. 

Si done on a developpe la fonction perturbatrice suivant les 
deux anomalies excentriques (ou vraies) et qu’on veuille ensuite 
tout exprimer a Faide de la variable independante qui sera Fano- 
malie excentrique (ou vraie) de Fune des deux planetes, il faut 
remplacer, dans le developpement, la seconde anomalie excen- 
trique (ou vraie) par son expression en fonction de la premiere, 
et cela donne lieu a des difficultes de calcul dont nous dirons 
quelques mots. 

220. Dans les equations canoniques, ce n’est pas la fonction 
perturbatrice qui figure directement, mais ses derivees partielles 
du premier ordre. C 7 est done le developpement de ces derivees 
qui serait surtout interessant. 

Quand le developpement de la fonction perturbatrice est donne 
sous la forme analytique, et surtout sous la forme d’une serie pro- 
cedant suivant les puissances des excentricites et des inclinaisons, 
on suivant celles des £ et des yj, on passe immediatement d’un 
developpement a Fautre. 

Mais il n’en est pas de meme quand on s’est borne k determiner 
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la valeur numerique des coefficients, comme nous Favons explique 
au n ,J 214; il faat alors recommencer le travail pour chacune des 
deri\ ees. 

C'est ce qui a amene certains astronomes a preferer aux equa- 
tions canoniques on meme a celles de Lagrange une autre forme 
d'equations. Ils ont remarque en effet que ces derivees partielles 
sont au nombre de six ; on peut au contra ire former des equations 
dans les seconds membres desquelles figurent, non les six derivees 
partielles de la fonction perturbatrice, mais les trots composantes 
de la force perturbatrice, soil suivant les trois axes de coordonnees, 
soit suivant le rayon vecteur et deux perpendiculaires a ce rayon, 
Pune dans le plan de Forbite, l'autre perpendiculaire a ce plan. 
On n'a plus alors a faire que trois developpements au lieu de six. 

C’est qu'en effet les six derivees partielles ne sont pas indepen- 
danles; il y a, entre elles et la fonction perturbatrice elles-meme, 
certaines relations, il y en a encore entre ces derivees et les trois 
composantes de la force, prises de l’une des deux manieres que je 
viens de dire. 

On peut done se proposer de former ces relations et de montrer 
comment on peut s’en servir pour deduire tous les developpe- 
ments de l’un d'entre eux, par exemple de celui de la fonction 
perturbatrice. 

Cette revue rapide montre sous combien de faces differentes 
peut se presenter le probleme du developpement de la fonction 
perturbatrice qui va faire Fobjet des Cbapitres suivants. 
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221. La premiere chose a faire est d’exprimer les coordonn6es 
rectangulaires ou polaires des pianetes en fonctions de Pun quel- 
conque des systemes d’elements osculateurs canoniques ou ellap— 
tiques, definis aux n os o9 et213. 

Nous avons vu, au n° 60 et aux numeros suivants, que les coor- 
donnees sont developpables suivantles puissances des £ et des 7j et 
nous avons etudie quelques-unes des proprieties de ces develop- 
pements. Mais il y a lieu de pousser cette etude plus loin et nous 
commencerons par former les developpements des coordonn^es en 
fonctions des elements elliptiques 

(0 *> e, h ^-+-0, 

Nous supposerons meme d’abord que Ton a pris pour plan 
des x A x. 2 le plan de Porbite, pour axe des x { le gi'and axe; de 
telle fagon que Pinclinaison i soit nulle, de meme que la longitude 
du perihelie g -j- 0, et que la longitude moyenne \ soit egalc 
a Panomalie moyenne l. Dans ces conditions nos coordonnees ne 
dependent plus de la longitude du noeud G, de sorte que nous 
n’avons plus qu’a exprimer ces coordonnees en fonctions de 

^ L — X. 

Dans ces conditions on a en introduisant Panomalie excen- 
trique u : 

^! = a(cosw — e), a? 2 =#siiiw/f — e *, x 3 — o, 

( 2 ) l — u — e sin u. 

II s’agit maintenant d’exprimer x, et x 2 en fonctions de mais 
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pour cela nous a\ons besoin de rappeler d'abord les proprietes des 
fonctions de Bessel dont nous aliens etre amends a faire usage. 


222. Considerons i'expression 

F = E' v * in w . 

oil j'ai represente par E la limite de ( i -4- ^ J pour m = oc, afin 

d'eviter toute confusion avec e qui represente bexcentricite. 

Cette fonction F est une fonction periodique de u , developpabie 
d apres la methode de Fourier en serie proeedant suivant les 
exponentielles 


ou m est un entier positif ou negatif. Les coefficients sont des 
functions de x , de sorte que je puis ecrire 

CS) £i.»*sin«- = y' 


Les fonctions J^cr) s’appellent les fonctions de Bessel. Le 
premier membre de (3) est le produit des deux exponentielles 

-20)’^ 

d J ou, par multiplication, 

Pour avoir J,,,, il faut chercher dans le second membre les 
termes en E' m “, e’est-a-dire faire a = {3 -+- m. Nous trouvons ainsi 

(— i)P /af\»t+*p 


<4) 


j »(*> =2sf 




P ! (3 -+- m ! 

Si m est positif ou nul, il faut donner a (3 les valeurs 

O, I, 2 , ... ? 

et appliquer la formule en supposant 


o ! = i ! = i. 
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On voit alors que . est line fonction entiere de .r, divisible 

par x m : la serie (4 ) converge en effet tres rapidement pour toutes 
les valeurs de x\ 

Si m est negatif, il faut donner a 3 les valeurs 
— m, — m — i . — m •>, 

de facon que a soil positif ou nuL Alors S m ( x ) est encore une 
fonction entiere de x divisible cette fois par x~ m . 

Observons d'ailleurs que le developpement (4) ne contient que 
des termes de degre pair si ni est pair, de degre impair si m est 
impair. On a done 

( 5 ) 1 m ( X ) = ( I ) m J f /i [ X ). 

Si d’autre part, dans la formule (3), nous changeons x en — x 
et u en — u , elle devient 

(3 bis) B* r sin “ = ^ J m (— x ) 

En identifiant les deux developpements (3) et (3 bis) de E za ‘ sinw , 
on trouve 

J —m(x) = J/«( — X)i 

d’ou 

(6) J-/«= (— 

formule qui permet de ne considerer que des fonctions de Bessel 
d’indice positif. 


223. DifFerentions l’equation (3) par rapport a u , il vient 

ix cos w E* r sin " = ^ im 3 /n E * rnu 
ou 

~ ( E ,w -h E~ iu ) =2 m J m E imu, 

En identifiant les coefficients de E imu dans les deux membres il 
vient 

(7) ^ (Jm-i-H- J^-4-i) = ^mJ m , 

relation de recurrence entre trois fonctions de Bessel consdcutives. 
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Differentions maintenant l’equation ^3) par rapport a x , il 
viendra 


i sin«E z * s,n " = ^P,T,„| x !& mu 


(E lu — J/« E imu = ‘2 yV m 

ou en identifiant les deux membres 

( 9 ) ^ ^ in ~ J m— I J m-r-l j 

formule qui permet de caiculer la derivee de J m . 

224. Differentions maintenant la formule (3 j deux fois par rap- 
port a u y ou bien encore deux fois par rapport a x , nous obtien- 
drons ainsi les deux formules 

( 10 ) — ix sin^E** sJn * r sin « = — ^ m 2 J m , 

(n) — sin*ME'-* : » ta "=2j5*^ ar ) E,w “- 


Ajoutons maintenant les quatre formules (3), (8), ( 10 ), ( 11 ) 
apres les avoir respectivement multiplies par 

x 2 , -h x, 1, x- ; 

il viendra 

xy m -*r(x 2 — m 2 )l m ]E imu = o; 

d’ou 

(12) -h (x 2 — m?) J m = o, 

equation diff&rentielle lin^aire du second ordre a laquelle satisfait 
la fonction de Bessel 


225. Bien que la serie(4) converge tres rapidement pour toutes 
les valeurs de x , il peut j avoir iniret a recbercher une valeur 
approclie de la fonction J m pour x tres grand. Nous avons par la 
formule de Fourier 

/ 27T 

P. — II. 2 
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ou ea posant 

(\i) 


sin // = 3 , 


71 J, 


I E~ 


dz 


/r 


Supposons x reel, positif et tres grand. \u lieu de faire Pinte- 
gration le long de la droite qui joint le point r* = — i an 
point s = i, c'est-a-dire en donnant a ^ des valeurs reelles, nous 
pouvons la faire le long d'un autre chernin ayant raemes extre- 
mites. 

Nous choisirons un chernin tel que la partie imaginaire de ; 
soit constamment positive, sauf bien entendu aux deux extre- 
mites 5=±i. Dans ces conditions ixz a sa partie reelle nega- 
tive et tres grande, et E lxz est tres petit. Les seules parties du 
chernin d’integration qui pourront donner des termes sensibles 
sont done les parties voisines des deux extremites, parce qu'aux 
extremites, z devient reel et que la partie reelle de ixz devient 
nulle. 

Mais pres de Pextremite z = i : nous avons sensiblement 

7 C 

— ~IW— 

w = - , E~^|«=E v /i — z 1 — \J‘i ( l — Z). 

*2 


Pres de Pextremite £ — — i ; nous avons sensiblement 




im 

E-itnii — E 


7 T 


/l — 2 2 =/'dl + 3 ). 


Mais il importe de voir quel signe il convient de donner aux 
radicaux. Nous devons supposer que le signe de yj i — z- a ete 
choisi de fagon que, pour z reel et compris entre — i et -j~ i , le 
radical soit reel et positif. Nous pouvons supposer cPautre part 
que le chernin cP integration aboutit a ses deux extrdmit^s ± i par 
deux petits segments de droite de longueur e, perpendiculaires 
a Paxe des quantites reelles. Ge sera d’ailleurs les seules parties de 
ce chernin que nous conserverons. 

L’argument de 1 — sera a l° rs ^ l e long du segment aboutis- 

sant k z = — i, et — ^ le long du segment aboutissant a z = x. 


Le longdu premier segment nous pourrons prendre argy/sH-i = - ^ 
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et le long du second argv's — 1 = nous aurons done 

/ 1 — ^ — «V‘ 4 — 1 ; /i — z-— v ' H- 

d’ou 

4) “j/«= 


< I 


//r 

E'«E '-~= 
y-i v/ 1 — 




E 


spi\J z — \ 


Mais nous pouvons rein placer Ies limites superieures dz 1 -f- sf 
des deux integrates par Finfini ; car nous ajoutons ainsi aux che- 
mins d integration des lignes le long desquelles la partie imagi- 
naire de z est positive et notable de sorte que E ixz est negligeable. 
Or 


Done 

d’oii 

on enfin 

U?) 


+ 

E i**dz /4 

/ = . = E lx/l 4 / -E 4 , 

Jh \/z — h V x 

T a / '* l miZ ' K \ 

J ' rt= \/ — cos(*- — --J. 


226. Vojons maintenant comment les transcendantes de Bessel 
peuvent etre employees an developpement des coordonnees du 
mouvement elliptique. Considerons l’exponentieile 

ou p est un entier; e’est une fonction periodique de u et par con- 
sequent aussi une fonction periodique de l. Elle est done develop- 
pable en s6rie de Fourier sous la forme 

il s’agit de calculer le coefficient A OT ; la formule de Fourier nous 
donne 

2izA m — I EP iu E~ mi/ dL 
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Or 

Q-nui dl = — dE~ mil \ dQP lu = ipEP^du. 
m 

On trouve done en integrant par parties 

2T*A. m =£- f EP lu E- imI du, 

171 J 0 


ou a cause de l’equation de Kepler 

27 T A, n = — / E^P- m ^ u E £mesitlu du. 

mj 

Mais l’integrale e’est evidemment au facteur 2 7c pi'&s le coefficient 
de E im ~P ]iu dans le deveJoppement de E imesinu. e jj e est (j onc egale a 

2.715 m —j, (me), * 

d’ou la formule cliercliee 

( 16 ) A m = — me ). 


227. Cette demonstration ne s’applique plus pour m — o et 
d’ailleurs la formule devient illusoire. On a alors 


27C An 


EP^dl 


OU 


ou 


-/ 

2tcA 0 = J Ep*' m (i — e cos u)du 


2TC A ( 


-/ 






ce qui montre que A 0 est <%al a i pour p = o,k — ^ pour p =zt i r 
et a o pour toutes les autres valeurs de p. 


228. Proposons-nous maintenant de d^velopper suivant les. 
exponentielles E‘ ml une fonction p^riodique quelconque de u 
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et de la mettre ainsi sous la forme 

Fuo=yA„ ? E iml . 

Les principes des deux numeros precedents nous donneront 

(17) — ~ ^m-p ( h 

pour m^o et 

(18) A 0 = B ^ 
pour m = o. 

Inutile cPajouter qu’en partant de Pidentite 
E i r m E ic i u = 

en y remplacant les trois exponentielles par leurs de\eloppements 
et en identifiant les deux developpements, on arrive a un grand 
nombre de relations ou le premier membre est une fonction de 
Bessel et le second membre une serie dont tous les termes sont 
des produits de fonctions de Bessel. 

229 . Grace aux formules ('17} et (18) du numero precedent, le 
developpement d'une fonction quelconque de u suivant les puis- 
sances de se deduit aisement du ddveloppement suivant les 
puissances de E*“; or la plupart des fonctions interessantes rela- 
tives au mouvement elliptique s’expriment tres simplement al’aide 
de E iu . 

Supposons d’abord que nous prenions pour plan des x { le 
plan de Porbite et pour axe des x K le grand axe, nous aurons 

x t = a{ cos u — e ) — — ae -4- — E /w -+- — 

22 

a \J 1 — e 2 sin u = ^ \f 1 — e 1 ( E fw — E ~ £u ), 

x z = o, 

r = \Jx f -b a?§ = a ( 1 — e cos u ). 

IPappli cation des formules (17) et(t8) nous donne 
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A f)i — - [ J m — i ( JY16 ) J in~\r 1 ( XL6 )J ; 


Kn = (me) -h J,n+i( me)]. 


a; = o. 


Les formules de recurrence (7) et (9) nous permettent d’ecrire 


A m — “ J m ( me )> 


A'„= ■ 6 - J „(me). 

lem 


Nous trouverions de meme 


. y jklggj e *«. 

Jmd 771 


Sous le signe l’indice m prend toutes les valeurs entieres 

positives et negatives, a Fexception de la valeur o. Cette formule 
se deduit imm^diatement de la relation 


r = a(i — e- ) — ex t . 


Nous pourrions ecrire egalement 


jsx^iast E /tt , / -\ E / / -n 

+ _ ^ , — h v/ 1 e l ) - — (1 


#1 £# 2 


E-za . £*« / , N 

— — (1+/1- «*)-+- — (x— /i — e 2 ), 


d ? ou Pon deduirait par les formules (17) et (18) les develop- 
pements de .r* dr ix 2 . 

230 . Reportons-nous maintenant au n° 63 et aux quantity que 
nous j avons appel^es X et Y ; nous voyons que Ton aura 

X — h* i Y = ( x \ — )— i x 2 ) E~ z/ , 

et par consequent, par application des formules (17) et (18), 


X-+- iY 


(21) 


'-E- ti -h -([-+■ y E «,„_ 

2 x v y 2/n 


2 /n 
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On en deduirait le develop pement de X — i Y en ehangeant i 
en — L et par consequent ies developpements de X et de Y. 

Supposons inaintenant que nous rapportions notre svsteme 
a des axes queleonques : les quantiles anxiliaires X et Y seront 
toujours deludes coinrae an 11 “ 63 el leur deveioppement en func- 
tion des variables a. e et l ne change ra pas. Quant aux eoordon- 
nees rectangulaires rapportees aux nouveaux axes, elles seront 
bees a X et \ par les equations ( 12 ) du n“ 63, t. 1, p. 79 . 

Ces formules pement aussi s'errire 


( 9 / 2 ) 


^ x 1 — ix-i = cos- i I X -J- i V ) E z * sin 2 ^ < X — i Y ) E 
| ;r^ — partie imag. <le 1 ^ -r- i 


Lai remplace la lettre i des formules du n° 63, qui reprd- 
sentait I'inclinaison, par la lettre J 7 alin d'eviter toute confusion 
avec i = \J — 1 ; on obtiendrait Texpression de x K — ix 2 en chan- 
geant i en — i. 

En rapprochanl les formules ( 21 ) et ( 22 ) on trouve 


X\ ~ — ix. 1 — a cos 


(*23) 


ns* £ if E‘V+0 

sin 4 i r~ — 

2 L -»• 

^2d 4^ (Mm-t-Mm-H )J- 


J'ai ecrit pour abreger £ t et s 2 an lieu de 1 — E- y/ 1 — e 1 et 
de 1 — \j\ — e-. L’argument des fonctions de Bessel i m ~\ et J m+f 
est egai a me. 

On trouve de meme 


(23 bis) x% — 


a sin J 


[-H 


S1U g 


sin ( ml - 


4 m 


“ ” ( £ l J m— 1 — s 2 J )1 


231. Les developpements obtenus par la methode du n° 228 ne 
sont pas les seuls que l’on puisse imaginer. Supposons par 
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exemple que Ton veuille developper 

du __ i 
~dl i — e cos u 

On pourraii evidemment developper ' t ” en sene de Fourier 

et appliquer ensuite la methode du n° 228; mais il j a quelque 
chose de beaucoup plus simple. 

INous avons 

u = / 4- e sin w, 

ou, en se reportant an developpement de au n ° 229 et a la for- 
mule ( 20 ), 

K==/H -8fc2 J " i(we)E "“ 7> 

ou, en differentiant, 

04 ) § = . + 2 J " i(7 " e)E "“ 7 ' 

Ce developpement nous donne en meme temps celui de 

1 _ r 

7* a(i — e cos w) 

Plus generalement, cela nous donne le mojen de developper les 
expressions de la forme 

F (u) 

1 — e cos a 

Si, en effet, Ton a 

F(k)=2 C '' E,> “’ 

on aura 

f(«) 

t — e cos u dl iv 

On d^veloppera alors V -jf & ipu en s ^ r * e p^ocedant suivant les 

exponentielles E /m/ et Ton differentiera par rapport a /. 

Cette methode peut etre plus avantageuse que celle du n° 228 
si le developpement de F(w) suivant les exponentielles E l ? u est 

plus simple que celui de — — 7 si par exemple il se reduit 
1 r u 1 — e cos a r r 
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a un nombre fini de lermes. £n particulier nous tromerons les 

j , i cos u sin u . . , 

dev eloppements de — — ? — - - ? expressions qui sont bees a cose 

et a sine par des relations lineaires (e est Panomalie vraie). 

232. Nous pourrions aussi diflerentier par rapport a e\ il \ient 

, chi 

(i — ecosn) -r- = sinw. 
de 


Du developpement d'une fonetion periodique quelconque <£ ( u ) 
nous deduirons done immediatement celui de 


d< P 
du 


i — e cos u 


ee qui pourra etre utile si le developpement de ^ est plus simple 

que celui de si nw~* En particulier on peut calculer de cette 
fagon 

sinw 

r 


qui ne differe de sinc> que par un faeteur constant. 

Mais nous pouvons egalement diflerentier deux ou plusieurs 
fois par rapport a /, ou au contraire integrer par rapport a L 

Par exemple les equations du mouvement keplerien nous 
donnent 


d t x l _____ a z Xt 

dfi r 


Or nous connaissons les developpements des coordonnees 
ou le coefficient de chaque terme ne depend que d’un nombre fini 
de fonctions de Bessel; Pdquation ( 20 ) nous donnera done celui 

de chaque terme ne dependant que d’un nombre fini de fonc- 
tions de Bessel. En particulier, si nous prenons les axes du n°229, 
nous connaitrons les ddveloppements de 

x x cos 9 x % _ sinp 

2 y 3 


v 6tant Panomalie vraie. 
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Si nous multiplions chacune des equations (20) par xi et que 
nous ajoutions, nous trouverons 



Comme le second terme du premier membre de (26) est une 
fonction lineaire de j en vertu de Fequation des forces vives, 

^equation (26) nous apprend qu’il y a une relation lineaire 

1 cl^r ' 1 
entre - et -nr- 
r dl l 

Or Tequation (24) nous donne le developpement de nous 

d 2 r~ . 

possedons done celui de d’ou, par une double integration, 

celui de r 2 . Nous avons done les developpements de r, ^ et r 2 . 
L’equation de Fellipse en coordonnees polaires nous donne 

a(i — e 2 ) 

e cos p = — 1, 

r 

er 2 cos p = a ( 1 — e-)r r 2 . 

Des developpements de r et r ' 2 nous deduirons done ceux 

de cos v et r-cosp. Nous avions d’ailleurs explique plus haut la 
facon de developper cosp. 

Tous les developpements dont il vient d’etre question sont tels 
que chaque terme ne depend que d’un nombre fini des fonctions 
de Bessel. Inutile d’ajouter que, par l’emploi des relations (7), (9) 
et (12), on peut s’arranger pour que chaque terme ne contienne 
plus que 

J m(me) et V m {me). 

233 . Nous pourrions nous proposer de developper les coor- 
donn^es en fonction des Aments canoniques 

L, p, to, 7). 

Si nous nous reportons aux formules ( 23 ) et (23 bis) nous 
verrons qu’il suffit pour cela de developper les exponentielles 
les expressions 

a cos 2 ^ 1-0) ) tang 2 - tang^e*^, 


(27) 
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et enfin la fonction 


K 


m 


-1 Tw— I =>i J /n-r-i 


*7 


Le< exponentielle^ \L mu - sont directement exprimees par ies ele- 
ments canoniques. Ies expressions (*27) se reduisent respecti- 
\ ement a 


L* 



i;i± if, 1 i£ 2 ±:zr l± ) 2 




en faisant abstraction <le facteurs constants. 

Quant a R //; il est developpable suivant les puissances crois- 
sanles de 



et le developpetnenl se deduit aisement de eelui des fonc lions de 
Bessel. Je n insisterai pas da\antage sur ces points, me bornant 
a renvojer aux n 0> 60 a 69 du Tome T 1 . 


23 i. INous observerons maintenanl que si e est tres petit Ie pre- 
mier terme de ebaque fonction de Bessel sera plus considerable 
que tous les autres. Nous pourrons nous borner alors a ce terme 
si nous nous placons au point de vue du n° 216. En partant des 
formules (id’) et remplacant chaque fonction de Bessel par son 
developpement, nous trouvons 


(**8> 



P i—ift 

m 3 ! m — p -r- 


i me y«-AHrsp 

pi \T / 




et en reduisant chaque fonction de Bessel a son premier terme 


(»9) 


/ me \ 

m -t> w ml , p 

( me \ 

p ~ m (~ 1 )p~ m 

\~) 

m — p ! ‘ jmd m 

\ a / 

p — ml 


E On/ m 


Sous le premier signe ^ on doit donner a m toutes les valeurs 
entieres telles que m — p soit positif ou nul, et sous le second 
signe ^ toutes les valeurs telles que m — p soit n^gatif. 

Mais nous devons observer que l’approximation relative avec 
laquelle chaque fonction de Bessel est ainsi represent^ est d’au- 
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tant plus faible que l’ordre de cette fonction est plus eleve. Et en 
eflet le rapport du second terme au premier s'ecrit ( m — p etant 
positif, par exemple) 



11 est tres petit si e est tres petit et m fini ; mais, quel que soil e, 
il eroit indefiniment a\ec m. Pour les termes d’ordre un pen eleve, 
on obtiendrait une bien meilleure approximation en remplacant les 
fonctions de Bessel par leur valeur approchee du n° 225. 

235. Le developpement (28) et les developpements analogues 
procedent suivant les puissances de e et les exponentielles E iml , Si 
nous les ordonnons d’abord suivant les exponentielles E imi , le 
coefficient de chaque terme, etant une fonction de Bessel, pourra 
etre developpe suivant les puissances de e en une serie toujours 
convergente; de plus, des qu’on aura calculi les coefficients des 
differents termes de la serie de Fourier, cette serie elle-meme sera 
toujoui's convergente. 

Avec cette fa^on d’ordonner les termes, la convergence est done 
assuree, mais il n’en est plus de meme si Pon ordonne d’abord 
suivant les puissances de e, le coefficient de chaque terme 6tant 
lui -meme une fonction de L Nous sommes done amenes a nous 
poser la question suivante : 

Nous avons une fonction quelconque de w, dependant par con- 
sequent de e et de /, nous la developpons selon les puissances de e\ 
quel est le rayon de convergence de la serie, e’est-a-dire la valeur 
maximum de | e | pour laquelle la serie converge? Cette valeur 
depend evidemment de l et je puis la designer par <p(7). Si alors M 
est la plus petite valeur de cp( /), quand l prend toutes les valeurs 
reelles possibles, la convergence de la serie sera assuree tant quee 
ne surpassera pas M. 

Pour calculer f(l) et M, nous allons appliquer le th£or&me de 
Cauchy. Pour que la s^rie cesse d’etre convergente, il faut et il 
suffit que la fonction cesse d’etre holomorphe. Or nous avons 

l — u — e sin w, 

d’ou 

du( i — e cos u) = dl h- sin u de. 
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ce qui montre que u (et en general toute fonetion de u) cesse 
d'etre une fonetion uniforme de e et de / pour 

( 3o ) e — — 

eos u 

E 11 combinant ^equation (3o) avee celle de Kepler, on trouve 
{ 3 1 ) u — tangtt = l. 

Les equations (3o) et (3i) definissent e et sa valeur absolue | e\ 
en fonetion de l et la plus petite des determinations de cette 
valeur absolue est o(l)- 

Remarquons que <p(/) est le meme en general pour toute fonc- 
tion F(&); il n’j aurait d’exception que si la derivee s’annulait 

precisement pour la valeur de u qui satisfait a l’equation ( 3 i). 
Cela nous permettra de raisonner sur une fonetion F(&) quel- 
conque, il nous suffira que la condition d’exception ne soil pas 
remplie. 


236. Il s’agit maintenant de determiner la valeur de l pour 
laquelle <p(/) atteint son minimum M. A cet effet je reprends la 
formule ( 28 ) et j’obtiens en la difierentiant par rapport a l 


_ 2sr 


(-i)P 

film — p fi ! 


m~ p-+-2^ 


E iml— <*>(/). 


Cette fonetion ne se trouve pas dans le cas d’exception signale 
a la fin du num^ro precedent, car, quand on a 


1 — e cos u = o, 


sa d^riv^e ne s’annule pas et devient au contraire infinie. 

Les Equations (3o) et (3i) ne changent pas quand on y change 
/ et e en u + tt, L -f- it, — e; il suit de la que les valeurs cri- 
tiques de e qui correspondent a / et a L -f- tc ont meme valeur 
absolue ©(/) — <d(1 + tz) 7 mais sont egales et de signe contraire. 

Supposons done que pour une valeur r^elle l K de l et pour une 
valeur quelconque de e que j’appelle'e* la serie (32) diverge, il 
en sera de m£me quand nous changerons l K en tc, et il en sera 
encore de m&me de la somme 
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Considerons, en effet, <!>(/') comine uae function de e el de / et 
faisons-y /— ! h : cette fonction presentera un point singulier pour 

e = e\ \e' | = ci l x i < e t 

et n’en aura pas pour e = — e ! : car, si Ton change e en — e dans les 
equations ( 3o \ et i3i ), / se change en t: ± / qui est different de l 

sauf pour / = ^ • 

La fonction <E> ( / 4 -t— tz ) presentera de meme un point singulier 
pour e = — e et n*en aura pas pour e = e' et la somme 

( / ! > — J— < t > < / 1 -T- 71 I 


presentera deux points singuliers e = e\ e = — et, ces points 
singuliers etant distincts, les sing ularites correspondantes ne pour- 
ront se detruire. Done cette somme divergera pour e = e K . 

Je puis done ecrire 


( 33 } 




oil l*on ne donne a m que des \aleurs paires . 

11 resulte de ce qui precede que la serie *F(^, e { ) diverge. 
Passons de cette serie a la suivante : 


*(- 1 ’ =*(- 


ou Ton remplace Fanomalie moyenne par ^ et Fexcentricite par i 

multiplie par la valeur absolue de e K qui est essentiellement reelle 
et positive. 

Si dans le dernier inembre de Pequation (33) nous faisons cette 
substitution, la valeur absolue de chaque terme demeurera la 
meme; mais tous les lermes vont devenir positifs ou toutau moins 
avoir meme argument (en supposant que le nombre p soit pair, ce 
que nous pouvons faire). En effet, le denominateur est essentiel- 
lement positif, Pargument de ( — i)P est (3 t r; celui de 

/ m \ m-jH- 

W 

est nul si m et p sont pairs comine nous le supposons; celui de 
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sera 

celui de E imi sera 


I m — /> ) - — 3 r. 


— />z 


a 


LTargument resultant sera done 

-j 5 - + 'i 3 tt 


(ou encore — p ~ ■> en supprimant un multiple de 27 tj; il sera done 
constant et independant de m et de p. c. q. f. d. 


Si done la serie diverge, il en sera de meme, a fortiori, 

de W ^ ~ j dont les termes ont meme valeur absolue, mais avec 

un argument constant, au lieu d’avoir un argument variable. 
Done 


di\erge; il doit done en etre ainsi au rnoins de Fune des deux 
series <I> ^ ^ j * ^ ( i) ’ ^ ais nous venons de voir que 

©(/ -h 1Z) s= ©(/), 


e’est-a-dire que <f>(/-|- tc) diverge si <£(/) diverge, et inversement. 
Si done Fune des deux series diverge, elles divergent toutes deux. 
Done 


diverge. Done 


*(r f »“i) 

\e l \ > z(j'y 


Or, tout ce que nous avons suppose sur e K e’est que 

I | > q(7); 

cette seconde inegalite entraine done la prdc^deute, e’est-^-dire 
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que 

Done le minimum de 'p(Z) a lieu pour / = ~ * 

237. Pour determiner M il faut done etudier Fequation 
(34) a — tang a = 

ou en posant ^ — u = s 

( 34 Z >£$ ) S "T" cot £ = o. 

Je me propose de demontrer que ^equation (34 bis ) ne peut 
avoir que des racines purement imaginaires. Posons, en effet, 


il vient 


9 = cosep, 
d 29 


a?p 2 


- £2 9 = 0 

avec [en tenant compte de (34 bis)] 

Dour p = 1 , 
pour p = o. 


(35) 


do _ 
do ~~ 


do 

d~ 9 ~°’ 


Supposons que Pequation (34 bis ) ait deux racines e ets 1 ; a ces 
deux racines correspondront deux fonctions (pet et l’on aura 


do 1 


<*?v 


r 1 o ? 2 9 ! <^ 2 9 \ ■. 

Le second membre s’annule, tant pour p = 0 que pour p = 1 , a 
cause des equations (35), et il reste 

(s 2 — sf) f 99^ = 0. 

d ft 


Si liquation (35 bis) admet une racine e, elle admettra la 
racine imaginaire conjug'uee e, ; les deux fonctions cp et cp ( seront 
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imaginaires conjuguees, leur produit sera essentiellement positif, 
de sorte qu'il restera 

s*— sf = °- 

Cela ne peut avoir lieu que si s 2 est reel. Si £ 2 est reel positif, 
- est reel. Mais les racines reelles ne sanraient convenir a la ques- 
tion, car Inequation (3oj conduirait a une \aleur de e plus grande 
que i. 

Si = 2 est reel negatif, s est purement imaginaire. II s 7 agit done 
de rechercher la plus petite racine purement imaginaire del'equa- 
tion ( 34 bis j ou, ce qui revient au meme, en changeant s en i s, la 
plus petite racine reelle de 

( 36 ) & ( E - — E- j — h ( E~ e — {— E E } = o i 

on en deduira la valeur de M par la formuie 

On trouve ainsi 

M = 0,6627, 

ce qui prouve que nos series convergent toutes les fois que I’ex- 
centricite est inferieure a 0 , 6627 ; l’equation (36) n’a d’ailleurs 
qu’une seule racine positive. 


P, — II. 


3 



CHAPITRE XVI. 

PROPRIETIES GENERALES DE LA FONCTION PERTURBATRIGE. 


238. II s’agit maintenant de passer au ddveloppement de la 
fonction perturbatrice elle-meme et je voudrais exposer d’abord 
quelques considerations generates sur la facon dont se pose ce 
probleme. 

Au n° 212 nous avons rappele sous quelles formes diverses se 
presente cette fonction. Nous devons d’abord nous occuper du 
developpement de l’expression (5) du n° 212 

— m-i m k 

y/( a/* — x\ ) 2 -+- ( x ' s — x'z ) 2 H- ( x' 6 — V 3 ) 2 

c’est-a-dire de ce qu’on appelle la partie principale de la fonction 
perturbatrice; c’esl le probleme le plus difficile, et, une fois qu’il 
sera resolu, le developpement des autres parties de la fonction 
perturbatrice s’en deduira presque immediatement. Supposons, 
en effel, qu’on ait developpe cette expression (5) en fonction des 
elements elliptiques osculateurs de Ja premiere planete 

hi hi H 6i, 0i, 

et de ceux de la seconde 

1 e 2i hi hi ( 1 )i 

et que l’on ait trouve ainsi 


(i) 






( l ) Cependant pour Economiser les indices je designerai quelquefois ces Ele- 
ments par a , a .... 
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ou f est une fonction des 12 elements elliptiques osculateurs. On 
en deduira 


( 2 i 


F— ( x\ — hx j )- 


= f(ha u 


h etant un coefficient constant quelconque. Si, en effet, on change 
a { en ha ,, ies autres elements e t , f,, Z f , «*i et 6, ne changeant pas, 
non plus que Ies elements osculateurs de la seconde planete, les 
coordonnees de la premiere planete, x\ 7 x 21 x Z7 se changeront en 
hx\, hx\ r hx'^ et celles de la seconde planete x \ , x[ y7 x ' e , ne chan- 
geront pas. 

Or, nous avons \aau n u 36 que la fonction perturbatrice, si Ton 
adopte les variables dtt n° 30, prend la forme 


avec 


m 1 m w 


(Jl L 

\ BD AB, 


rrix m 


*(■>■ 


<m)- 


BD 5 = x} -r- x’-- ~h x e - 


I 


~'S 


ab- 2 ^*- 

BC* = 2f^ + 


m- t 

m-i 

m l 

mi -+■ m~ t 



II s’agit d’obtenir les developpements de 


1 1 1 

IT bc’ bd' 


Or, Fapplication de la formule ( 2 ) nous donne immddiatement 

1 j./ 7 n n \ 

AB [m, -+-/», au ■"/’ 

I r! ~ m \ \ 

BC =A^^ a, ’ ai ’ 

g^=/(o, «,, ...)• 


Cette derniere expression /*(o, a 2l . . .) ne peut evidemment 
dependre que des elements de la seconde planete; on Fobtiendra 
imm^diatement si Ton observe qu’elle n’est autre chose que Fex- 

pression ~ relative a la seconde planete, et que nous avons appris 

plus haul au n° 231 a developper 
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Supposons maintenant que nous ayons adopte, au lieu des 
variables du n° 30, celles du n 1J 26 ou celles du ri' 44. II faudra 
alors developper la partie complementaire de la fonction pertur- 
batrice. qui. ainsi que nous 1‘avons rappele au n° 213, peut prendre 
Tune des trois formes suivantes 


"B C*" 






i 

m- 




Reprenons la formule ( 2 ) et developpons les deux membres 
suivant les puissances de h en negligeant les termes en h 2 et les 
termes suivants. En observant que, avec les variables des n os 26 et 
44, on a 

bc*= yv 4 * 

et non plus 

BD’-=yv 4 *, 

il viendra 


1 ^x\x\ 

BG : BC 3 


= /( o, a a , . . .) -t- ha x 



On en deduit 
(3; 


2x'x[ df 

"BG ~ = ai dar 


Dans il faut faire a, = o apres la differentiation. 
On trouverait de m&me 


(3 bis) 


2x[ x\ df 


et ici encore, dans il faudrait faire a 2 = o apres la differentia- 
tion. 


239. Il reste a ddvelopper 

y* ■ 


Il nous suffira, pour rattacher cette expression a la fonction 

/(<*!> «!, 
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de la rattaeher aux expressions 

£ JT , x \ £ x \ 

BC 3 * AC 3 

liees a / par les relations | 3 ) et ( 3 bis). 

A cet efiet, je les relierai les unes et les autres a la fonction 

^ =2^ *+• 


Les equations du mouvement keplerien nous donnent 
d-x\ M.r t 

dr /* 3 ? 

M etant la masse centrale attirante et r le rayon vecteur, ou bien 

, d l x t Mar, 

dr- ~ r 3 ' 


n designant le moyen mouvement, et / 1’anomalie moyenne. 

Nous pouvons appliquer cette formule a notre seconde plan&te 
Active puisque les relations entre les coordonnees de la planete et 
les elements osculateurs sont les memes que dans le mouvement 
keplerien d’apres la definition meme des elements osculateurs. 
Nous devons done changer 

x u r, l 7 7i 2 , M 
en 

x\, BC, U, n|=“, M, 

d’oii 

d- od r 4 _ a | x\ 

dl\ ~~~ 


Nous aurions deux equations que Von d^duirait de la prec^dente 
en ehangeant x\ en x\ ou en x\. Ajoutons ces trois equations 
apr£s les avoir multiplies par x \ , x{ 21 il viendra 




d % od k 

~dif 


2 x\ x[ 
1 BC 3 


ou en se reportant a la formule (3) et en se rappelant que x\,x[ v 
x ! t ne dependent pas de L 2 , mais seulement des diluents oscula- 
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teurs de la premiere planete, 

dr XT 


(4 > 


on trouverait de meme 


dl% 


= — a\a x ——— j 
da { 


(4 bis) 


d*W 


dl\~ a ' a ~da , / 
Occupons-nous maintenant de Fexpression 


df 


Cette expression s’introduit quand on adopte les variables du 
n° 26; dans ce cas les masses m\ et m \ attribuees auxdeux plane tes 
ficlives ont pour valeurs 


m x m 7 
ni\ m 7 


n\\ m- 
m x -h m 7 


(cf. n° 43); on a alors, dans le mouvement keplerien, 

(5) 

ou bien 

(6) y\ = m\ n x 


/ , dx\ dx' w 


dx\ 


dl x 


y\ — m’ k m 


dx ! . 


dl o 


n { et n 2 etant les deux mojens mouvexnents. Les formules (5) ne 
seraient plus vraies dans le mouvement trouble; mais les foi'- 
mules (6) le seront encore, puisqu’elles expriment une relation 
entre les x r , les y r et les elements osculateurs et que ces illations 
sont les memes dans le mouvement trouble et dans le mouvement 
keplerien d’apres la definition meme des elements osculateurs. 

On aura done 


( 7 ) 




m\ m\ n x n 2 


2 


dx\ 

dl x 


dx\ 

dU 


d 2 XJT 


Ainsi y\ se trouve rattach6e a? et par la k f \ 


240. Consid^rons maintenant une fonction p^riodique quel- 
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conque 


F { u . u ) 


des deux anomalies excentriques u et id) elle peut se developper 
en serie de Fourier sous la forme 


( 8 > 




les p et Ies p ! etant des entiers positifs ou negatifs, mais on peut 
aussi la developper en serie de Fourier proeedant survant les ano- 
malies moyennes / et d sous la forme 


( 9 > 




II s’agit de savoir quelles relations il y a entre les coefficients 
des deux series ( 8 ) et(9). Ces coefficients nous sont donnes Tun 
et Pautre par des integrales definies 


(10) 

47 t 5 Bp P -_ f 

et 


(ID 

4 K*A., nm ' = J 1 


dldl' , 


Les integrales doivent etre prises entre les limites o et 2ir, 
tant pour u et u! que pour l et i J . 

Nous transformerons la formuie (i i) par une double integration 
par parties. 

Nous pouvons poser : 

d*<$> 


£— Uml-t-m'/') — 


dl dd ’ 


avec 


— L ff'-Uml-hm'/')' 

mm! 


Nous trouvons alors, en integrant par parties par rapport a u 

ff e £h d,dr —ff d £§ d “ d, \ 

et, en integrant encore par parties par rapport & u 1 et a 

ff e sm d “ lr =ff £h* d “ M ’ 



4o 

il \ient done 
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d* F 
du du 


Q-unl-r-m'}') du dll'. 


Or, du deseloppement de F sous ia forme ( 8 ), nous pourrons 
deduire 


du du' 


= — ^ pp'B PP ’E"P u -rp' u,t , 


d’ou, en tenant compte de l'equation de Kepler, 

X/nm^^pp'Bpjy J jW- du du\ 


ou 


A = (p — m)u -+-(// — m' )u', 
Q = /(me sin w -j- mV sinz^'), 


et ou e et e } designent bien entendu les deux excentricites. 

Le coefficient de pp'^pp' est une integrate double qui peut se 
decomposer en le produit de deux integraies simples 


p 


•/ 


g2(p--;«)M j£//««sum du j EHp'—ni^n' }£im’e' sinn 1 du! 


c ? est-a-dire 

4 TT 2 J m—p ( MIG ) J m r —p' (ffl e )• 

Nous avons done la formule 

(l2) ~37~7 & pp> J m —p ( TY16 ) J rn'—p'jm & )? 

Jmm 11XTYI 

tout a fait analogue a la formule ( 17 ) du Chapitre precedent. Cette 
formule est due a M. Baillaut. 


241. Cette formule se trouve en defaut quand m ou m! sont 
nuls. Mais elle peut dans ce cas etre modifi^e et transform^ en une 
autre analogue k la formule ( 18 ) du Chapitre prudent. 

Nous avons 

F(u , u r ) E Hpu-hp'*'\ 

d’autre part nous avons 
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formule dans laquelle on doit, pour m — o, remplacer le coeffi- 
cient ~ J m _p par i pour p = o, par — ~ pour p = ± i , par o dans 
tous les autres cas. De meme on a 

E*p ' u *=2 m’e' 

avec la meme convention pour le cas de/w' = o. 

On en deduit, en remplacant EV“ et E l P' uf par ces valeurs dans 
le deveioppement de F( a, //'), 

Ffw. «')==V-^7 

mm 1 

Nous retrouvons ainsi la formule ( 12 ) pour le cas ou m elm 1 sont 
differents de zero. Cette formule peut s’ecrire sous forme de pro- 
duit symbolique 

A-mm ■= £ ^ />("*« )Bpj j , 

en convenant de remplacer le produit symbolique par B^'. 

On trouve de meme sous forme de produits symboliques 

1 A 0 m' = ^B 0 — - ( Bj -H B — 1 )j ^ ~7 ’l m'~p' B p'J, 

(IS) A.*, =(2m J '"- pB ") [ Bi_ 7(B'. + B -i)]> 

( A„„ =[B 0 -f(B 1 + B_,)J [Bi-^Bi+Bl,)], 

formulas analogues a la formule ( 18 ) du Chapitre precedent. 

242. Supposonspar exemple qu’il s’agisse de developper la partie 
principale de la fonction perturbatrice, c’est-a-dire la fonction 

1 __ 1 

v / S (^-<)2 " A 

du n° 238. 

Nous observons que 
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se presente sous la forme d’un polynome du deuxieme degre par 
rapport aux exponentielles 

E iu , E-'«, E iu ', E~ lu \ 

u etant l'anomalie excentrique relative a la premiere planete et u 1 
l'anomalie excentrique relative a la deuxieme planete. Et, en effet, 
x\ , x ! 2 , x\ 2 sont des polynomes du premier degre en E*“, E w “ 
et x\, x\ sont des polynomes du premier degre en E* w/ , E~ iu \ 
En examinant de plus pres, nous voyons que ce polynome cou- 
tiendra seulement un terme tout connu et des termes en 

E= s *’“, E= /w , E- 2/w \ E =, * tt ', E 

en tout treize termes. Nous etudierons plus loin ce polynome plus 
en detail; mais, pour le moment, observons que si Ton fait 

x = E /w , y = 

F expression x % y- A 2 deviendra un polynome entier du sixieme 
degre en x ety que Ton peut appeler R(x, y), d’ou 

x\y$ A 2 = K (x,y). 


La formule (io), si Ton y fait F = I, devient alors 


04 ) 




dx dy 


xPyP' /R(#, y) 


ce qui montre que \S pp > est exprime par une integrale double 
definie dependant de la racine carree du polynome entier R. Seule- 
ment les valeurs que l’on doit donner a x et a y sont imaginaires ; 
on doit faire varier Pune et 1 ’autre de ces deux variables le long 
d’un cercle de rayon i ayant pour centre l’origine. 

On peut aussi exprimer le coefficient A mm ' par une integrale 
definie, mais cette integrale est plus compliqu^e. Nous avons 
trouve en effet 

P C d 2 F 

-4n*WA / „.= / / B-iiml+mT) du du', 

et 1’ integrate double peut s’^crire 




dx dy 


eq 


dx dy 

gpffty/l 
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oil Q ;i le meme sens qu'au n" 240 et s'eeril 


ii 




La derhee seconde de serai t encore egale a une fonction ration- 

nelle de w et de r divisee par y Rf x. y ), mais il \aut mieux partir 
directement de la formule (i i) en remarquant que 


1 — y—n 

On trou \e ainsi 


(n i 

avee 


{- S A.„ 


QE Q dr dy 


Q = 


= ff- r— 


Seulement ici Pintegrale (io) est beaucoup plus compliquee que 
Pintegrale (i4) parce que la fonction sous le signe j n'est plus 

algebrique, mais contient Pexponentielle E Q . On voit pour quelle 
raison le developpement suiyant les anomalies mojennes est plus 
compliqu^ que le developpement suivant les anomalies excentriques. 


213. 11 y a deux cas particulars sur lesquels nous reviendrons 
plus loin en detail, mais dont je voudrais des maintenant dire 
quelques mots. 

C’est d’abord celui ou les excentricites sont nulles. Dans ce cas 
il n’y a plus a faire de distinctions entre les anomalies mojennes 
et les anomalies excentriques, et Ton a Q— o; de plus Porigine 
de ces anomalies devient arbitraire et nous pouvons les compter 
a partir de la ligne des nceuds. Nous trouvons alors 

A 2 = a- -+- a 2 — i ad cos <r, 

oik tj designe Paagfe des deux rayons vecteurs, ou bien 

J J 

A*= a?-*- d 1 — %ad cos 2 T cos(w — u!)— saa'sm 2 - cos(w-h u l ) 9 
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(l6) 


\ H < x, y ) = xy j^( a- *4- a' 3 ) xy — aa' 


COS" - (X -r r) 
2 


sin 2 - f a : 2 }' 2 -i- i ) 


]■ 


244. Le second cas est celui ou Tinclinaison mutuelle des deux 
orbites est nulle; on peut alors choisir les axes de coordonnees de 
telle sorte que 

X j — ■ X g — o* 

d : ou 

A 2 = [( x\ — x' x ) -f- i( x 2 — x - )] [(x\ — x'i ) — i{ — :r' 5 )]. 

On voit ainsi que R(^, v) est le produit de deux polynomes entiers 
du troisieme degre 

R(^ y) = Ri(#,y) R2O ,y)- 


Le premier de ces deux polynomes R 2 (^, y) contient seulement 
des termes en 

x-y, xy*, xy, x et y. 


Les deux courbes du troisieme degre 

Hj = o, R-2 " o 


passent par 1’origine et ont des points a l’infini communs ; elles se 
coupent en outre en 4 points dont il est aise d’interpreter la signi- 
fication. 

A chaque point M de la premiere orbite correspond une valeur 
de u et par consequent une valeur de x\ a chaque point M' de 
la deuxieme orbite correspond une valeur de a 1 et par consequent 
une valeur de y. inequation R* = 0 exprime que la droite MM/ a 
pour coefficient angulaire i] Fequation R 2 =o signifie que MM' 
a pour coefficient angulaire — L 

Les quatre intersections des deux courbes du troisieme degr^ 
R^mRo—o correspondent aux quatre intersections (g6n<5rale- 
ment imaginaires) des deux orbites elliptiques. 

Si Ton suppose que y et par consequent M' soient donn^es, 
I’dquation R* = o est une Equation du deuxi&me degr£ en x \ cela 
s’explique parce que la droite de coefficient angulaire i qui passe 
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par M' coupe la premiere orbite elliptique en deux points M et M|. 
Les deux points M et M* se confondent, de sorte que Fequation 
— o a deux racines egales et que Ton a 

dR x 

dx ~ °’ 

toutes les fois que la droite MM' est tangente a la premiere orbite 
elliptique; cette tangente a Fellipse ayant pour coefficient angu- 
laire -f- i passe alors par Tun des deux foyers de la deuxieme orbite 
elliptique* 

De meme, si la droite MM' passe par Tun des deux foyers de la 
deuxieme orbite elliptique, on a 

-~r =o. 
dy 

Mais les deux ellipses ont un foyer commun qui est Forigine. 
Pour la droite MM' correspondante, on a 

dRi _ __ 

dy ~ dy- ~ ° 5 

ce qui (en regardant x etyjcomme les coordonnees rectangulaires 
d’un point dans un plan) montre que la courbe du troisieme degre 
R, = o (de meme que la courbe R 2 = o) a un point double. Pour 
R 4 = o, ce point double est 

e e 

7 — = » y — - ==== ? 

i + y r — e* I H- y i — e 2 

e e r 

7==’ y= — — 

i — y / 1 — e 2 i — / x — e 1 

245. Nous avons vu au n u 238 que le calcul de lapartie comple- 
mentaire de la fonction perturbatrice se d^duisait aisement de 
celui de la partie principale. Mais il y a mieux a faire, car le pre- 
mier de ces calculs est beaucoup plus ais 6 quele second, de sorte 
qu’il est preferable de le faire directement. 

Nous avons vu au n° 239 que, si l’on pose 


X = 

et pour R 2 = o 

X — 
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les trois formes <4e la partie complements! re, auxquelles on est 
conduit si Ton adopte les variables des n 0 * 2B et 44, sont respecti- 
vement proportionnelles aux trois derivees secondes de ^ . 

Gherchons done a developper W et pour commencer etudions le 
developpement de W suivant les exponentielles E ipu+-ip'u'. ii sera 
aise d’en deduire ensuite par le moyen de la formule ( 12 ) le 
developpement procedant suivant les exponentielles YJ mlJrlTn ' l \ 

Or x\ , x. 21 x’ :i sont des polynomes du premier degre en 
a., sont des polynomes du premier degre en HL~ IU ' . Done 
W est un polynome de premier degre en E ~ IU d’une part, en E — lu ' 
d’autre part, de sorte que le developpement de W suivant les expo- 
nentielles Etpu+ip'u' se composera de neuf termes seulement. 

II resulte de la que, si Ton applique la formule ( 1 2 ) au developpe- 
ment de suivant les exponentielles chaque coeffi- 

cient A mm ’ dependra seulement d’un nombre fini de fonctions 
de Bessel. On peut done, en se servant des relations de recurrence 
entre ces fonctions de Bessel, ramener ce coefficient a ne plus 
dependre que des transcendantes 

5 *!///' 6 ) } J / n '( fYl C )» 

246- Hansen prend pour variable independante Vanomalie 
excentrique u de Tune des planetes; supposons alors qu’il ait 
developpe la function perturbatrice, on Tune de ses derives, ou 
l’une des composantes de la force perturbatrice sous la forme 

F = ^ B pp > e*<p u +p' u '), 

nous avons 

L = u — e sin w, 

= u r -4- e* sin u’ . 

D’autre part, dans la fonction perturbatrice dont tous les termes 
contiennent en facteur la masse perturbatrice, nous pouvons, en 
negligeant le carre de cette masse, appliquer les formules du mou- 
vement kepl^rien. Alors l et l 1 sont des fonctions lin^aires du 
temps et sont par consequent liees par une relation Iin^aire. J’dcris 

l' = kl e, 

A etant le rapport des moyens mouvements et e une constante. 
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J 5 ai alors 
ou en posant 

et enfin 


V = ku -r- s — he sin u. 

u l = ku ~~ z, 

V = u\ — he sin u , 


ii\ — u’ — e' sin u -h ke sin u. 
Proposons-nous de developper F sous la forme 


11 viendra 


— 4^ 2 A,, 


. = f du dui. 


on en integrant par parties par rapport a u { et a u c’est-a-dire 
comine si u etait ane constante, 


-i- t\~- 


du du'. 


2A •— f f — — e- 

mm J J rn' du' L 

— 4 7I 2 A mm ' = y* B p// dw afo'. 

Mais le produit des deux exponentielles sous le signe j* peut 


Mais 

d’ou 


s ecrire 


mUi ]£—iAm , esinu mOw' Jgwi'e'siiu^ 

On trouve done finalement 


W) 


A/iirn ' — ~“7 Mpp'^m—pi — ^ /n'—ptm e )• 


247. Gyld^n cherche & developper la fonction perturbatrice 
non plus en fonction des anomalies moyennes, ni en fonction 
des anomalies excentriques, mais en fonction des anomalies vraies. 
Si v et v f sont les deux anomalies vraies et qu’on veuille developper 
la fonction F sous la forme 


F =2 C mm > 
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ies coefficients de la serie seront donnes par la formule 
— 4-rr-C mm'— f f FE -O/h'+wV) 

Si Ton pose $L iv -=x 7 E iV = cette mtegrale double se transforme 
en une autre ou la fonction sous le signe f* est une fonction ration- 

nelle de de y et de A, et ou A lui-meme est la racine carree 
d’une fonction rationnelle de x et de y. G’est ce qui arrivait deja 
dans le developpement suivant les anomalies excentriques. Les 
deux developpements presen tent done a peu pres le meme degre de 
difficulty. 

On pourrait chercher d’ailleurs des formules analogues a la for- 
mule (12) et qui permettraient de passer de Fun a Fautre. 

II faut ensuite tout exprimer en fonction de la variable indepen- 
dante choisie qui est Fanomalie vraie de Tune des planetes. Pour 
cela, il faut se livrer a un calcul analogue a celui du numero pre- 
cedent, mais pour lequel il n’existe malbeureusement pas de for- 
mule aussi simple. 

Nous avons une relation entre L et v que je puis £crire 

1= p -h cp(t>, e), 

<p(c, e) etant une fonction pdriodique de c, et de inline 

r=v'~ c'), 

et enfm 

L' — kl s, 

d’ou, en posant 

Ci = kv -4- £, 

(17) ci — c'-H <p(c', e r ) — k cp(c, e). 

11 s'agit ensuite, a Faide de l’equation (17), de developper l’expo- 
nentielle 

en serie proc^dant suivant les exponentielles 

L’analogie avec le calcul du numero pr^c^dent est dvidente. 
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248. Nous commencerons par le cas ou les excentri cites sont 
toutes deux nulles, ainsi que Pi ncli liaison mutuelle des orbites. 
L’expression A 2 prend aiors la forme 

A 2 = a 2 -\- a' 2 — 2 aal cos( l — 

coinme nous l’avons vu au n° 243. II s’agit de developper 

i 

S’ 


mais nous ne nous bornerons pas la, pour les raisons exposees au 
n° 215, et nous nous efforcerons de calculer le developpement de 
A" 2 *, 2 s etant un entier impair quelconque. 

Nous observerons d’abord que A 2 est homogene du deuxieme 
degre par rapport aa eta et d : autre part que A 2 ne depend des 
angles l et l ! que par l’interinediaire de 

cos ( l — V ) = cos ( u — u') — — ( ~ -h— V 

2 \y at/ 

Je rappelle que, quand les excentricit^s sont nulles, les anoma- 
lies moyennes ne se distinguent pas des anomalies excentriques. 

Nous sommes ainsi conduits a poser 


d’ou 


a x 

a = —,9 ~ — z; 

a y 


A-26' — a'~- s i + a 2 — a ( z -h 


■i)r 


Nous pouvons toujours supposer a < i ; il suffxt pour cela de 


P. - II. 


4 
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supposer que Ton a designe par a le plus petit des deux grands 
axes. 

Nous sommes done conduits a de\elopper 



suivant les puissances entieres positives et negatives de ^ sous la 
forme 

Jm md jmmm 

Les coefficients bfi ont recu le nom de coefficients de Laplace. 
249. Si nous posons 

F =(i — nz) (i — ■)» 

nous voyons que F ne change pas, ni par consequent F~ s quand 
on change : en ^ Or cela revient a changer k en — k ; on a done 

(i ) = b'f Ic) . 

D’autre part changer k en — k cela revient a changer i en — i] 
Fegalile precedente montre que les coefficients ne changent pas 
quand on change i en — i, nos coefficients sont done reels, el nous 
pottvons ecrire 

( a ) F -4 = bl - 4 - 2 ^ b { ^ cos k(l — V ). 


II y a entre les coefficients de Laplace cerlaines relations de re- 
currence qui en facilitent le calcul et qu’il est aise de deduire de 
Tidentite 

(3) F-* =2 *«*>**. 


Nous avons identiquement 


(4) 


<fF 0 

5 -r F-*-+- 
dz 


d F-s 
dz 


F = o 


ou 


sa (? ~ *) S b * k) z> k S k zk ~ 1 [ i+a2— = °> 
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d’ou en egalant a zero le coefficient de z k ~* : 
o = s'jl ( b/ l+l> — b/ c ~ 1 ' j -h (' i -f- a 2 ) A' b { /' — a [( k -+- i ) b 1 /^ 1 > -4- ( k — 1 ) 
d'ou la relation de recurrence 

( 5 ) a n ( s - k - 1 ) -+- a > (— s — k -M ) 4- ( 1 -f- a 2 ) k b s ® = o. 

Nous avons ensuite 

F-=[, — (. t-i)*- ]F-. 

ou 

2 ^-**= [»-»(*+ z) + *>] 
ou en egalant a zero le coefficient de z k : 

( 6 ) b { g k) = (1 + a 2 ) ^ — a ( 1} ^+I n )i 

relation de recurrence qai permettrait de passer des b s+{ aux b s ; 
mais il serai t preferable de chercher une relation de recurrence 
permettant de passer des b s aux Nous pouvons prendre deux 
relations (6) consecutives qui nous donnei’ont deux equations 
entre bf\ bf +i) et les quatre inconnues b*~f , , 6^ ,} , 

D’autre part les relations (5) nous fournissent deux autres rela- 
tions entre ces quatre inconnues. Nous pouvons done determiner 
ces quatre inconnues en fonction de b^ +i) . 

Pour arriver directement au meme rdsultat, nous partirons de 
l’identite 

(,) ;FP + *. S Q_,, 

ouP et Q sont des polynomes du premier degre; il est facile d’eta- 
blir cette identity et de determiner les polynomes P et Q. Cela 
pose, nous pouvons observer que les coefficients de la sdrie (3) 
peuvent 6tre etablis par la formule de Fourier : 

( 8 ) 2 LTzb^ = J* F~~ s jg^—i dz, 

L’intdgration doit £tre prise, par i^apport a l — depuis o 
jusqu’a 27z et par consequent par rapport & z le long d’un cercle de 
rayon ( 1 ). 
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En \ertu cle 1 id en tl te ( 'j », nous pouxons multiplier la fonction 
sous le signe J* par le premier membre cle cette identity, ce qui 
donne 





P dz 



F-* s z k+l Q dz. 


Xous pouvons transformer la seeonde integrale par une inte- 
gration par parties, en observant cjue, Integration se faisant le 
long cbune ligne fermee, nous pouvons laisser de cote la partie qui 

sort du signe / et qui prend la meme valeur aux deux limites. 

Xotre integi'ale devient ainsi : 


— L_ j" Fi -5 + i)- A Q + z k + l dz. 


Mais il est clair cjue 


P 



[(*■+■ 0Q + 



est un polynome du premier degr£ en z\ soit 


? z - 


ce polynome, notre equation devient 

2 — f F^-^p^ 4-1 y z k ) dz, 

d’ou 


(8') b<*>=$b£p> i). 

C’est la relation de incurrence cherchee. 

Grace aux relations de recurrence (5), (6) et (8), le calcul dcs 
coefficients de Laplace peut £tre ramene a celui de deux qucl- 
conques d’entre eux, par exemple 

AO u 1 
U U 

2 2 


230. Consid^rons maintenant la d^riv^e 
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En differentiant la formule (8) sous le signe 



\ient : 


d'oii 

(9) 


. db^ 

2l~ — j 

c/a 




dz, 


i db ( *> 
s d'j. 


/,(*+!> , Uk- 1) _ 

u s + 1 ^ u s + 1 2 a u s + 1 1 


ce qui permet d’exprimer les derivees des b s en fonction dcs b s+ 1 
et par consequent en fonction des b s . 

Nous pourrons done finalement, par ce mojen, exprimer la de- 
d/A h) 

rivee — en fonction de b {) ., 6!, sous la forme 

d'i ± 7 J 7 


( f 0) 


dbp 

d<x 


= P*£+Q*1, 

2 2 


P et Q etant des fonctions rationnelles en a. 

En differentiant la relation (to), nons trouvons : 


dtb'p 

dct. 2 


b 



-nS+ p 


db]_ 

"2? 


db\ 


> 


de sorte que la derivee seconde est exprimee en fonction de b\ 

et de leurs derivees premi&res, et peut par consequent, par une 
nouvelle application de la formule (io), £tre exprimee en fonction 
de b\ et b\ seulement. 

IT 2 

On op^rerait de raeme pour les d£riv6es d’ordre sup^rieur. 


2ol. Nous avons 


.,>-.(,-1)- 


Or la formule du binome nous donne 

/ , . S (5 4 *l) „ „ 

(i — az )-~ s = i4- su.z H ol*z-- 


ou plus gen^ralement 

. n e V 1 r (s-+-p) 

(l — LL. 




olPzp, 
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OU 


u F repi'esente le fonction eulerienne. De raeme 


Vi s -f- q ) 




^ ^ y ^T(s)Y(q-^\)' 

ou, par multiplication, 

F-^y T(* + p)?(s+_q )_ %P + qzP -a 

^P(5)r(^ + i)r( ? + i) 


ce qui nous donne pour le coefficient de z k : 


<") 


b , k) __ + ^ 4 - k) 

' s -^*1 r 2 ( s ) r ( q -r- i)T (q -t- k -f- I ) 


X Am “ ? 7, 


formule qui donne le developpement des coefficients de Laplace 
sui\ant les puissances croissantes de a. Nous remarquerons tout 
d'abord que b (k) est divisible par a*, et que c’est une fonction paire 
ou impaire de a suivant que k est pair ou impair. 

Rapprochons la serie (i i) de la serie hypergeom^trique de Gauss. 
Cette serie s’ecrit : 


F(A, B, C, X )=^ 


r(A + g )r(B + g )r(Cj 
r(?+i)r(C + ? )r(A)r(B) 




La comparaison nous donne immediatement 


b^ = 


T(s -f- k) 

rJTYrl ITT) 


a /c F(s, s -h /c, £-+-i, a 2 ). 


232. On sait que la serie hypergeometrique de Gauss satisfait b 
une equation lineaire du second ordre, il doit en £tre de mfine 
de b (k) . C’est en effet ce qu’il est aise de verifier. Nous voyons en 
effet que dans la serie (i i) le rapport du terme general au prdcd- 
dent est 

(s-hq — i)(s^-q~\~k-~i) a 
q(q-^k) “ ' 

Si done j’appelle C 7 le coefficient de x h+ -</ : nous aurons 
q(q +/f)C ? =(5 + ^ — r)(^-i-^~f-/c — i) C 7 ~i, 

d’ou 

( 12 ) ^ “+■ k) C ? a2?+*=^(s -+- q - 1 ) (s -+- q -+- k — 1 ) C^-i «*?+*. 
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Or, si nous observons que 


b =Vc ? ^- 2 + /:) c '/ 22 ' 7+/ ‘’ 

“ 2 ^t =2( 2? ' hA ’) + k-i)C q &v+i> 


55 


et que de meme 


a 2 


6 = 2 C7_iaI ’ H ’*’ =^(2? + ^ — 2)C ? _,a2'?+*, 

a'*^T + k — •>■)(?■ q -+* a- — 3 ) Cj-i ***+■*; 


si nous remarquons de plus que les coefficients du premier et du se- 
cond meinbrede(ia)quisont<jr(y + £), (5 + q — 1) (.? -+- q -h A~ — 1) 
sont des poljnomes du second degre en q et par consequent peuvent 
s’cxprimer lineairement, le premier a Paide de 

I, ‘iq k, {rtq H- /c) (2 q -h /l — 1), 
le second a faide de 


iq-\-k — 2, (2^ -h £ — 9.) (a^ -4- A- — 3 ), 


les coefficients ne dependant que de £ et de nous verrons qu’il 
y a une relation lineaire entre les six quantites 


7 db d*b , „db 
b ’ a Ta } *-d*’ * 3 d~*’ 


a 4 


d*b 
dtf 9 


relation dont les coefficients ne dependent que de k et s . C’est 
Pequation diffdrentielle cherchee. 

J’ai supprime pour abreger les indices s et k de bf\ 

Cette equation sMcrit 

(> 3 ) + l a “( 4 « + i)a 3 ]^ ~[ 4 « 2 a 2 + k*(i- a 2 )] b= 0. 


Elle peut nous indiquer comment se comporte la serie (11) 
pour les valeurs de a voisines de 1. Les methodes de Fuchs nous 
apprennent en elfet que les intdgrales de cette Equation ne peuvent 

presenter de singularity que quand le coefficient de s’annule, 
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c’est-a-dire pour 

z = o, a = i . 

Pour a = o, les racines de Inequation determinante sont ±: Aq 
ce qui veut dire que Tequation admet une integral e particuliere 
developpable suivant les puissances de a et commencant par un 
terme en y. k (cette integrale iPest autre chose que la function b ( s A) 
qui nous occupe) et que Pintegrale generale est de la forme 

S -+• hb { s L) logs, 

h etant une constante quelconque etS une serie procedant suivant 
les puissances entieres croissantes positives ou negatives de a et 
commencant par un terme en a ~ k . 

Restent les points singuliers a = rfci; je me contenterai d’exa- 
miner le point 2 = + 1; puisque Pequation differentielle ne change 
pas quand on change a en — a. 

Pour a = +i, les racines de liquation determinante sont 0 
et 1 — 2 s ; ce qui montre que, outre une integrale particuliere deve- 
loppable suivant les puissances de a — 1, Pintegrale generale est 
de la forme 

S h- P log(cc — 1), 

ou P et S sont developpables suivant les puissances entieres erois- 
santes de a — i, le developpement commencant pour P par un 
terme de degre o, pour S par un terme de degre 1 — ■ 2$. 

Pour $ = ~ » onai — 25 = 0 etS ne devient pas infinie, c 7 esl le 

second terme qui est preponderant, de sorte que bf ] devient infini 

pour a = 1 de la meme maniere que log(a — 1). Pour s = 

le terme S devient infini et preponderant de sorte que bf ] devient 
infini de la meme maniere que 

(ct — I ) 1 - 25 . 

On voit par la que la serie (1 1) converge pour 

M<i 

et diverge pour |a|>i. Les resuitats precedents se deduisent 
d’ailleurs immediatement des proprietes connues de U serie hypcr- 
geometrique. 
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233. On a propose nn tres grand norabre de procedes de calcul 
pour les coefficients de Laplace, mais il j en a un qui est tres su- 
perieur a tous les autres, e’est celui qui est fonde sur Temploi des 
fonctions elliptiques. 

Nous Scions que la function p(u) de Weierstrass est definie par 
l’equation 

Lp'0)P = 4P 2 (*0 — ffip(u) — grs— 4(P — e t ) (p — e % ) (p - e 3 ), 
ou e i: e 2: c% sont trois constantes liees par la relation 


avec la condition 


<?, -h e- 2 -+- e z — o 
p(o) = co. 


Nous savons en outre qu’elle satisfait aux conditions 

p( w ) = p(— M)=pftt-h^to 1 ) = p(^H-2to 2 ) = p(i^-+- *20)3); 
p(u 1 )=e ] , p(w 2 ) = e 2 , p(co 3 ) = e 3 . 

D’autre part la fonction £(«), qui est telle que 

= — PO) 

jouit de la propriety 

till) = — £(— u), Z(u-+-‘2.u t ) = Z(u)-\-'Lr li , 

Tll= ?(«!),), 71l + ^ + ^= O- 

Reprenons rinlegrale 

z k dz 


ou 


b\* = 


i r £ k 

J i / z( \ — a 




et plus generalement la suivante : 

* z k ~ x dz 1 


bsp 


i r z k ~ x dz _ i r 
o.iiz J F s ~~ 2 itc J [z( 1 


z /fc~l-4-2s d z 


■ ctz) ( z — a)] s 


Pour identifier aux formules des fonctions elliptiques, nous de- 
vons poser : 

Z = p ( u) — 63, a = 65 — 63 , — — C\ — £3, 6 1 -f- 6J2 -h 63 ” 0 , 

Of 


d’ofi 


■a) = — P'* 2 0), 

4 


— xz) (z 
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d'ou 

01 ) (jT~ j *i~b* k = f (p — e*) A ~ l ~ hi *P' 1 ~ ss d u 

et, en particular, 

(i5j Try/a ^ = /* (j> — ^ 3 )^-i-+-2,? 

2 

231 . Rappelons une propriete importante des fonctions double- 
ment periodiqties : c’est la possibility de decomposer ces fonctions 
en elements simples. 

Soit F(«) une fonction doublement periodique, soient a,, 
flr 2 , . . ., a u ses inlinis distincts; bien entendu je ne regarde pas 
comme distincts deux infinis qnand ils ne different que par des 
multiples des periodes 2 a >/. 

Soit ay Fun de ces inFmis qui sera, je suppose, d’ordrc /r ; deve- 
loppons F ( u) suivant les puissances croissantes de u — a J: et soient 

A a- A/ L — ] ^ A 2 ^ Ai 

( a — aj) k (u — ctj Y‘~ l ( u — aj ) 2 u — a y 

Fensemble des termes de ce developpement qui ont un exposant 
negatif. Posons 

± ?<*“»(« — aj), 

ou Q k) (u) represente la derivee k ihwc de %(u) par rapport a n. 

De cette fagon la difference 

F(m) — *j(U) 

reste finie pour u = ay. Nous aurons alors 
(16) 

C etant une constante. C’est cette formule bien connue (16) qui 
represente la decomposition de F(«) en 4l4ments simples. 

Revenons aux integrates (14) et (i 5 ) ; elles doivent 4tre prises 

de o k 2(0,. Si nous d 4 signons la fonction sous le signe 
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F(w), nous decomposerons FuF) en elements simples, et nous 
integrerons separement ehacun de ces elements. Parmi ces ele- 
ments il y en a qui nous donneront zero, ce sont les elements 


f — aj)du , 


oil k^. 2 ; Pintegrale indefinie nous donne 

nk-i)( u — ay), 

qui est (au signe pres) la fonction p{u — aj) si k — 2 et une de 
ses derivees si k > 2. Dans tous les cas ce sera une fonction dou- 
hlement periodique qui reprendra la meme valeur aux deux 
limites o et 20, ; Pintegrale definie est done nulle. 

Si k = 1 


I ~ K'( il — a j) du. = X^{u — aj)^ 


et l’integrale definie est 


C(aci>i — «/)“ ?(—«/) = ^1. 


Si k = o, on a 


J* t^(u — aj) du = logs' ( u — cij), 


et Pintegrale definie est 


logG^S CO! — dj ) 

Jogtf( — aj) 


— iiz -+- ‘2 7)! (cOj — a,j). 


INous n’aurons pas d’ailleurs a faire duplication de la formule 
(17) dans ce Chapitre. 

Reste enfin la conslante C qui nous donne 

f G du = ^GtOt. 


Dans le cas qui nous occupe, et qui est celui des integrates (i 4 ) 
et (i 5 ), la fonction F (u) ne peut devenir infinie qu’avec p(u) ou 
p'(u ), e’est-a-dire pour 


o>j, co 2 ou co 3 . 
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J'ajoute que, pour s ^ [formule (i 5 )], elle devient infinie 
seulement pour a = o. 

De plus, Pexposant i — 2 s est pair, et, com me p' 2 (tt) est une 
fonction rationnelle de p(u)* il en sera de me me de F(w). Done 
F(w ] sera une fonction paire de aq ainsi que de u — o>/ 

F(u) = F(— «), F( u ) = F( 2to, — m)* 

Done le developpement de F(^) suivant les puissances crois- 
santes de u ou de u — <x>i ne contiendra que des termes d’ordre 
pair; il ne contiendra done pas de ter me de degre — 1 [qui aurait 
pudonner un element simple en ^(u) ou — o> 4 *)] ; et e’estpour 
cette raison que nous n’aurons pas a faire usage dela formule (17). 
Nous n’avons pas a nous inquieter des termes de degre — 4 i 
— 6, .... qui, comme nous Pavons vu, donneraient une integrale 
nulle, mais seulement des termes de degre — 2, ainsi que de la 
constante C; nous trouvons ainsi pour notre integrale 

(18) ft lie (^2) & ( S A) = ftCcdi — 

V 

^ Ao etant la sojmxne des coefficients des termes de degrd — 2 
dans les developpements relatifs aux quatre infinis 0, co 1? to 2 , 003* 

2oo. Ces coefficients C et A 2 sont des fonctions ralionnelles 
de a. En effet le developpement de p(u) ou de p f {u) suivant les 
puissances croissantes de u ou de u. — c oi a pour coefficients des 
fonctions rationnelles de e i} c*, c 3 et, par consequent, de a; il en 
est done de meme du developpement de 

F O) =(p — 61 ) 6 - 1 +** 

Done tons nos coefficients A 2 sont rationnels en a. 

Parlous mainlenant de C. Nous observerons que F non seule- 
ment ne peut devenir infinie que pour a = 0, (o l7 u> 2 , co 3 , mais ne 
peut non plus s’annuler que pour Pune de ces quatre valeurs de u. 

Il en resulte que F(u) ne peut pas admettre les quatre infinis 
mais s seulement une partie d ? entre eux, ce qu’il est cPailleurs ais^ 
de verifier; pour celles de ces quatre valeurs pour lesquelles elle 
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ne devient pas infinie elle s’annulera de sorte qae l’on aura 

(19 ) O = G — 1 ^ ^ /• ' A r ! — «■>/) 

pour u = o, o) t , co 2 , o> 3 ; et pour w* = o, to,, co 2 , to,,. D’aiileurs on 
n'aura pas u = to M puisque w est une des valeurs pour lesqueiles F 
s'annule, et co/ une des valeurs pour lesqueiles F est infinie. Done 

u — co 2 ou o)j (a une periode pres). 

Mais Aa est rationnelle en a. II en sera de meme de £W(co,), 
^(u>D, s (A, ( w di q u i sont des derivees de p(u) pour u = to,, to 2 
ou to 3 , et nous venons de rappeler que ces derivees sont des fonc- 
tions rationnelles de e t , <? 2 , e mi et, par consequent, de a. L’equation 
(19) montre done que C est rationnel en a. c. q. f. d. 

En particular nous avons pour 

* = k = °'< F ( U ) =I > 

~~ r 
2 IZyOL 

k = i, F (w) = p(a) — e 3 = — e. } — £'(«)> 

- a „+| Wl ( a+ i) 


Pour un coefficient b [ s h] quelconque nous trouverions de meme 
a(ao _ P-ni-*- Q^i 

ou P et Q sont des fonctions rationnelles dc a. On ddterminerait 
aisement ces fonctions rationnelles a 1 ? aide des relations de recur- 
rence du n° 249 . 

256 . 11 reste a calculer to, et *e, d’oii b° { et b \ ; pour cela je 
renverrai aux formules et propositions pour l’emploi des fonctions 


d’oi 


pour 

d’ou 


2 ’ 
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elliptiques redigees d'apres W eiers trass par M. Schwarz et tra- 
cluiles de 1 'allemand par M. Pade (Paris, Gauthier-A illars, i 8 g 4 )* 
Portons-nous a la page 61. Deux cas sont a distinguer, et d’abord 
celui ou 

, , 1 1 

e ± — £3 < ei — e->, cest-a-dire a < a. a<— • 

a V /l 2 

Dans ce cas nous poserons 

t/ e , — e 2 1 — . y r — a 2 

1 = — - — = — - — 1 

\C\-~6% -+* \/^I — C-2 I v 1 K 2 

et alors, en supposant 

‘TKO.tt 

A = E W i , 


ce qui est le q de Jacobi, nous aurons 

l — ^ + 2/^+... 

1 H- 2 A 4 --}- 2A 16 -f-. . . ? 

d’ou 

( 20 ) A = - -h 2 ( - 

2 \ 2 

On calculera A par la formule (20) et Ton Lrouvera ensuite 


i5 


‘'v.so a 


121) 


\J a/>l= Y~~r- V* = 2 (a* - f- &*-+- A 4 + ...), 
v/ 3 !"" [/' =‘ + 2A + 2A<4-uA» + ... ) 
y/(i \/i -« = i-aA + aA»-a/i» + .... 


Ce sont les formules (6), (7), (8) de Schwarz, mais nous ferons 
usage pour le calcul de b'\ de la formule 


(22) \/&? = - JL. (i-4-2A*H-2AI6h-...), 

qui est la formule (4) de Schwarz, et pour le calcul de b\ de la 
formule 2 

1271,0)! I — 3 3 /l 2 H- 5 5 A 6 — . . . 

^ i — 3 A " -4- 5 /( ® — . . . ’ 


qui est la formule ( 5 ) de Schwarz. 
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Le second cas est celui ou a > — • On se ser\ira alors des 

V/2 

formules (18) et suivantes de Schwarz et Ton posera 


/ — y/gt— e- A — \J e 2 — gj _ i — - y/q 

V e i~ e i-r- y e *. — e z i -f- \J^ 


h { = E «> 3 , 


et Ton retrouvera la formule 


(20 bis ) / 

analogue a (20). 

Nous aurons ensuite 

/ y - v , A w 3 

(22 £>£$) 4 / - 


t A / A 
/2-1 = h 2 ~ 

2 V 2 


■(i + aA} + ...), 


b\=™L = ('*\ (log * ) * 

2 tt / a \™/ \ J ti y a 


formules (19) de Schwarz. Quant a yi, et y] 3 et par consequent b\ 
on les deduira des formules 


(23 bis) 


n* 1 — 33 h\ 

2 r i3 o) 3 = — i , 

b 1 — -3/if-H... 

w 3 — ( 0 ! y ) 3 ■= 7 T i» 


On remarquera que co 3 et r j;{ sont purement imaginaires. 


257 . 11 importe de se rendre compte de la rapidity de la con- 
vergence; celte rapidite est extreme; en eflet, si a< on a 


fa-i 


h<E-K, 


et si a — ? on a 


7 ^ ^2 — 1 

y *>< v= — ; 

* y 2 *+- I 


hi < E- 71 . 


Suivant le cas, L ou A sera <0,08; tandis que h ou h K sera 
< o,o 4 - Or les series qui procedent suivant les puissances de h et 
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qui ne sont autres que les series 0 de Jacobi ou des series ana- 
logues, et qui ne contiennent par exemple que des tennes ou 
Fexposant de h est un carre parfail, convergent tres rapidemenl. 

Les series (20) et (20 bis), qui donnent h et /?,, convergent 
aussi tres rapidemenl, puisque dans le cas le plus defavorable les 
tennes successifs sont respectivement plus petits que 

1 r 1 1 

io ? io 6 ’ io 10? 10 15 


Et meme en negligeant 2 A 5 , c’est-a-dire .; 0 ^ o l>0 , 011 a pour 

_ 

\^i= — ==, 

V 2 H - y l — a 2 


et pour a > 


A 


6 ? = 


TZs/cL 


2 . 2 -+- 2 \J1l 

— log 


(tft-K 


[ — /a 
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LEb POLYNOMES DE TISSEHAND. 


5o8. Nous allons examiner maintenant le cas ou, les excentri- 
cites etant toujours nulles, l'inclinaison mutuelle des orbites est 
quelconque. Dans ce cas, nous pouvons ne pas faire de distinction 
entre les anomalies vraies, moyennes ou excentriques et nous 
avons trouve au n° 543 

A 2 = 61 *-+- a' 2 — <iad cosar, 

J , . J 

cos a 1 — cos 2 - cos( a — u ) 4- sin 2 - cos( a -r u ). 

*2 7 2 7 

II s’agit de developper 

i 

et nous trouvons d’abord, en appliquantla formule ( 2 ) du Chapitre 
precedent, 

( 1 ) ^ ~ ^ 4 - cos k 7. 

II reste a developper cos A* a-; c’est ce qu’a fait Tisserand en 
employant les artifices suivants : 


2o9. Posons 


cos 



sin 2 - = v, 
2 


COSd = (JL COs£ 4- v COS 7 ), 


£ = U — u! . 7 ] u- 1- 

# = y ~ E*“\ 


z ~ , w — xy , 2 cos £ = £.4- 2 cosy] ,= tv 4 - mp -1 , 

Z = F -6 ' = (1 — a}2 z — afx z ~ l — otv w — av w p ~ 1 4 - a 2 )~ 5 . 

P. —II. ,5 
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La formule du poiynome, generalisation de la formule da 
binome, nous donne 

Z=Ya<- ot[iz \ a (— aa - 1 )P(— av w ) c ( av w~ l )<? (a 2 )* 

avec 

r ( i — s ) 

<2) A -ni-fl-i- c -(;-e-s)r(a+ijr(/)+i)r(c+[)r( ? +i)r(e+i)’ 

#, o, /?, £ etant des entiers positifs quelconques et F etant la 
fonction eulerienne, ou bien encore 

(3 ) Z = 2 A ( - I za+p-hc+q-+-2e ^a-hp yC+-q z a-p W C-q % 

Nous voulons developper Z, d’une part, suivant les puissances 
de a et, d’autre part, suivant les cosinus et les sinus des multiples 
des anomalies ou, ce qui revient au meme, suivant les puissances 
positives et negatives de z et de w. Nous chercherons done le 
coefficient de 

a m z h w k % 


ou m est entier positif, h et k entiers positifs ou negatifs. Nous 
ferons done 

m = a -h p H- c q -+- a <?, 
h — a — p, k = c — q. 

L’exposant a-j-p de u est plus grand en valeur absolue que h 
et en differe d’un nombre pair; Fexposant c + q de v est plus 
grand en valeur absolue que k et en differe d’un nombre pair. La 
somme de ces deux exposants est plus petite que m et en differe 
d’un nombre pair. Le coefficient cherclffi est done un poiynome 
entier en p 2 et v 2 de degre 

™-\h | — 1 *| 

2 

multiplid par le facteur 

( — I 'jh-’rk jjJAI v^l. 

Nous observerons en effet que 

( — i — ( — i )a-p^-c—q _ ( j^A-wfc, 
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Supposons d'abord h et k positifs de facon que h = | h !, k = ] k ] 
et formons le polynome P en question de sorte que le coefficient 
cherch£ soit 

( — i)h+kyhyk'p ' 

Faisons done 


d ; ou 


m — h—k — 2g , m -t Ji k = 2/, 
a -T- p = h -h ‘ip, c + q = k + 2^, 


m + h + k 

a- 4 -p-hc- J r q-be = h p H- q, 


m — h — k 

e = p-q\ 


il viendra, en se reportant aux formules ( 2 ) et (3), 


( 4 ) 




T( 1 — s) 


r{j-/—s—p—q)r(/i-hp-i-i)r(k~hq-hi)r(p-hi)r(q~hi)r(i-hp—p~q) 


Nous observerons que m et h 4 - k sont de meme parite et, par 
consequent, que /, £*, a -f s et [3 sont entiers. Cela va nous per- 
mettre de transformer la formule; et, en effet, si l est un entier, 
on a 


( 5 ) 


r( i — s) r (* + i) 
Y{i-s-l) Y(s) 1 


Cette formule resulte de la relation connue 


r(Or(i-*) = 


% 

sin sit 


et du fait que, si l est entier, on a 

• sin (s -h /)/? = (—- 1 ) 7 sinsTt. 


260. Le polynome P n’est autre chose, a un facteur constant 
pres, que Fun des cas particuliers de la serie hyperg^ometrique a 
deux variables ^tudiee par M. Appell. Introduisons, en effet, la 
notation de M. Appell, en posant 


K 


m) =! 


T W 


r(i--a) 


(- 0 *- 


La premiere expression se pr^sente sous une forme indeter- 
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minee quand a et a-j~ m sont des entiers negatifs; mais il est aise 
de \oir que la vraie valeur est alors le produit 

( a -+- m — i ^ ( a -h m — 2 j . . . ( a -f- 1 ) a 

et d’ailieurs la deuxieme expression est alors determinee. 

Nous aurons alors 


(A -HI, p) = 


= T(p-hi'), (<,q) = T(g-r-i), 

r( h -f* p -t- 1 ) 


d’ou 

( 6 ) 

et 

(7) 


r(A + i) 

i—g-p + q) = (— 1 ) /1+9 


. t V (k -p q 1 ) 

(A + 1 '^ = - f(A- + ,)- " J 

r(i— /— s) _ 

Td—f—s—p — q)’ 

r (i + i'i 


p = 


r(i - 1 -g — p — qY 
r(i — i)Pi 


p ,=2 


r ( 1 — / — s) T( h -+- 1 ) r(/c-Hi) T( g 1 ) 




[JL 2 ^V 2 7 . 


On reconnait, dans Pexpression ( 7 ), la serie hypergdometrique 
a deux variables de M. Appell. 

Nous avons suppose dans ce qui precede h et k positifs, mais 
nous n’avons pas ainsi restreint la g<£n£ralite et, en effet, Z ne 
change pas quand on change, soit sen r 1 , soit w en ; d’oii il 
suit que le coefficient de 

c X m %h yyk 


est le m£me que celui de 

CL m Z~k cc Ttoj&h w~ OL fn Z~h- W~k . 


261. On sait que la fonction de M. Appell satisfait k deux 
Equations aux derives partielles; on doit s’attendre 4 ce qu’il en 
soit ainsi de notre poljnome P qui n’en difffere que par un facteur 
constant; c’est ce qu’il est aisd de verifier. Soient, en effet, Z' et 
Z" les ddriv^es premieres et secondes de Z par rapport k 

F = a* — 2a(n cosj; + v cost)) -I- i, 
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il viendra 




(8; 


dZ 

JJL COS; V COST, ) ; 


— j — = 2 Z ( a — 
da 

(9) 

dZ 

= — 2 Z' a COS;, 

dZ 

d\ x 

— T- = 2 Z a COST, 

dv 1 

(jo) 

\ 

d' Z 

-y— = 4 a 2 a 2 Z" 

«r 

sin 2 £ 4 - 2 sc{jlZ' cos£, 

j 

Z 

sin 2 T] 4 - 2 av Z' cost], 


1 

^ =4 * 2v2Z 


(II) 


(12) 


( 1 3 ) 


d*Z 


d* Z 


| = 4 = s Z"cos^. ^= 4^0085 008 ,, 

i d^L 

( = ^ V COS 2 1), 


d~ Z 

-^—5- = 4Z ff ('a — cos cr j- —{— aZ', 


^Z 
c£a d\j. 
d* Z 
rfv 


= — 4 ^Z"(a — cos<T) cos£ — 2Z' cos 


= 4 ^ Z" ( OC — COSO-) COS 7] — 2 Z' COST), 


ce qui nous monlre que ces onze derivees partielles de Z peuvent 
s’exprimer lineairement en fonctions des neuf quantitdssuivantes : 

/ Z', Z'cos£, Z' cost], 

(i4) | Z", Z ff cos£, Z"cost], 

( Z"C0S‘ 2 £, Z” COS £ COST], Z"C0S 2 7], 


les coefficients £tant des polynomes entiers en a, jx, v. 

D’ailleurs, ces quantity (i4) ne sont pas ind6pendantes, elles 
sont li^es paries relations 

/ (i + a 2 )Z ff — 2 «|jlZ"cos£ 
l —2 avZ ff cosT) = — (5 + i)Z r , 

J (i 4- a 2 )Z" cos£ — 2oc(j. Z n cos 2 £ 

( 1 5 ) ( 

j — 2 av Z" cos £ cost) = — (s 4- i)Z* cos£. 

f (l 4 “ CC 2 jZ ;/ COST] — ' 20 L\xV COS^ COST] 

\ — 2avZ"cos 2 7] = — (s 4 - 1 ) Z' cost). 


li est ais^ de comprendre i’origine de ces trois relations. 
De 
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nous deduisons 

FZ" = -U-i- \)Z\ 

En remplacant F par sa valeur, on aura la premiere relation (i 5 ) 
et l’on aura la deuxieme ou la troisieme en multipliant la premiere 
par cos^ ou par cost,. 

Mais nous pouvons aller plus loin; les relations (8) a (i 3 ) nous 
montrent que les onze derivees 



dZ 

dZ 

dZ 

d' 1 Z 

d' 1 Z 


1 


dff 

dv’ 

dP 

dr { 2 ’ 


j d*-Z 

d*Z 

d- Z 

d* Z 


d* Z 

d' 1 Z 


d\i dv ’ 

dr - ? 

a “ dv.' 1 

> a 

dv dff 

a dv d'> 


s’expriment lineairement en fonctions des dix quantites 

/ Z'a 2 , Z'a cos £, Z'a cos iq, 

( 17 ) | Z"a 2 , Z"a 2 cos 2 £, Z" a 2 cos $ cost), Z"a 2 cos 

( Z"a 4 , Z r a*cus5, Z'^cost], 

les coefficients etant des polynomes entiers en p et en v inclepen- 
dants de a. D’autre part, la premiere relation (10), multipliee 
par a 2 , sera une relation lineaire entre les dix quantites (17); 
nous avons done douze relations lineaires entre les vingt-deux 
quantites (16) et (17); en eliminant les dix quantites (17) il nous 
restera deux relations entre les quantites (16). 

II y a done , entre les onze derivees (16), deux relations 
lineaires dont les coefficients sont des polynomes entiers en p 
et en v. 

Or nous pouvons poser 

Z = 2 u 0 L m z h w* =5 2 u E 

avec (si par exemple h et k sont positifs) 

U = ( — t) h + k phvkp. 

On voit alors que, dans les onze quantity (16), le coefficient de 
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se reduit respectivement a 
m U , 


dU 
du. ’ 


dU 

d't 9 


• A 2 U, - ^U, 
dU 


, » T d\J dU 

— - > -7 — ~t- • — rr> m(m — i)U, m-r-, m— r-« 

ajjt 2 ap. av av 2 ap av 


Si done nous prenons Tune des deux relations lineaires dont il 
vient d’etre question et qui lient les derivees (16) et que, dans 
cette relation, j’egale a zero le coefficient de 


j’obtiendrai une relation lineaire entre les derivees 
. TT dU dU d 2 U d 2 U d* U 

1 ^ 5 dp . 7 dv 9 dp 2 y d]±dv’ dv 2 

Done le coefficient U satis fait a deux equations lineaires aux 
derivees partielles du deuxietne ordre, dont les coefficients 
sont des polynomes entiers en pi et v. 

Ce sont les equations de M. Appell. 

262 . Je rappelle la forme de ces Equations telles que M. Appell 
les a dtablies dans son Mdmoire du Journal de Liouville (p. 182, 
1882). 

(x — x' 2 )r — y-t — 1 xy s 

+ [7 + (a + (3 •+■ i)x]p — (an- (3 -+- 1 )yq — a fiz = o, 

(y — y 2 )t — x 2 r — *ixys 

[7' — (a+P + i )y] q — (a+f + \)xp — <t$z = o. 

Dans ces equations x et y sont les variables; z est la fonction 
inconnue,/> et q sont ses deux derivees du premier ordre; r, s , t 
ses trois ddriv^es du deuxieme ordre. 

Pour passer de ces Equations a celles auxquelles doit satisfaire U, 
il suffit de faire 



x = p‘ 2 ; y — v 2 , z == jjl ti v tc Uj 

«=/+«, p =— i, 

Y = h -b 1 , 7' = k -+- 1 . 

Mais une premiere observation se pr^sente. Les Equations li« 



CHAPITRE XVIII. 


7‘ 2 

neaires simultanees ng) ne comportent pas une solution unique; 
elles component quatre solutions lineairement independantes, 
ainsi que Fa montre M. Appell. Quelle est alors. dans le cas qui 
nous occupy, la signification de ces di\ erses solutions? 

Reportons-nous a ce que nous avons dit au n° 242 , et appli- 
quons des principes analogues au probleme qui nous occupe; nous 
verrons que le coefficient de z h w k dans le developpement de Z 
sera egal a 


(20) 


zf.r* 


hdz dw 


nk-ri 


L'integration doit etre prise par rapport a ^ le long d’un cercle 
avant pour centre forigine et pour rayon l’unite et par rapport 
a w le long d’un cercle analogue trace dans le plan des w. La 
combinaison de ces deux cercles definit un contour ferine a deux 
dimensions le long duquel doit etre prise l’integrale double. En 
d'autres termes, le coefficient cherche est une periode de dinte - 
grale double (20). 

Pour avoir U, il faut developper ce coefficient, qui depend de a, 
suivantles puissances croissantes de a et conserver le coefficient 
de a' 72 . 

Mais Fintegrale double (20) possede d’autres periodes que celle 
que nous considerons; une quelconque de ces periodes est fonc- 
tion dea et peut etre developpde suivant les puissances croissantes 
de a. Le coefficient de cr. m dans ce developpement satisfera aux 
m£mes equations diflerentielles de U. La multiplicity des solu- 
tions des equations (19) s’explique done par la multiplicity 
des periodes de Cinte grale (20). 


263 . La serie hjpergeometrique de M. Appell 

> ■ — i — x ’ n Y fl 

m ) (t i n ) (>» m ) ( 1 , ») J 

peat £tre expriinde par le mojen d’une intdgrale ddlinie, mais je 
me bornerai sur ce point a quelques brdves indications. Conside- 
rons l’mtegrale 


ut 1-1 ( 1 — w)v-i da 



LES POLYXOMES DE TISSERAND. 


73 

que nous prendrons le long d'une ligne allant de Tinfini a 1 ’infini 
en passant entre les deux points o et 1 . 

On pourra faire. par exemple, u = - j- in' et faire varier u f par 

valeurs reelles depuis — 00 jusqu'a -f- oc. Cette integrale sera finie 
et determinee pourvu que 

P + 9 < 

Dans ces conditions, & un facteur pres C facile a determiner 
et qui ne change pas quand p on q augmente d’un nombre 
entier , notre integrale est egale a 

T(p)T(q) 
r(/>H \-q) ' 

D’autre part, notre serie hypergeometrique (a un facteur cons- 
tant pres, independant de x et de y comme de m et de n) est 
egale a 

V"’ r (' 1 — y — m ) r ( 7 — a — 11 ) r < 1 — *(' — n ) T [ y — p — m ) 

2 mk r ( 1 — a m — n) T ( 1 — p — m — n) 

( 1 -h j 3 — Y* m ) xm ( 1 -+- a — y , n ) v 11 
U ,m) ( 1 , n ) 


On voit que chaque terme sous le signe ^ est decompose en 
quatre facteurs; le premier facteur c’est Pint^grale 


/ 


U~l~ m ( I — U du, 


k un facteur pr&s C qui est independant de m et de n, puisque m 
et n sont entiers (et que le facteur C, ainsi que nous venons de le 
remarquer, ne change pas quand Fun des exposants augmente d’un 
nombre entier). 

De m£me le second facteur, a un facteur constant pr&s, c’est 
1 ’int^grale 


/ 


p-y'-n 


(1 — )r^ 1 P m dv. 


Remarquons que toutes nos int^grales sont fmies, au mo ins 
pour les valeurs de a, |J, y^ y 7 , qui satisfont k certaines megaliths. 
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Le troisieme facteur est un terme du developpement de 
( i — x jT'-p " 1 = V k m x m 

et le quatrieme est un terme du developpement de 
(i — v )Y— ^B,^. 

Nous trouvons done 

^ ^ f du dv u~y~ m l i - u) r-i-a-* - e A OT tf OT B*/«, 

ou bien 

2 // rf “ * “~ r( 1 - “ ^~ x p_r - ^ r " _p A "' [siTTTj ] m B " [wT^j] 

ou bien 

f jdu dv K-Y(i - «)Y-« p-r(. - P)T’-Hl[i T7j T_H [' -^rJ T_1 ' 

ou enfin 

(21) J J* dudv u a v b (u~uv — x) c (v — uv — y) d 


a= i-H P — Y — £ = i + a — y — y', 

C = Y r — P — I , <3? = y — i — a. 

L’integrale doit £tre prise tant par rapport a u que par rapport 
a 9 le long d’une ligne allant de Finfini a Finfini en passant entre 
0 et 1. Les autres solutions des equations (19) pourraient etre 
represents par cette m£me integrate (21) prise le long d’autres 
contours convenablement choisis. 


264 . Envisageons toutes les series 


Myu 


(an m + n)( ft t , m + n) f 


ou les eonstantes a f , p 4 , y, , y' sont 6gales k a, (J, y, 8 augments 
dun nombre entier. Je dis que toutes ces series pourront s’exprimer 
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lineairement a Faide de quatre d'entre elles, les coefficients de 
Fexpression lineaire etant des fonctions rationnelles de x et de y. 

Soient, en effet, z x , z*, z 3 , quatre solutions des equations (19); 
la solution generate de ces equations sera une combinaison lineaire 
de ces quatre solutions; quand x et y reviendront a leurs valeurs 
primitives, apres avoir varie d’une maniere continue, il pourra se 
faire que ces solutions ne reviennent pas a leur valeur primiti\e, 
si les variables x et y ont tourne autour d'un point singulier. Mais 
alors les solutions z u z 2: 5 3 , auront subi une transformation 
lineaire a coefficients constants; c’est tout a fait analogue a ce 
qui se passe dans le cas des equations differentielles lineaires 
ordinaires. 

Cela pose, considerons quatre systemes d’equations analogues 
aux equations (19), mais ou les constantes a, [ 3 , y, y' aient des 
valeurs differentes, et supposons que la difference des valeurs de a 
pour deux de ces equations soit un nombre entier (et de meme 
pour j 3 , y, y'); soient 




-r(l) 

* 

z\ l >. 


4 ”, 

*(2) 

4 2> 

r(3) 

1 1 

4 ”. 

5j(3 ) 

Z 3 7 

4 *’: 

*(4) 
Z 1 l 

4 *\ 

z (*> 

3 > 

Z l /> 


les quatre solutions fondamentales de cliacune de ces quatre equa- 
tions. De ce que les differences entre les a, ... sont supposees 
dtre des nombres entiers, il r£sulte ceci, ainsi que le montrerait 
aisdment Fetude des equations (19), c’est que, si x et y ddcrivent 
un contour ferme, 

rO) Z 0) *(0 Z U) 

subiront La m&me transformation lineaire que 

Z U Zt L^ z Zi Z k' 

Done les determinants contenus dans la matrice 

z t Z Z z k 

z ^ -s‘ 2 l) .. 



8 ) 
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seront multiplies par un meme facteur. Les rapports de ces deter- 
minants sent done des fonetions uniformes et Ton verrait aisement 
(puisque nous n'aurions pas de point singulier transeendant ) que 
ces fonetions uniformes sont rationnelles. Nos determinants sont 
done proportionnels a des fonetions 

R, Rj, R 2 , R 3 , R 4 

qui sont des polynomes entiers en x et en y; de sorte que nous 
aurons la relation 

( *22 ) R-',-rR,^U R*2^1 2 ~ R 3 -l 35 -+- R4^ 4) = 0 , 

ce qui montre la possibility d'exprimer lineairement toutes nos 
fonetions a Paide de quatre d’entre elles. 

On arriverait facilement au meme resultat en partant des inte- 
grates (21), ee qui permettrait de former effectivement ces rela- 
tions ( 22). Je reviendrai plus loin sur ce point. 

II semble d'abord que ce resultat, important pour la theorie des 
series hypergeometriques les plus generates, soit sans inter£t dans 
le cas particular qui nous occupe; cas ou ces series se reduisenl 
a des polynomes entiers. Mais il faut observer que ces polynomes 
peuvent £tre de degre tres eleve, car le nombre m du n° 259 d’oii 
depend ce degr^ peut etre tres grand; tandis que les degres des 
polynomes R de la relation (22) restent limites si les differences 
entre les a (ou entre les { 3 , les y, les v') restent des entiers limites, 
surtout si ces differences sont toutes egales a o ou a ± 1; e’est ce 
qui permettrait d’etablir entre les polynomes hypergeometriques 
de M. Appeli un systeme de relations de recurrence qui en facili- 
terail considerablement le calcul. 

Ces relations seraient analogues aux relations de recurrence 
etablies par Gauss entre ses series hypergeometriques et que Ton 
trouvera a la page i 3 o du Tome III de ses OEuvres completes , 
editees par la Societe de Gottingen. 

265 . Quels sont les points singuliers de la serie de M. Appeli, 
consideree comme fonction de x et de y. On peut les trouver de 
deux manteres differentes, en partant des equations (19). Je sup- 
pose que Ton ait difterentie ces deux Equations par rapport aux 
deux variables x et y; onobtiendra quatre equations, otifigureront 
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lineairement les quatre derivees troisiemes de Les coefficients 
de ces quatre derivees dans ces quatre equations seront 


X — x- 

— 2 xy 

-r 2 

0 . 

0 

X— x 2 

— 2 XV 

— r'b 

— x- 

— ixy 

y—y- 

0 , 

0 

— X’ 2 

— 2 xy 

y —y- 


dont le determinant est egal a 

or’ 2 v- ( x 2 -r -y 2 ■+- 1 — ixy — 2X ~ a y l 

Les points singuliers nous sont alors donnes par x = o, y = o 
et 

x- -I- y- — - i — 'ixy — 2X — 2 y — o, 

c’est-a-dire 

(23 ) \/x i/y = i . 

On pourrait arriver au memeresultaten partantdebintegrale(2i); 
les points singuliers de la fonction definie par cette integrale 
s’obtiendront en envisagent les quatre facteurs de la fonction sous 

ie signe j 

w, p, u — uv — x , v — uv — y. 

Si Ton regarde un instant a et pcommedes coordonnees rectan- 
gulaires et qu’on egaie ces facteurs a z^ro, on obtiendra les 
equations de deux droites et celies de deux hyperboles. II y aura 
un point singulier, si les deux hyperboles se touchent ou si trois 
de ces facteurs s’annulent a la fois. 

Quand la variable x tourne autour de zero, deux des solutions 
particulieres des Equations (19) ne changent pas et deux autres 
sont multipliees par un facteur constant; il en est de meme quand 
la variable y tourne autour de z^ro. Quand x eiy tournent autour 
d’un point singulier satisfaisant a liquation ( 23 ), il y a trois 
solutions qui ne changent pas et une qui est multiple par un 
facteur constant. On v^rifierait aisement ces resultats etron"dt$ter- 
minerait ces facteurs par IMtude des Equations (19). 

Dans le cas particulier qui nous occupe, la s^rie de M. Appell 
se r^duit b un polynome, elle n’a done aucun point singulier et 
les points singuliers que nous venons de determiner appartiennent 
seulement aux solutions des Equations (19). 
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D'autre part, dans ces cas particulars, nous avons 

> . J , . , J 

x — u 2 = cos + -7 v = v- = sin 4 - } 

“ \> 2 

d'ou 

sj x -4- \! y = i. 

Nous retombons done justement sur Fun des points singuliers 
definis par Fequation ( 23 ). Cela iFa pas d’inconvenient puisque 
ces points singuliers, comme je viens de le dire, n'appartiennent 
pas a notre polynome, mais seulement aux autres solutions des 
equations ( 19). Nous verrons meine bientot qu’il resulte de la une 
simplification. 

266 . Jusqu'ici nous avons raisonne comme si p. etv etaient deux 
variables independantes. Nous savons qu’elles sont li6es par la 
relation 

fJL -4- V = I 

de sorte que le polynome P, qui etait entier en jjl- etv 2 , devientun 
polynome entier en v seulement. M. Appell a transforme les equa- 
tions (19) en y remplacant y par v 2 et x par (1 — v) 2 ; ii montre 
ainsi ( Journal de LiouviUe , 1 884 ), par un calcul que je ne repro- 
duirai pas, que la solution generate z de ces equations (19) et par 
consequent P satisfait a une equation du troisieme ordre. II forme 
cette equation (p. 4 1 8 ) et Fecrit : 

(24) (v — v 2 ) 2 >s /// ~h (v — v 2 ) (a-4-6v) z u -1- (c -4- dv -+- ev 2 )z' -y- ( Iv -4- p)z =0, 

z " 7 ... sont les derives successives de 5 par rapport av; les 
constantes a 7 6, ... ont pour valeurs 

a — A — y -t- 2 Y — -4- 3 k -4- 2 -+- s, 
b — — 2 A — y — Y = — 3 (\h -+- k -h 1) — is, 
c = (2y — 0 (A — *y) = (l/c-h 1) (k + «), 

d*= — 2 (aB + 2 Ay' — y') = — 4B-M(A+-i) ^A-hX-t-s — 

e = 4B-i-(2A — 1) (y -h y') 4B + 2 (h-hk-hs — i) (h k - 1-2), 

^ == 4 B (y ^ — 1) = 4 B (A 4- ,/c h~ i), 

P = 2B(i — iy') = — 2B(2£-4- I); 

A = a -h p -4- 1 , B = a $), 
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Comment se fait-il que nos equations (19), qui admettaient quatre 
solutions independantes, se trouvent transformees an une equation 
du troisieme ordre qui nen admet plus que trois? Cela est aise a 
expliquer. En effet. la relation p. + v = 1 equivaut a la relation ( 23 ) 
qui definit un des points singuliers. Quand on tourne autour de 
ce point singulier, trois solutions restent holomorpb.es, tandis que 
la quatrieme est multipliee par un facteur constant, elle devient 
done nulle ou infinie (suivant les valeurs de a, etc.) au point sin- 
gulier; done, quand on fait ;jl = 1 — v, Tune des quatre solutions 
disparait et il n’en reste plus que trois. 

267 . Dans le cas partieulier de s = il se produit une grande 

simplification, comme l’avait montre M. Tisserand et comme l’a 
retrouve depuis M. Appell par un calcul que nous allons resumer. 
Designons toujours, en effet, par les lettres p, q , r, s, t les derivees 
premieres et secondes de ^ par rapport a x et a y, et faisons 

# = fx 2 = (i — v) 2 , r = v2 > 

il viendra 


K ' zS) j d*z 

f 5^ = 4 ( r H a — -tv*)-h‘l(p-irq). 

En Eliminant r, s et t entre ces deux equations et les deux Equa- 
tions (19), et se rappelant que p. -f- v = 1 , on trouve 

(25 bis) v (1 — v) + (Av + B) ~ +C z = 2 D (/)[x4*^v) 
ou A, B ; C, D sont des constantes, et ou 


enremplagant a, { 3 , y, y f par leurs valeurs 

/-f- 5, — g , h -+- 1, k H- 1, 

et nous rappelant que 

/— <r = h+k, 
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nous trou\erons 



Done, pour jf = - ■» noire ineonnue z satisfait a 1 equation difle- 

rentielle lineaire du second ordre cjue 1 on obtient en egalant a 
zero le premier membre de (so bis ). On reconnaLt la foime des 
equations auxquelles satisfont les series hypergeometriques de 
Gauss a une variable. 

Done dans le cas de s = noire polynome P se reduit d une 

serie hyper geo me trique de Gauss ct une seule variable . II arrive 
alors qu 7 une seconde solution des equations (19) s’annule iden- 
tiquement pour u-{-v= 1. L’ elimination de r, s : t entre les 
equations (19) et ( 20) presente quelques particularity. Si Ton 
ajoute les equations (19) et seconde equation (26) respectixement 
multipliees par 



on verra disparaitre a la fois les termes en r 1 s et t si 1 on suppose 

fJt-f- v = t: 

il nous res tera done apres Telimination non pas une, mais deux 
equations dont la combinaison fournira (20 bis). 

268 . Supposons maintenant 5= 1. II semble d'abord que cela 
n’ait pas d 7 inter£t puisque s doit toujours £tre la moiti£ d’un 
nombre impair. Mais on verra bientdt que la simplification que 
nous allons obtenir aura une application importanle. 

Dans le cas de $— 1, le polynome P devient, pour p = 1 — v, 
non plus un polynome hyperg^om^trique de Gauss a une variable, 

ainsi que cela arrivait pour s = mais le carre d’un pared poly- 
nome. 

En effet : 

La s 6 rie hyperg^o me trique de Gauss satisfait k une Equation 
differentielle du deuxieme ordre, son carr^ satisfait done k une 
Equation du troisieme ordre qu’il est ais 6 de former. Si nous cher- 
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chons a identifier cette equation avec Fequation (24}, nousverrons 
que Fidentification est possible dans le cas de s — 1. 

Ce resultat du d’abord a Tisserand, a ete demontre par Stieltjes 
d’une facon fort elegante. Reprenons les equations (19) etfaisons-y 
un changement de variables en posant 

^ = 0 — p) (* — ?')> 7 = pp'; 


nous verrons que, dans le cas de s= 1, les equations (19) se 
ramenent a deux equations lineaires du deuxieme ordre, dont l’une 

ne contient que p, z, tandis quePautre ne contiendra que 

o', 5, “p? O n passe d’ailleurs d’une de ces equations a Pautre 

en accentuant partout la lettre p. 

La premiere de ces equations est line equation difierentielle 
lineaire ordinaire entre z et p, et Fon voit facilement que c’esl 
precisement celle qui definit une serie hypergeometrique ordinaire 
de Gauss. Si done s = F(p) est une solution de cette equation, on 
satisfera au systeme des deux equations enfaisant 


: = F(p) F(p'). 


On conclura de la que notre polynome hypergeometrique a deux 
variables est a un facteur constant pres le produit de deux poly- 
noines hypergeometriques a une variable, Fun en p, l’autre en p'. 
Faisons maintenant 

JJL V — I , 

e’est-a-dire 

p = p'; 

d’ou * 

jk = p 2 = v 2 , a? = (i — P)*= H-*; p = V, 

il viendra 

= [ F 0)P- 


On voit que, quand on tiendra compte de la relation 


P + v = I, 


notre polynome hypergeometrique a deux variables deviendra 
d un facteur constant pres le carre d’un polynome hypergio- 
metrique d une variable e/iv. c. q. f. d. 

P. - II. 6 



82 


CHAPITRE XVIII. 


Je renverraij pour le detail du calcul, au Traite de Mecaniqite 
celeste de Tisserand, t. I, p. 4^8 et suivantes. 


269. Pour comprendre Implication possible du resultat pre- 
cedent, reportons-nous au commencement du Chapitre; soit 

F ~s ~ = (i — a a cos a* -+- a-)” 5 . 

Nous avons trouve 

F“^ — P a /rt ( — %f l w>k p/i 

P etant le polynome hypergeometrique a deux variables que nous 
venons d’^tudier. Dans le cas de s = i , nous venons de voir que ce 
polynome est a un facteur constant pres le carre d’un polynome 
hypergeometrique a une variable. Mais dans ce cas on a 

i _ 1 1“ E* g 1 

h S ~ (I — aE zcr ) (I — aE-'ff) " E ia — E - l< * [ i — i— aE-teJ 


11 est aise de developper le dernier membre suivant les puissances 
de a et Ton trouve 



sin (m n- i) a* 
sin <t 


Done le polynome P n’est autre chose que le coefficient de 
( — z n)/* (’ — w v) A 

dans 

sin (m -+- i)a 
sin a 

en supposant 

2 coscr — ~h z- 1 ) + v(w+ w ~ 1 ). 


Ce coefficient s’exprime done a l’aide des polynomes de Gauss 
(polynomes hyperg^omdtriques a une variable). 

Supposons maintenants quelconqueelrevenons k la formule (i), 
elle peut s’^crire 


(26) + 2 


*</■> 


sin (/c h- i)cr 




sin (k — i)a 


?mcr 


sill <r 



LES POLYSOMES DE TISSERAND. 


83 


Or nous venons de voir que les expressions de la forme 

sin (k — {— i )<7 
sin cr 

pouvaient se developper a l’aide des polynomes de Gauss; il en 
est done de meme de 


270. On peut encore operer d'une autre maniere. Ecrivons 
F-- ? = ( r -b a 2 )-*(i — 2 {3 cosa)-*, 


ou 


il reste a developper 


8 = — 

P i -4- a 2 

(i — 2)3 COS ff)-*. 


Il est evident que nous obtiendrons ce developpement en nous 
reportant aux formules ( 2 ) et (3), en y remplacant a par j3, et en 
n’y conservant que les termes ou e = o. 

Or, ces termes sei'ont ceux ou 

m = h -+- k -r* p -b q, 

e’est-a-dire ceux ou Fexposant de j3 (qui remplace a) est la somme 
deceuxde |j.etdev; il suffira donede conserver, dans lepolynome P, 
les termes du degre le plus 6lev^ en u et v. 

Nous avons trouve plus haut 


(27) Pa"*( — ]xz) h ( — vw) A 

et nous trouvons maintenant 


(28) 


F-* = (1 + p ° f**) A ( — 


V w) k 


P 0 (Slant ce que devient P quand on n'y conserve que les terznes 
du degr6 le plus elev^ en p. et en v; or, nous avons, a un facteur 
constant pres, 


P 


V J jtl. lPt9)(b,p + q) 
^{0,p) (fi\ q) ( 1 , p)(l, q) ^ ’ 


a, b , c, c’ 4tant des constantes et b entier n^gatif. Il faut conserver 
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les termes du degre le plus eleve, c’est-a-dire faire 


d’ou 

Le facteur 


p + q-= — b\ 


(b,p + q) =±(-b)\ 


(a, p -h q) = (a,— b) 


est egalement constant. 

(c',q) = (c', —b — p) = (—l)l J 

= ( — O' 1 


r(i -c') 
r ( I — c' -+- b -+■ p ) 

rd - c 1 ) 

H i — c' 4- b ) (i — d -h 6, p) 


de sorte que (e', q ) est en raison inverse de 

(i — c'+b :P ). 


De meme (i, q) est en raison inverse de (6,/?); nous avons 
done, a im facteur constant pres, 


■ 2 b V (i -c r -+-b,p)(b,p) 



c’est-a-dire que P 0 est egal a y~~~ b multiplid par un poljnome hyper- 

u . 2 J 

geometrique de Gauss en ~ = cot 4 -• 

Si nous rapprochons les relations (26), (27) et (28), nous pou- 
vons resumer nos resultats comme il suit : 

Nous voulons developper A~ 2s et nous cherchons le coefficient 
d’tme exponentielie quelconque 


Cette exponentielie peut etre toujours mise sous la forme 

g/i wk — 1 £ t Hh% + hf[) ■ 

ce coefficient 6tant divisible par v* je Fappellerai 

M {j/* v* ; 

M depend de a, de p. et de v. Je puis le developper : 

i° Suivantles puissances de a, le coefficient de a" 1 est alors un 
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poljnome hjpergeometrique d'Appell a deux variables en cos 4 

sin 4 ^7 qui se reduit a un poljnome de Gauss en sin 2 ^pour s = ~ 

2" Suivant les coefficients de Laplace b s k qui sont desfonctions 
de a; le coefficient de b ! s A) est alors la difference des carres de deux 

poljnomes de Gauss en sin 2 
3° Suivant les quantites 

a /« 

(1 -f- a + s ’ 

le coefficient de chacune de ces quantites est alors egal a 

s i n ^ O/i — ^ — A* 5 — 

•2 

multiplie par un poljnome de Gauss en cot 4 ^ • 

Dans les trois cas le calcul des poljnomes de Gauss peut etre 
simplifie par les relations de recurrence dontj’ai parle plus hautet 
que Gauss a demontrees dans le Tome III de ses OEuvres completes > 
p. 1 3o et suivantes. 
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27 L Nous allons maintenant supposer que les excentricites ne 
sont pas nulles et nous proposer de developper la partie principale 
de cette fonction perturba trice suivant les puissances de ces excen- 
tricites. La meilleure methode est celle de Newcomb qui repose 
sur l’emploi de certains operateurs ( Astronomical Papers, vol. 3). 

Pour bien faire comprendre l’esprit de cette methode. je vais 
prendre d’abord un exemple simple. Soit F(a?) une fonction quel- 
conque de x ; changeons x en 

x -+- ai h -+- a 2 A 2 -+* a 3 A 3 = x -+- 1 

et cherchons a developper F(x -+- oi i h -f- cu A 2 + a 3 A 3 ) suivant les 
puissances croissantes de A; je n’aurai pour cela qu’a appliquer la 
formule de Taylor 

F( " ?+e)= 25 F, " ! '’ 

a remplacer par (a, A + a 2 A 2 + a 3 A 3 )™, a developper cette 
expression et a ordonner par rapport a A. Le coefficient de h k 
seraun polynome de la forme 

b \ F h o F w — h . . . — j— b k F ? 

ou les b sont des coefficients constants et les F w les derives 
successives de F; ceci peut s’ecrire 

(^>iD b% D 2 -i- . . 6 / c D^)F, 

ou D k est Tindice d’une differentiation repetee k fois. Cette 
expression 

b 1 D -+- b/ c D^‘ 
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sera ce qu’on appelle un operateur ; c’esl un symbole cooperation. 
Si j’appelle 11^- cet operateur, je pourrai ecrire 

(i) F(#-^s) = F(#)-f-An 1 F-i-A 2 n2F -r- 


Cela nous indique deja Pesprit de la metkode; mais nous allons 
voir comment on peut former ces operateurs. A. cel effet, je fais 

F(a?) = E^: 

d’ou 

F(x -h s) = EA*E* e , 

U k Ete = E**( b \ X -4- . . .-+- A&X^*). 

En faisant done F(x) = F lx dans la formule (i), nous trouvons 


E^=y A*n* 


en convenant de remplacer T> m par \ m et de faire II 0 = i . 

En d’autres termes, on obtiendra Foperateur ET* en develop- 
pant E D£ suivant les puissances de A comme si D £tait une quantity 
ordinaire et en conservant le coefficient de A*. 

Nous pourrions done regarder E D£ comme un operateur et 
ecrire 

(a) FO-t-s) = E D£ (F). 

Envisageons maintenant une fonction de deux variables F(#, j ) 
et cherchons a d^velopper F(x + e, y -+- e r ) avec 

s = otj A — h 0C2 A 2 — f— 1X3 A 3 — . . . , 
e' = A h- ag A 5 -h 0C3 A 3 -4- 

Nous pourrohs operer de la meme mani&re, nous developperons 
F(^ + e, suivant les puissances de e et de e' par la formule 

de Taylor; nous developperons ensuite e m e fn suivant les puissances 
de A et nous ordonnerons par rapport k A, nous trouverons ainsi 

(ibis) F(a? -b e, y H- e r ) A A n*F, 

ou 11* sera un operateur qui prendra la forme d’un polynome entier 
en D et en D', avec cette convention que 


DpD'? 
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est le signe d'une operation qui consiste a diflerentier p fois par 
rapport a x et q fois par rapport a y. 

Pour determiner ce polynome XI*, il suffira de prendre pour 
F (x,y) la fonction particuliere ce qui, par une analyse 

toute pareille a celle qui precede, donnera 

EAS-rlLE' — }i k H/j 

en remplacant D et D' par A et a dans le second membre. 

Nous pourrons done ecrire symboliquement 

(2 bis) F(x-^£ 1 f^£ , ) = E^^'(F) 

en regardant E D£+D ' £ ' comme un operateur. 

272. On pourrait supposer qne s et s' dependent, non d’une 
seule variable A, mais de deux variables A et A' ou d’un plus grand 
nombre, et qu’on veut developper F^ + SjK + s') suivant les 
puissances de A et de A'; il n’y aurait rien a changer a ce qui 
precede. 

On retro uverait la formule (a bis) sans aucune modification; 
quant a la formule (i bis) il faudrait Pecrire comme il suit : 

F (x -f- s, y -b e' ) = ^ h} til U/j F. 

Examinons de plus pres l’operateur U/y, je dis qu’il y aura un 
cas ou Pon aura 

xifj = n /0 tTo/j 

de sorte que notre operateur a double indice pourra etre regarde 
comme le produit de deux operateurs a simple indice. G’est le cas 
ou s est egal a une fonction de A seulement, plus une fonction 
de hi seulement, et ou il en est de meme de s'. Supposons, en effet, 

£ = £i 4- s 2 , 

. e' = £j 4- s 2 , 

s, et s'* dependant de A seulement, tandis que s 2 et s^ dependent 
de A' seulement. On aura alors 


EDsh-D's'—. ED£ 1 -f-D'5iEDs 2 -+-D'£^ 
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de telle fagon que le coefficient de h £ h f j dans le developpement 
deE Ds+I>,£ ', c’est-a-dire Il/y, sera le produit dn coefficient de h l dans 
le developpement de E D£r4_I> % c’est-a-dire de LI ; 0 par le coefficient 
de h’J dans le developpement de E D£3_hl) ' £ s c ? est-a-dire par FI 0 y. 

Ce resultat s’etendrait evidemment au cas ou le nombre des 
variables serait plus grand. 

Remarquons que d’apres leur definition meme nos operateurs 
sont per mu tables, puisque les operateurs simples D et D 7 le sont 
eux-memes. 


1273. Appliquons ceci au developpement de 

Nous prendrons pour origine des longitudes dans les plans des 
deux orbites la ligne commune des noeuds; la longitude moyenne 
est egale a l’anomalie moyenne plus la longitude du perihelie 
(qui ici se reduit a la distance du perihelie au noeud). Mais, au lieu 
de prendre pour variables la longitude moyenne et la longitude du 
perihelie, ou bien encore l’anomalie moyenne et la longitude du 
perihelie comme on le fait d’ordinaire, nous pouvons evidemment 
prendre la longitude moyenne et l’anomalie moyenne. 

Nous consid^rerons done ~ comme une fonction de 9 variables 

qui seront, outre l’inclinaison mutuelle des orbites : i° les loga- 
rithmes des grands axes; les excentricites ; 3° les longitudes 
moyennes; 4 ° les anomalies moyennes. 

Nous aurons done, apr£s developpement, 


(3) 



e m Q'tn' ^apt^pX'^rql^-q'n 


ou m : tn r , />, p' : q , q 7 sont des entiers, ou l et V sont les anomalies 
moyennes, £ et £ 7 les longitudes moyennes^ et ou A est une fonc-^ 
tion de Finclinaison et des logarithmes des grands axes. 

Si les excentricites dtaient nulles, les longitudes moyennes se 
confondraient avec les longitudes vraies, les grands axes avec les 
rayons vecteurs; et la distance ne depenefrait plus des anomalies 
moyennes l et t! (mais seulement des grands axes et de £, ^)- 
Done dans notre developpement (3) disparaitront tous les termes 
dependant de e, e 7 , /, / 7 , e’est-a-dire tous ceux ou les'entiers /n, 
q, q f ne sont pas nuls. 
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Mais, dans le cas des excentricites nulles, le developpement (3) 
est connu, c’est celui que nous avons obtenu dans le Chapitre pre- 
cedent; et il s'agit d’en deduire Je developpement dans le cas 
general. 

Or il est clair que A ne peut dependre que des longitudes vraies 
et des rayons vecteurs. On obtxendra done le developpement dans 
le cas general, en partant du developpement pour les excentricites 
nulles et en y remplacant les logarithmes des grands axes 

loga, loga' 

et les longitudes moyennes 

par les logarithmes des rayons vecteurs 

r' 

log/’ = log a -+* log > logr' = loga'-t- log 
etles longitudes vraies 

= £ + (o'- £')• 

Il suffit d’appliquer les principes du numero precedent, en 
faisantjouer a j le role deF(#, y) ) a 

loga, log a', C, 5', 

le rdle de x et y, k 

logr, logr', v, p', 
celui de x + e, y + £ ; , et par consequent a 

(4) 1 lo g^, o-s, 

celui de £ et s'. En effet, ces quatre quantity (4) sont d^velop- 
pables suivant Ie§ puissances de 

eE", eE~ lV , e # E"\ e r E~ £r 

qui vont jouer le r61e de h et h! et ne dependent d’ailleurs pas 
de a, a\ 'Q. 

Soit done F le developpement pour les excentricites nulles et 
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r\ le developpement dans le cas general; soient 

D, D lf D', d; 

des indices de differentiation se rapportant a 
loga, £, loga', C; 

nous aurons pour l’expression analogue a De-fD's' 

Dlog^ +D, (*-?)-*- D'JogjH-Di(p'-C'), 


de sorte qu’on aura symboliquement par la formule (2 bis) : 


(5j 


F 1 = 



D' 




>F. 


L’expression 

(6) 



D' 




peut se developper suivant les puissances de 

eE ±il , e'E ±il \ 

en faisant le calcul comme si D, D', D , , D' etaient des quantiles 
ordinaires; les coefficients du developpement sont alors des poly- 
nornes entiers par rapport a D, D', D<, D' . Soit alors 

(e W )« ( e E-")P ( e' E* 7 ' )“' ( e' E-* 7 ')?' 


le developpement ainsi obtenu; ce queje puis ecrire : 


2 


We m e' m ' WW+Q'M 


en ^crivant m, m! , q : q' au lieu de a (3, a'-f- (3', a — (3, a — (3'. 
Le coefficient II sera un operateur qui depend des quatre nombres 
m, m! , q , q\ ce que nous mettrons en Evidence en ^crivant 

t -tmm' 

U <W" 

Nous aurons alors 

F, = y ene'n'EW+tfnnffiF. 



9 2 


CIIAPITKE XIK. 


Maintenant F c'est le developpement dans le cas des excentri- 
cites nulles, tel quhl a ete obtenu an Chapitre precedent. Soit 


2 




ce developpement. Nous sommes conduits a calculer 

n^'fAE 

pour appliquer cet operateur a la fonction AE 1 (pZ+pV) ilfautfaire 
subir a cette fonction diverses differentiations par rapport a 

loga, 1 oga', C; 

mais, pour la differentier par rapport a £ ou a il suffit evidem- 
ment de la multiplier par ip ou ip 1 . Nous pouvons done Pecrire 

n^’fAE w+p'O] = E^/^t)n;^'(A), 

Dans le premier membre II™'?' a le meine sens c[ue plus haut; 
e’est un polynome entier par rapport a D, D', D 1? D' ; dans 
le second membre, e’est le meme polynome, mais oii les sym- 
boles D(, D' f sont remplaces par les quantites ip et ip ' . II ne con- 
tient done plus que deux symboles d’operation D et D' et il est 
applique a la fonction A qui ne depend plus que de loga et log a'. 
Il reste finalement 


P- = Fi = "V e'w' A. 

A JmmU ( W 


Le coefficient du terme en 

e m e ' m r £ i ( p £ 4 -// £ '■ ■+■ q /-(- (/ 7 ’ ) 


est done 


TT/n/n' \ 

A V/ A > 


A designant le coefficient du terme en 


lequel, ne dependant pas des excentricites, a At 6 d^termin^ au 
Chapitre precedent. 

274. Nous avons Vu au n° 240 qu’ii est ais6 de passer du 
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developpement procedant suivant les anomalies excentriques an 
developpement procedant suivant les anomalies moyennes. Nous 
pouvons done nous proposer, au lieu de former directement ce 
dernier developpement, d’y parvenir indirectement et par Pinter- 
mediaire du premier. G’est ce qu’a fait M. Newcomb. 

A cet effet nous prendrons comme variables Pinclinaison mu- 
tuelle des orbites el, pour chaque planete, le logarithms du grand 
axe, Pexcentricite, Panomalie excentrique u et ce qu’on pourrait 
appeler la longitude excentrique , c 5 est-a-dire Panomalie excen- 
trique, plus la longitude du perihelie, toujours comptee a partir 
du noeud. Nous chercherons alors a former le developpement 


(3 bis) 



m g in' Jgz(/>Y]-h. p'ri -‘rqu-v-q’ it’) 


oii u et u f designent les deux anomalies excentriques, 7 ] et r/ les 
deux longitudes excentriques. 


27o. M. Newcomb introduit encore une autre simplification en 
posant 


1 -+- e' 2 


II cherche ensuite a developper, non plus suivant les puissances 
de e et e ! , mais suivant celles de e et e' de fagon a obtenir le deve- 
loppement 


(3 ter ) 


1 

A 


2 


QUprrhp'-q'+qu+q'u') % 


11 est clair qu’il est ais 6 de passer du developpement (3 ter ) 
developpement (3 bis ) et de la au developpement (3). 

Si nous posons 


il vient 


e = sm<p, 


au 


cp 

s = tang A. 


Remarquons en passant que les elements canoniques et r\ { 
du n° 08 sont proportionnels a sin |* 


276. 11 y a des terxnes des developpements (3), (3 bis ) et 
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(3 ter) qui sont connus; ce sont ceux ou 

m =r m = q — q' = o, 

les seuls qui subsistent quand les excentricites sont nulles; nous 
les avons determines au Chapitre precedent. Ils sont d’ailleurs les 
memes pour les trois developpements, sauf la substitution de r t etr/ 
a s et 

De ces termes connus, il s’agit de deduire les autres; et, pour 
cela, que nous adoptions les variables du n° 274 ou celles du 
n° 275, il nous suffira d’operer tout a fait de la meme maniere 
qu’au n° 273; on n’aura qu’a faire jouer a 

e , e', 7] , r/, w, it' 


le role que jouaient au n° 273 les quantites 

* c, r, l'- 


On trouvera comme au n° 273 que le coefficient de 


ou celui de 
sera 


qui Q r m' g/( /> r,-t- p ' r t '-4-/7 « -4-7 'u') 

rttnm' \ 
u r/r/' 


A designant le coefficient du terme en 


tandis que II™™' est un operateur qui n’est autre chose que le coef- 
ficient de 


ou de 


Qtn q’ in' j£i(qu-hq f u') 


£in ^ m' u r ) 


dans le d^veloppement de Fexpression symbolique 

(7) (0 D (?) DE ‘ 7,,P “ r ‘ ,+, ' //< ^'’’ 

suivant les puissances de eE~‘ u , e'E~ iu ' ou celles de eE ±1 “, s'Er'*'. 
On remarqueraque, dans l’expression ( 7 ), j’ai remplac< 5 D, etD', 
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par ip et ip ainsi que nous avons demontre qn 5 on devaitle fairc 
a la fin du n° 273. 

Seulement les operateurs II sont beaucoup plus simples quand 
on adopte les variables du n° 274- et surtout quand on adopte 
celles du n° 27o. 


277. II j a d’abord une remarque que nous devons faire et qui 
restera vraie que l’on adopte les variables des n ,,s 273, 274 ou 27o. 
L’expression symbolique (6) est le produit de deux autres 

( 8 ) 


II en est de merne de Pexpression (7), a la condition bien 
entendu de changer £ et 'Q en r\ et 7/; quant a D, et D' ils peuvent 
etre remplaces par ip et ip’ } ainsi que nous Pavons dit. 

La premiere des deux expressions (8) ne depend que de eE 1 ' 7 , 
ou de eE ±lu , ou de s E ±lu : la seconde au contraire ne depend que 
de e r E ±l1 ' ou de e r E ±lu \ ou de b ! E~ IU ' . 

Appliquant en consequence une remarque faite au n° 272, nous 
voyons que Poperateur II a quatre indices peut etre considere 
comme le produit de deux operateurs a deux indices 

(9) n-r = n-°n 0 7 . 


Le premier facteur depend seulement de D et de D*, c’est-a-dire 
de D et de /?; le second de D'etdeD', c’est-a-dire de D 1 et de p l . 
Mais de plus nous avons une relation entre D et Db En efiet, 


et de meme 


DF = 


dF 

d loga 


dF 

a ~r~ 

da 


JD'F = a 


t dF 
da’ ’ 


Mais la fonction l - est homogene et de degrd — 1 en a et en a! . 


Celte propria n’est alter^e ni par Pop^ration D, ni par Pop^ra- 
tion D', elle appartiendra done a toutes les fonctions auxquelles 
nous aurons ik appliquer soit Poperateur D, soit Poperateur D'; 
de sorte qu’on aura 


dF , dF 
a da + a da! 
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D + D' = — i . 

Cela nous permet d'exprimer D' en fonction de D de sorte que 
nos operateurs II seront des poljnomes entiers par rapport a Pope- 
rateur simple D et par rapport aux entiers p et p ! . 


278. On rencontre de plusgrandes simplifications encore quand 
on adopte les variables du n° 27o. On a alors, en effet, 


d*ou 


1' I — 2£ cos w -+- 

— =i — e cos u = 

a i -+- s 2 


^ j D = (i -f- e*)- D (, — sE*“) D (i — eE-^) d . 


Reste le facteur 


E°i( 


Onobservera que, en vertu del’equation des aires, 




d’autre part 
d’ou 


r 2 dv 
a 2 dl 


= /i-e=, 


du 


i — e cos u = — j 
a 


dv _ \J \ — e 2 ___ i — £- _ i zE~ £u 

du ~~ i — e cos u ~~ i — 2£ cos u-h z- ~ i — e E iu i — eE-* w ’ 

dv 


— l & mj£miu _j_ 

et en integrant 
On voit que 




mi 

Ei>i ( ^ — *n ) 


> E 77z g-w/a 

d mi 


se decomposera en deux facteurs, le premier dependant seulement 

' r \ D 


de e E m comme le deuxi&me facteur de (- j , le deuxi£me depen- 
dant seulement de sE"~*“ comme le troisieme facteur de Mais 
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il vaut mieux operer comme il suit. Considerons 1‘expression 

- E iv . 
a 


Si la longitude du perihelie est, pour un instant, supposee nulle, 
elle est egale a 


cos u — e -i - i sin usj \ — e~ = 


K lu \ i - 


s- E - £u ) 


ou a 


E /Z£ , „ . 

— — (i — s h~ lu )‘ 2 . 


On aura done, quelle que soit la longitude du perihelie, 


- = 

a 



Reprenons l’expression symbolique 

E D i i ^ i = (^~ j ; E 



elle pourra s’ecrire 


(i -h £ 2 )"’ D (l — £E 2M ) D -^(l — sE" tup-bp. 


Or, on a par la formule du binome 


(i_6E*«)»-p = ^ 
(i — s E- f «) 1> -*-/ > = ^ 


(D — *-M, k)_,„ 
kl 

( D — p — 1, k’) 

k’\ 

(DH-p — F-4 -i ? F) 
k"\ 


( — z) k ' 


(— e)** 


la notation (D, k) ayant m£me sens qu’au Chapitre precedent. 
Ces formules nous donnent imm^diatement nos op^rateurs II. Si, 
en effet, nous tenons compte seulement d’abord des deux derniers 
facteurs et que nous posions 

(i — e E*“) D+ P(i — 


P. - II. 


7 



il \iendra 
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, _ m -+- a m -h q \ / n 

D— /> — l —2, — — 1 




H™ = (— n«- 


et ensuite 


n'"» = H"* 


/» — c/ l m-q l 
2 " 2 


1 r ) fJrn—2 , ( ^ 1 7 2 ) JJW-4 

71 n 7 + o' V 


on en mettant en evidence le terme general 

ttoto _ X? |j/«-aA* 
n 70 — 7^ /-t 11 </ 

On voit que toas nos operateurs se presentent sous la forme de 
poljnomes entiers en D. 

M. Chessin prefere, au lieu de passer par les intermediaires des 
developpements (3 ter) et (3 bis), former directement les opera- 
teurs qui donnent le developpement (3). Ces operateurs sont 
encore des poljnomes en D, mais dont les coefficients sont beau- 
coups plus compliques; mais M. Chessin a donne dans V Astrono- 
mical Journal des relations de recurrence qui en facilitent beau- 
coup le calcul. 


279. Une fois les operateurs formes, il est ais£ d’obtenir les 
developpements eux-memes. Au Chapitre prt§c£dent, nous avons 

developpe i en supposant les excentricites nulles et cela sous plu- 
sieurs formes differentes. Nous avons trouv6 d’abord 


l = V E am 

A a' m ~ hi 


E HpZ+pV), 


B etant 6gal a un polynome de Gauss en v = sin-^? multiple par 
un facteur num^rique et par des puissances de p = cos 2 i et 

v==sin 2 -- C’esl la formule (27) du Chapitre prudent. 11 s’agit 
alors de voir ce que c’estque 


D 


a m 

a r m-\-i 7 


D' 


a m 

a ’ m -h ’ 
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on trouve 


a m 

a ' 9 


— an — 


a m 


II suffit done de remplacer D, D', D, D' par m, — (m-f-i), 
ip , ip ' ; notre expression (6) deiient tout simplement 

/• in f W ' 

-I (-; ^ipw-’O^-ip-w-D 

a 1 \ a j 


■et le developpement (3 ) est ramene au suivant : 


i=y B ii 

A Jmmd r WH-1 




formule d’ailleurs evidente sans tant de detours. 

Nous avons trouve aussi une autre forme de developpement 


i = V B -L b\ k E iW+p’W, 
A a y 5 


P+P' P~P' 

•ou B est ^gal a p 2 v 2 , multiplie par la difference des carres de 
deux polynomes de Gauss. C J est la formule (26) du Chapitre pre- 
cedent. Cette fois Pemploi des op^rateurs sejustifie. II faut appli- 
quer nos opera teurs II a 


•et pour cela ii faut savoir appliquer a cette expression l’op^rateur 
simple D une ou plusieurs fois; or, on trouve 


•et les expressions 



d b < 

T_ 

dx ? 


a 


bf 


.s’exprimeraient de m£me tr£s simplement a Paide des d^riv^es 
sup&rieures du coefficient de Laplace Or nous avons appris 

au n° 250 a calculer par des relations de recurrence les d£riv£es 
successives des coefficients de Laplace. Le probl&me peut done 
^tre consid&r^ comme enti&rement rdsolu. 
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CONVERGENCE DES SERIES. 


i80. Nous a\ons vu au n" 240 que la partie principale de la 
fonction perturbatrice *- pouvait se developper, soit sous la forme 

soit sous la forme 

N A w E 


et au n° 342 que les coefficients de ces deux, series pouvaient s’ex- 
primer par des integrales definies 


(*) 

_ C C dr civ 

i\7Z B /n/tl ' — ^ / , y— 

J J x ni y ni v/R(.r, V) 

et 


( 2 ) 

r r Q dr dy 

4 TC - A.ffim' == / / , , > 

J J x m y ,n 

ou 

x = E /w , y — E* w ', R = .r 2 j^A 2 , 




2Q 


— me ( x ) ~r - m' e' ( y — — 

\ * ) V r 


et ou Pintegration doit etre dtendue, tant par rapport k x que par 
rapport a y 7 le long d’une circonference ayant pour centre l’ori- 
gine et pour rayon Punite. 

Ces coefficients B mm , sont des fonctions de Pinclinaison, 
des excentricites, des longitudes des p^rihedies complies a partir 
du noeud et enfin des grands axes. Nous avons appris dans les deux 
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Chapitres precedents a de\elopper ces fonctions suivant les puis- 
sances des excentricites et des inclinaisons et il s'agit maintenant 
de savoir quelles sont les conditions de convergence de ces series. 

Nous naurons pour cela qu a appliquer la methode de Cauchy 
et a chercher les points singuliers de ces fonctions. 

281. Nous sommes done ainsi amenes a expliquer comment on 
trouve les points singuliers des fonctions representees par des 
integrates definies et d’abord par des integrales definies simples. 
Soit 

j R (x,z)dx, 

une integrale definie prise par rapport a x le long d’un certain 
contour : cette integrale sera alors fonction du parametre z. 

Pour que pour cette fonction une valeur de s soit critique, il 
faut d’abord que Tun des points singuliers de R(hr, z) 1 consideree 
comme fonction de x , se trouve sur le contour d’integration. 

Mais, comme on peul deformer ce contour d‘une maniere con- 
tinue, on peut le faire fuir de\ant le point singulier, et I’on n’est 
arrete que lorsque ce contour se trouve pris entre deux points 
singuliers et ne peut plus fuir. 

On obtienclra done toutes les valeurs critiques de 5 , en expri- 
mant que deux des points singuliers de R, consideree comme 
fonction de x , se confondent. Mais toutes les valeurs critiques 
ainsi trouv^es ne conviennent pas; il faut, en effet, que les deux 
points singuliers qui se confondent ainsi soient, avant de s’etre 
confondus, de part et d’autre du contour; e’est seulement a cette 
condition que le contour, pris entre deux feux, ne peut plus fuir. 

Soit done 

<p(.r, z) = o 

liquation qui exprime que la fonction R(#, z) presente une sin- 
gularity, et supposons qu’elle se decompose en un certain nombre 
d’equations ind^pendantes, trois par exemple : 

<PjO, 3)= 0, © 2 (#* s) = o, cp 3 (ar. 3 ) = o. 

On obtiendra les valeurs critiques de 5 de l’une des deux ma- 
nures suivantes : 

i° En annulant deux des trois fonctions <p 2 , ©t en ^cri- 
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vant, par exemple, 

.r, z) = 'f 2 * x, "*i = o. 

En eliminant et resolvant par rapport a 3, on aura une valeur 
critique de 

2 ° Eel annulant Fane de ces trois functions et sa derivee et ecri- 
vant, par exemple, 

, doi 

On aura encore une valeur critique de s, en eliminant x et resol- 
vant par rapport a 

282. Considerons maintenant une integrate double 


R( 5 t,j, z) dxdy , 


envisagee comme fonction du parametre ^ et etendue a un champ 
quelconque. 

Le cas particulier le plus simple est celui ou Ton doit integrer 
par rapport a x, le long d’un contour fixe independant dejK, et 
par rapport a v, le long d’un contour fixe independant de x . C’est 
precisement ce cas particulier simple que Fon rencontre en ce 
qui concerne les integrates (i) et ( 2 ). 

Ici, en effet, les deux contours C x et C r ind^pendants, le premier 
de v , le second de x, sont deux eirconterenees de centre z^ro et 
de rayon 1 , le premier dans le plan des x, le second dans le plan 

desjK* 

Soient done C* et C r les deux contours fixes d’integration par 
rapport a x et a y. Soit 

d(y,z)— J R (x,y,z)dx, 

Fintegrale £tant prise le long de C#. 

Notre integrate double sera alors 

nU) = J 9 (/i z ) d r : 

Fintegrale 6tant prise le long de C r . 
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Nous n avons plus qu a appliquer les principes du numero pre- 
cedent aux deux integrales simples 


J R dx, Jo dy. 


Nous devons observer toutefois que nous risquons de n'obtenir 
ainsi que des conditions necessaires et pas suffisantes. En effet, 
nous avons observe plus baut que toutes les valeurs critiques ne 
conviennent pas. 

On obtiendra d’autres conditions egalement necessaires en chan- 
geant de coordonnees, et en particulier en permutant x et y, c’est- 
a-dire en intervertissant l’ordre des integrations. 

Soit 

= o 

Tequation qui exprime que la fonction R a une singularity 
Decomposons cette equation en equations irreductibles 

($) <?i(x,y,*) = o 9 = o, y*(v,y , -) = o. 

On obtiendra les singular! tes de la fonction 0 (y, z) : 

i° En annulant deux des fonctions cpo, cp 3 et ecrivant, par 
exemple, 

= o; 

2 ° En annulant une des trois fonctions et sa deriv^e el ecrivant 



o. 


Ges singularity peuvent ne pas toutes convenir, car il peut arriver 
que les deux points singuliers, qui se confondent, soient d’un 
meme cot6 du contour d’int^gration. 

Pour avoir les singularity de r i (z) f considerons maintenant 
celles de 9 (y, s), qui nous seront donates par des syst£mes d’equa- 
tions de Tune des deux formes 


? 1 = ¥- 2 = 0 , 




ffoi 

doc 


o, 


et cherchons les conditions pour que deyix de ces singularity se 
confondent. 
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Si nous considerons ^ comme un parametre, x et y comme les 
coordonnees d"un point dans un plan, les equations ( 3 ) repre- 
sentent un rertain nombre de courbes planes. 

Les singularities de H donnees par un systeme d’equations de la 
forme 

©1 = ©2 = O 

correspondront aux points d'intersection de ces courbes ; celles 
qui seront donnees par un systeme de la forme 


correspondront aux points de contact d’une de ces courbes avec 
une tangente parallele a Paxe des y. 

Pour que deux de ces singularites se confondent, il faut : 

i° Ou bien que trois des fonctions © s’annulent a la fois 

cp 4 = cp 2 = ©3 = 0 ; 

2 0 Ou bien que Ton ait a la fois 

c/cpi 

©1 = =0. ©0 = 0; 

1 ax * “ 


mais, comme la condition ne doit pas dependre du choix des coor- 
donnees, et qu’en particulier elJe doit subsister quand on per- 
mute xety, on devra avoir aussi 


et, par consequent, 


do x 

©1 = -7— = O. Cpo = O 

‘ dx ' T “ 




__ _ dfcpi __ ^ 


dx dy 


ce qui veut dire que Pune des courbes ( 3 ) aura un point double; 
3 ° Ou bien que les deux courbes 


©1 = 0, Cp2 = o 

soient tangentes Pune a Pautre ; 

4 ° Ou enfin que les deux courbes 
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soient tangentes Tune a l'autre. \lais cela entraine : ou bien 


c hl 

tlv 


= o. 


ou bien = o; mais, comme ia condition ne doit pas changer 
quand on permute x et r, on de\ra avoir dans tous les cas 


Cl = 


clz > ! 
~dx 


d r 3\ 

dr 


= o. 


En resume, les valeurs critiques de ^ sont : 

i° Celles pour lesqtielles trois des courbes (3) se coupent en 
un meme point; 

2 ° Celles pour lesquelles deux de ces courbes sont tangentes; 

3° Celles pour lesquelles une de ces courbes a un point double. 

Seulement toutes ces valeurs peuvent ne pas convenir. car les 
deux singularites qui se confondent peuvent etre situees d un 
meme cote du contour d'integration. 

Et d’ailleurs nous a\ons vu que nos conditions ne sont que neces- 
saires. 


283. Appliquons ces princ.ipes aux integrates (0 et ( 2 ). La 
fonction sous le signe / sera holomorphe tant par rapport a x et y 

que par rapport aux excentricites e et e' et a sin-> J etant Tincli- 
naison. 

11 j aura exception seulement : 

i° Si e = 1 et e'—i (condition oil x et y n’interviennent pas 
et sur laquelle nous reviendrons); 

2 0 Si x ou y est nul ou infini ; 

3° Si A s’annule, c’est-a-dire si 

R = A a? 2 jk 2 , 

qui est un polynome entier du sixteme ordre en x , y, est 6gal 
a z^ro. 

Cela est vrai quels que soient les entiers m et m! ; cela est vrai 
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d ailleurs aussi bien de l'integrale \ 2 ) que de 1 integrate ( i). car 

Q et E 11 

ne cessent d etre holomorphes que si x ou y est nul ou mfim . 

Les courbes qui correspondent aux courbes (3 ), c est-a-dire 

celles dontles equations expriment que ia fonction sous le signe ^ 

cesse d'etre holomorphe, se reduisent done, en ce qui concerne 
les integrales (i) et ( 2 ), aux quatre droites 

(4) ^ = o, y — o, X = =c : y = n 

et a la courbe du sixieme degre 
(4 bis) R = o. 

Cette derniere, dans le cas ou 1 inclinaison est nulle, se decom- 
pose en deux courbes du troisieme degre 

(4 ter) Ri=o, FU = °* 

Pour trouver les valeurs critiques des exeentricites ou des incli- 
naisons, nous devons done cbercher celles pour lesquelles trois 
des courbes (4) ou (4 bis) se coupent, ou pour lesquelles deux 
de ces courbes se toucheni, ou pour lesquelles une de ces courbes 
a un point double. 

Mais, comme ces courbes sont les memes pour les intdgrales (i) 
et(a), les valeurs pour lesquelles une de ces trois circonstances 
se presentera seront egalement les memes pour les integrales (x) 
et ( 2 ). 

Ainsi les valeurs critiques des excentricitds ou des inclinaisons 
sont les memes pour les deux integrates (i) et ( 2 ). 

Mais une question pourrait encore se poser; nous avons vu que 
toutes les valeurs critiques ne conviennent pas. Ne pourrait-il se 
faire qu’une de ces valeurs convint & ( 1 ) sans convenir k ( 2 ) ou 
inversement. 

La r^ponse doit £tre negative. Comment se fait-ii, en efFet, que 
certaines valeurs critiques conviennent et que d’autres ne con- 
viennent pas? Pour nous en rendre compte, Faisons varier d’une 
manure continue Tun de nos param&tres, par exemple l’excentri- 
cite e, et, en m£me temps, d^formons d’une mani&re continue les 
contours d’intdgration. Nous devons diriger cette deformation de 
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telle sorte quaucune des valeurs singulieres definies par les equa- 
tions (4; et (4 bis \ ne se trouve dans le champ dlntegration. 
Si, e variant d une maniere continue depuis zero jusqu'a e 0 , on 
peut s arranger de facon que cette condition ne cesse jamais d ? etre 
remplie, c est que e 0 n’est pas une veritable valeur critique; e 0 ne 
convient pa s-; si, au contraire, il est impossible de s’arranger pour 
que la condition soit remplie (parce que, corame je Texpliquais 
plus haut, le contour d’integration se trouve prisentre deux points 
singuliers), c’est que e 0 est une veritable valeur critique; e 0 con- 
vient. 

Faisons done varier e de zero a e 0 et, en meme temps, deformons 
nos contours d’integration en partant des contours initiaux 

[ jr I = 1 

qui sont les memes pour les integrales ( 1 ) et (2); si e 0 ne con- 
vient pas a Tintegrale (i), c’est que nous pouvons deformer nos 
contours de telle facon qu’a aucun moment une des valeurs sin- 
gulieres satisfaisant aux equations (4) et (4 bis) ne se trouve dans 
le champ dhntegration. Mais, si cette condition ne cesse jamais 
d'etre remplie en ce qui conoerne Tintegrale ( 1 ), elle ne cessera 
jamais non plus de l’etre en ce qui concerne Tintegrale (2), 
puisque les equations (4) et (4 bis) qui definissent les valeurs sin- 
gulieres sont les memes pour ( 1 ) et pour (2). Done e 0 ne con- 
viendra pas non plus a Tintdgrale (2). c. q . f . d . 

Void done un premier r^sultat. 

Les coefficients du developpement de la fonction per - 

turbatrice suivant les anomalies excentriques , ainsi que les 
coefficients du developpement suivant les anomalies 

moyennes sont eux-memes develop pables suivant les puissances 
des excentricites et des inclinaisons . 

Les limiles de convergence de ces nouveaux developpements 
sont les m&mes pour les coefficients et pour les coeffi- 

cients A m'm'l etles sont les memes pour tous ces coefficients 
quels que soient les entiers m et m 1 . 

284. Nous examinerons sp^cialement le cas oh les excentricites 
sont nulles et celui ou Tinclinaison est nulle. 
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Dans le cas ou les exeentri cites sont nulles, ii n y a plus a faire 
de distinction entre les anomalies excentriques et les anomalies 
moyennes, de sorte que 


A. //> m ' B/j 


Q = I, 

-ass 


a = O, 


dx d v 


yin y It ( X. y ) 


Cette derniere integrale peut encore s’ecrire 



d.r dr 

sp/tl-hl y III’ ~hl ^ 


ou bien en posant comme au n° 259 



nous pomons encore Pecrire 


1 

2 



dz dw 


1 r r dz dw 

2 J J zp wp’ A * 


en posant 



p'= 


m -4- in' -+- 9 
2 


Quand x ety decrivenl chacun dans leur plan des circonferences 
de rayon i, z et cc font de meme ; mais, a chaque couple de valeurs 
de £ et w situees sur ces cteconferences, correspondent deux 
couples de valeurs de x et de y\ on obtient done deux fois Pintd- 
grale, de sorte que, finalement, nous pourrons ecrire 


(5) 


— 4^ 2 A 



dz dw 
zP Wp' A ? 


Pinlegration etant dtendue, tant pours que pour (p, a une circon- 
faience de rayon i ayant son centre a Porigine, de telle fa§on que 

M = M =I - 

II sera plus commode d’opdrer sur Pintdgrale (5) que sur Pint6- 
grate (i). Nous avons alors 


X- — ct ~ H— ~t~ [xckz' 
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Nos courbes singulieres sont alors 

z = w = o, z = ZC, w = DC, A 2 = o. 

Nous devons alors chercher Ja condition pour que trois de ces 
courbes se coupent en an meme point, ou pour que deux d'entre 
elles se louchent, ou pour que l'une ait un point double. 

Nous n : avons pas a nous inquieter de la premiere condition, la 
courbe l' 1 — o passe to uj ours par Forigine; or, pourqu’une func- 
tion definie par une integrate definie presente un point singulier, 

il faut que deux points singuliers de la fonction sous le signe j* 

primitive merit separes viennent a se confondre; ctest*a-dire, par 
exemple, que trois courbes singulieres qui primitive merit ne pcis- 
sciient pas par un meme point viennent a passer par un meme 
point. Ce n’est pas ici le cas, puisque les trois courbes A 2 = o, 
z = o, w = o passent toujours par un meme point. 

De meme la courbe A 2 = o passe toujours par les points 

,3 = 0, w = cc; ,3 = oc, W — o ; Z = w = 00 , 

de sorte que nous n’avons pas a nous inquieter de cette condition. 
Pour que : = o ou : = x touche A 2 = o, il faut que 

fx aa! — o soit u = o. 

Pour que w = o, ou w = 00 touche A 2 = 0, il faut que 
vaa' = o soit v = o. 

0 ait un point double, il faut 

dtf __ dA 2 ^ 
dw ~~ dz °’ 

W = ±Z I, Z — ± I 

a 2 -h a '- dz 2 1 xaa'zt 2 vaa' = o. 

285. Examinons de plus pr£s ces r^sultats. Nous verroas d’abord 
que les solutions p = o ou v = o ne pouvaient convenir a la ques- 
tion. Supposons, en effet, d’abord que, remplagant p par i— v, 
on d^veloppe suivant les puissances dq v. Si la valeur v= o <itait 


Pour que A 2 = 

c’est-a-dire 

et 

( 6 ) 
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un veritable point singulier pour la branche de la fonction consi- 
deree, le developpement suivant les puissances de v serait toujours 
divergent; or nous savons qidil n*en est pas ainsi. 

Si la valeur a = o ou v == i etait un point singulier, nous 
demons en conclure que la serie diverge des que |v|>i, mais 
nous n'aurions pas a nous inquieter de cette consequence puisque v 
est toujours <C i . 

Supposons maintenant que Ton developpe suivant les puissances 
de p etv regardees comme des \ariables independantes ; si p= o, 
ou bien v = o etait un veritable point singulier pour la branche 
de fonction consideree, le developpement serait toujours diver- 
gent, et nous savons qu’il n’en est pas ainsi et que le developpe- 
ment converge pourvu que a et v soient assez petits. 

II est aise d’ailleurs de constater que, pour les petites valeurs 
de p et v, les points singuliers qui se confondent pour v = o, par 
exemple, sont tous deux interieurs, ou Lous deux exterieurs au 
contour dbntegration. On pourrait toutefois se demander s’il en 
serait encore de meme pour des valeurs plus grandes de p. Mais 
la reponse est aisee; si les deux points singuliers sont tous deux 
interieurs au contour pour p petit, puis que nous fassions varier p 
d’une maniere continue, les deux points singuliers A et B et le 
contour varieront aussi d’une maniere continue, mais ces deux 
points ne pourront cesser d’etre interieurs au contour que si l’un 
d’eux, A par exemple, franchitle contour; mais alors, nous pour- 
rons toujours fair e fuir le contour devant le point A, a moins que 
ce point A ne se confonde avec un autre point singulier C primi- 
tivement exterieur au contour, different par consequent du point B. 
Mais alors la convergence cesserait deja par suite de la coincidence 
des points A et C. On n’a done pas a s’inqui£ter de la coincidence 
possible des points A et B puisqu’elle ne pourrait amener une 
singularity de notre fonction, qtfapres que la convergence aurait 
deja cesse pour d’autres causes. 

Ainsi les points singuliers p =. o, v = o ne conviennent pas k la 
branche de fonction consideree; nous devons ajouter qu’ils pour- 
raient convenir k d’autres branches de la fonction qui seraient 
definies par Fint^grale (5) ytendue a d’autres contours que ceiui 
qui est d^fini par les deux Equations 
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286. Pour la branche etudiee. nous n'avons done a nous pre- 
occuper que de la condition i 6 L 

Posons v = sin 2 ^ et proposons-nous de developper suivant les 
puissances croissantes de v. La condition ( 5 ) de\ient 
a 2 -i- a’- -±- 2 ad [ dz < i — v)±v] = o. 
ce qui peut s ? ecrire de I‘une des deux manieres suivantes 
(5 bis) < a'zh a)-~ o ; a 2 -r- d-dz iaa\ i — av) = o. 

Dans la premiere de ces relations v n'entre pas ; la seconde 
donne 

(a'dza )- 
4 ad 


Le rayon de convergence sera done la plus petite des deux 
valeurs 

‘ ( ddz a)- ] 

4 ad 


c 7 est-a-dire 


(y~ a? 

4 cia 


Ainsi la condition de convergence du developpement sera 



(d— a ) 2 
kad 


La discussion des equations (5 bis) montrerait de meme que les 
coefficients sont developpables suivant les puissances de a, 

pourvu que 

a < d . 


Considerons maintenant p et v comme des variables ind^pen- 
dantes et d^veloppons suivant les puissances de a etv. 

La relation (6), qui definit les limites de convergence, s’6crit 
alors 

a 2 d- dz 2 ad jji d= a ad v = o 
ou 

, a 2 d- 


La condition de convergence du developpement est done 


lnl-HvK 


a*-h a ' 2 


%ad 



II 2 


CHAPlTRh XX. 


Cette condition est toujours remplie, car on a 


I = p; i v I 


j u . — 


= i < 


a 2 a ' 2 
a act 


jYos coefficients sont done ton jours deceloppcibles suivant les 

/ .,J - >J 

puissances de cos- ~ et sin- -• 


Nous avons donne, au n° 270 , un developpement qui procede 
suivant les puissances de 


\m va 

i — a - 9 i -i- a 2 5 


c est le developpement ^ 28 ) du Chapitre XVI1L D'apres ce qui pre- 
cede, ce developpement converge pourvu que 


pa 


va 

I -+■ 0 C“ 

— r 

1 a 2 


<1. 


Or, comme a, v et a sont essentiellement positifs et plus petits 

a , 1 • I 

que i, que ■ est, par consequent, plus petit que - et que 
jjl — f— v = 1 , on aura 



va 

1 -f- a 2 ] 

1 -r a 2 


( P *+■ v ) - 


< 


La condition de convergence sera done toujours remplie. 

Du developpement (28) du Chapitre VIII on passerait aisement 
au developpement (27) qui procede suivant les puissances de a, 
p. etv, il suffirait de developper chaque terme en 

/ [jloc ^ m j va \ » 

\ 1 a 2 J [i + a 2 / 


suivant les puissances de a, ce qui est possible pourvu que a «< 1 . 
Comme cetLe condition est toujours remplie, le developpement (27) 
converge toujours. 


287 . Consid^rons maintenanl le cas particulier ou rinclinaison 
est nulle et d^veloppons suivant les puissances des excentri- 
cit^s e et 
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Les courbes qui correspondent aux courbes < 3 i sont les 
courbes i 4 et < 4 ter i, a savoir 

x = o. y = *K x = y = ac, Rj = o. R, = o. 

i° Cherehons d abord la condition pour que R* = o ait un point 
double. Pour nous en rendre compte, eberchons la signification 
de ces equations. Soient C Porbite de la premiere planete, C celle 
de la seconde; d'apres notre hypothese ces deux coniques sont 
dans un meme plan. A chaque valeur de x correspond un point M 
de la conique C et a chaque valeur de y un point M' de la 
conique C'. 

Inequation R = o exprime alors que la distance MM' est nulle, 
ce qui peut ne pas vouloir dire que les deux points M et M' coin- 
cident puisque ces deux points peuvent £tre imaginaires. 

L’equation R ( == o exprime que le coefficient angulaire de la 
droite MM' est egal a — i et Pequation R 2 = o exprime que ce 
coefficient est egal a — y — i . 

Les equations x = o, x = zo expriment que le point M est & 
Pinfini sur Pune des deux branches infinies de la conique C; les 
equations y = o, y — co expriment que le point M' est a Pinfini 
sur la conique C'. 

Pour que la courbe R i = o ait un point double, il faut et il 
suffit que la droite MM' soit tangente a la fois aux deux coniques C 
et C'. Les deux coniques C et C' ont un foyer commun, le Soleil. 
La droite MM', qui est une tangente isotrope a C, doit passer par 
un des foyers de C, et pour la inline raison elle doit passer par un 
des foyers de C'. Soient alors S le foyer commun de'C et C', R le 
second foyer de C, R' celui de C'. 

Pour que les deux coniques admettent une tangente isotrope 
commune (en dehors de celle qui passe par S qui existe toujours 
et qui ne saurait convenir), il faut done que la droite MM' passe 
par R et R'; e’est-a-dire que la droite RR' ait pour coefficient 
angulaire \J — i . Cette condition s’exprime analytiquement comme 
il suit : 

ae E-^= a'e'E~ £ **'. 

Je repr^sente par w et m r les longitudes des p^rih^lies. 

C’est Ik la condition pour que = o ait un point double (en 
— II. 8 


P. 
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dehors de celui qui, correspondant a la tangente isotrope passant 
par S. existe toujours et ne saurait comenir). 

2 ° De ineme la condition pour que FU= o ait un point double 
s’ecrit 

ae E^= a'e' E 1 ® 1 . 


3° Pour que les deux courbes R 1 = o, Ro — o soient tangente s, 
il faut que les deux coniques C et O se touchent; car les points 
d’intersection a distance finie des deux courbes R i = o, Ro=o 
correspondent aux quatre points d’intersection des deux coniques. 

11 faut done que la distance des deux foyers R et R/ soit egale 
a la somme ou a la difference des grands axes; ce qui s’ecrit 
(aeEto—a'e'E**)(aeE-to— a’ e’ E-*®') = (a r ±«) 2 . 

4° II faut voir ensuite si Tun des points d’intersection des 
courbes (4 ter) ne peut pas se confondre avec l’un des points 
d’intersection de l’une de ces courbes avec l’une des droites (4) 
ou de deux de ces droites entre elles. 


Ces deux categories de points correspondent respectivement : 
aux quatre intersections de C et de C ; (M et M 7 confondus) et aux 
cas ou les deux points M et M ; sont a l’infini sur les deux coniques C 
et a. 

La condition est done que 1’une des quatre intersections de C 
et de CJ s’dloigne indefiniment, e’est-a-dire qu’une asymptote de C 
soit parallele a une asymptote de C'. 

Cela s’ecrit d’une des quatre manieres 


tang 2 |E 2 ^= tang 2 E 2 **', 
cot 2 2 = tang 2 ~ E 2 *^', 


en posant 


tang 2 2 E-*fcr = tang 2 2- E -2 ^', 
cot 2 2 E“ 2/CJ ~ tang 2 2- E~ 2 **', 
e = sin cp, e’ = sincp'. 


Je puis £crire aussi 


l rfc sj i — e~ 


l zh sj i — e'~ 
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en prenant dans chaque fraction le signe inferieur dans les deux 
termes ou le signe superieur dans les deux termes. 

\ oici done quelles sont les equations qui definissent les valeurs 
critiques 


(7 ) ae E-^= ci 

<8j aeE ir7 = a'e'E 1 ^'. 


( 9 ) 
{10 ) 


( ae E /gt — a' e'E l &' j { ae E~ /Gr — a 1 e E~^' ) = ( a f ±: a ) 2 , 


1 dz \l 1 — e - 
1 =F /i — 


E 2 *® = 


1 ±: \f 1 — e'- 
1 =F V / i — e 7 * 


E-* CT ', 


auxquelles il convient d J adjoindre 

{1 1 ) e = 1, e'= 1. 


Mais nous avons vu que toutes les valeurs ainsi trouvees peuvent 
ne pas convenir et il est merne aise de voir qu’elles ne peuvent pas 
toutes convenir. 

En effet, envisageons Fequation (7); elle est satisfaite pour 

e = e' = 0. 


Si done les valeurs critiques definies par cette equation conve- 
naient, notre fonction presenterait des singularites pour des 
valeurs tres petites dee ete'. Nos coefficients ne seraient pas deve- 
loppables suivant les puissances de e et de e\ meme pour des 
valeurs tres petites de ces quantites . Nous savons qu’il n’en est 
pas ainsi. Done les valeurs definies par liquation (7) ne sauraient 
convenir. 

Le meme raisonnement s’applique aux equations (8) et (10). Il 
importe de remarquer que Fequation (10), malgre la presence des 
doubles signes, ne repr^sente qu’une seule equation analytique 
irr^ductible qui prendrait sa forme algebrique si Ton chassait les 
radicaux. Cette Equation irreductible est satisfaite pour e = e 1 = 0. 

Restent les Equations (9) et (1 1). Pour trouver la limite de con- 
vergence d^finie par l’equation (9), remarquons que a, o', m et w 1 
sont des donndes de la question et sont r^elles, mais que e et e ! 
qui sont nos variables independantes pourront prendre des valeurs 
imaginaires. 

Donnons k e une s^rie de valeurs de module constant \ e |, mais 
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d argument variable; iaisons de meme pour d . Le maximum du 
module du premier membre de { 9 i est 

a- \ e- | -r- a' 2 | e' 2 1 2 ad | ee' cos(bj — ®')|; 

la condition de convergence est done 

a 2 1 e 2 | -f- a' 1 | e 2 \ <iaa’ \ee co$(m — ts’ )\ <( a' — ct) 2 . 

Nous sommes amsi conduits a supposer que les seules condi- 
tions de convergence de notre dei eloppement sont 

j |«l<i. 

{l2) j a i\ e i \-x- a '*\e 2 1 -+- iaa’\ee’ cosf t«j — zct' ) [ < (a — a ) 2 . 

Je ne voudrais pas entrer dans trop de details 5 je ne puis cepen- 
dant me contenter d un apercu : il faut done que j explique en 
quelques mots comment on peut donner a la demonstration toute 
sa rigueur. 

Considerons le domaine D defini par les inegalites (12). II ne 
contient pas de valeurs singulieres satisfaisant aux equations (9) 
et (1 1). Mais il en contient qui satisfont aux equations (7)5 (^ )? ( 1 o). 
Si Tune de ces valeurs convenait, il en serait de meme de toutes 
celles qui satisferaient a la meme equation et qu'on pourrait ren- 
contrer en faisant varier e et e 1 d’une maniere continue et sans 
cesser de satisfaire a cette equation. Cela serait vrai au moins en 
ce qui concerne la determination de la fonction que l’on attein- 
drait par cette variation continue. 

Or on verrait qu’on peut alteindre, par une variation continue, 
une quelconque des valeurs singulieres qui satisfont aux Equa- 
tions (7), (8), (10) et aux inegalites (12) en partant de la valeur 
e = o , e'—o et sans cesser de satisfaire a ces equations el sans 
sortir du domaine D. La valeur singuliEre e — 0 : e ! =o ne eon- 
venanl pas, aucune de ces valeurs ne convient. 

Les conditions (12) sont done les scales conditions de co mer- 
gence. 

La troisiEme condition (12) sera satisfaite quels que soient gj 
et iri, si l’on a 

| ae | ( 4- | aV | < a’ — <2, 
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c’est-a-dire si la distance perihelie de Tune des planetes est plus 
grande que la distance aphelie de 1'autre. 


288. Le theoreme que nous avons applique a la fin du numero 
precedent est analogue a celui que nous avons rencontre au n" 28o 
et pourrait s’etablir de la meme maniere. Mais il vaut mieux le 
rattacher a un theoreme d’ Analyse. 

Si, dans un domaine D simplement connexe, une fonction de 
deux variables F(x : y) est holomorphe, sauf peut-etre pour les 
points qui satisfont a une relation analytique 

^ = ?(y)> 

et si un de ces points 

y=y o, v = *(yo) 

est ordinaire, il en sera de meme de tous les autres. 

En effet, d’abord si ce point est ordinaire, il en sera de meme 
de tous les points voisins. 

Decrivons maintenant dans le plan des x un contour tres petit, 
autour du point douteux x = o(jy). Soit x = ‘i(y) un des points 
du contour. Nous pouvons supposer que le contour se d^forme 
d’une maniere continue quand y varie d’une maniere continue, et 
que ty(y) est une fonction analytique. 

Je dis d’abord que F(x,y) est une fonction uniforme. 

Supposons, en eflet, qu’elle ne le soit pas, que F \(x,y) el 
F 2 (#, y) soient deux de ses determinations, et que l’on passe dc 
Tune a 1’autre quand x d^crit le contour. Alors 

F i[<Krb y\ — Fd'Krbrl = *(y) 


serait une fonction analytique de y, mais cette fonction serait 
nulle poury=j' 0 et pour les valeurs voisines, puisque la fonc- 
tion F(x,y) est uniforme pour y=y 0 et pour les valeurs voisines^ 
La fonction &(y) est done nulle pour toutes les valeurs de y et, 
par consequent, F (x,y) est uniforme pour toutes les valeurs dc y. 

La fonction F ^tant uniforme, la condition pour que x = ®{y) 
soit un point ordinaire, e’est que 1’int^grale 


/ 


x m F (a?, y) dx , 
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prise le long du contour, soil nulle quel que soit l'entier positif m. 

Cette integrale est une fonction analytique de y\ elJe est nulle 
pour v = Vo et pour les valeurs \oisines puisque x — 'f(yo) est un 
point ordinaire ainsi que les points voisins. Elle est done identi- 
quement nulle, et le point x = 'f(yi est ordinaire quel que soit y. 

Le theoreme s'etendrait aisement au cas ou tous les points de D 
seraient ordinaires, sauf peut-etre ceux qui satisfont a Pane des 
n relations analytiques 

= x = ?-2(v), x = v, L (y), 

et oil l 7 on saurait, d'autre part, que les points 
y—yu ^ = 

r=y2, X = \y 2 ), 

y=yn, x = <? fl (yn) 


sont ordinaires. 

Nous nous bornerons aux cas particuliers simples traites plus 
baut, bien que les memes proced^s soient applicables au cas general. 
Nous mentionnerons toutefois que M. Feraud a traite d’autres cas 
particuliers dans le Tome X des Annates de L° Observatoire de 
Bordeaux . 
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RELATIONS DE RECURRENCE ET EQUATIONS DIFFERENTIELLES. 


289. Reprenons les equations du n 1 ' 280 


(1) 

( 2 ) 


— 4^ 2 B„ 


— 4«*A W 


-SS-. 


dec dv 


x m y m ' )/K(x,y ) 

~ C C doc dy 
J J cc m y m '\/R(x,y) 


qui definissent ies coefficients du developpement, soit suivant les 
anomalies excentriques, soit suivant les anomalies moyennes. 
Rappelons que 


R = &2y2 ^2 ? 


2 Q 


Q= [:_£(,+ !)] [, + 


Ces coefficients A et B sont des fonctions des elements, par 
exemple des grands axes, des excentricites et des inclinaisons et 
l’on peut se proposer Fdtude analytique de ces fonctions. Je me 
propose de montrer : 

i“ Que ces fonctions satisfont a des equations diff^rentielles li- 
n^aires, dont les coefficients sont des fonctions rationnelles des 
elements. 

2 ° Qu’il y a entre ces fonctions, ou du moins entre les B, des 
relations de recurrence lin^aires dont les coefficients sont des 
fonctions rationnelles des elements. 


290. A cet effet, je vais envisager une expression plus generale 
definie comjme les A et les B par une intdgrale double. 
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r r H E— cLt dv 

n = / J — 7TF r 1 - 


une integi ale doable prise le long de contours fermes. 

H, Q et F sont des polynomes entiers en x, —•> r, j dont le 

degre est respectixement k. to et f. Quant a s , c’est la moitie d un 
entier impair. 

Ax ant d’aller plus loin, il taut que j'explique ce quej'entends 
par le degre d'un polynome en y% ~ • Un pareil polynome 

est de la forme 

A jr a y b , 

a et b etant des entiers positifs ou negatifs. Je dirai qu’il est de 
degre m si 

| a j 1 m , !/>(_/??. 

Un polynome de degre m contient (2 m -h i)- coefficients arbi- 
trages, puisque a de meme que b pent prendre 2W + 1 valeurs. 

Geci pose, je suppose que 0 et F soient des polynomes donnes, 
mais que nous fassions varier H. II y a (2 h - 1 - i ) 2 polynomes H de 
degre h lineairement independants, mais pour quelques-uns 
d’entre eux la fonction II est nulle. En elfet, si P est un polynome, 
l’integrale 

/./ iyf^ d:cd} ' 

sera nulle. En effet, int^grons d’abord par rapport a x, nous trou- 
verons 

PE a 

y F* -1 * 

Comme nous integrons le long d'un contour ferine, celte quan- 
tity reprend la inline valeur aux deux limites et Eintegrale est 
nulle. 

On verrait de meme que, si Q est un polynome, I’intdgrale 

fj 

est nulle. II suffirait d’intygrer d’abord par rapport a y. Done 
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n sera nul si 

HE- _ _d_ PEQ _ d_ QE^ 
xy F*' dx dy 


e'est-a-dire si 


( 4 ; 


TT c/P „ d ( 2 d¥ 

H= X -J- F -f- X -7- PF -r- (1 5 ) X P 

dr dx dx 


.dQ 


dQ 


d¥ 




dy 


dy 


0 . 


Si P et Q sont de degre />, nous voyons que le deuxieme membre 
est de degre p f-r- to. Nous devons done avoir 

p = h — f — to. 


Ii y a done (2 A — 2 f — 2 to i ) 2 polynomes P distincts et au- 

tant de polynomes Q. Nous pouvons done former 

2 ( 2 h — 2 f — 2 to -T- 1 ) 2 

relations de la forme ( 4 ). Mais ces relations sont-elles distinctes; 
ne peut-il pas arriver que le deuxieme membre de i 4 ) soil identi- 
quement nul ? Cela arrivera si 

d_ PE- d_ QEQ _ 

dx y¥ Ji ~ l dy irF 5 * -1 °’ 

e’est-a-dire si 

QEQ dx _ PEQ dy 
F*- 1 x F 5 " 1 y 


est une differentielle exacte; je repr^sente cette differentielle par 



dU, 


ce qui donne 

1 

dS _ 

dQ 

x d¥ 



- s)x 7Tx 

( 5 ) 

dS _ 

dQ „ 

x d¥ 

1 


-•"-j 


Je dis que S est unpolynome entier 
que e’est une fonction uniforme. 


1 1 

en x. -? r, - 
J x J y 


et d’abord 
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i° L'integrale f dXJ ne presente pas de periodes cycliques. 

Supposons que I'on calcule cette integrale en regardant y comme 
une constante et en faisant decrire a x dans son plan une courbe 
fermee, a l’interieur de lacjuelle il j ait un nombre pair de va- 
leurs dea; telles que F = o. Il en resultera que quand on aura decrit 
cette courbe fermee \/F et F f_l reviendront a leur valeur initiale. 
L'integrale prise le long de ce contour serait une periode cjclique 

de f dU. 

Cette periode, si elle existe, sera une fonction dey. Je dis que 
cette fonction devra se reduire a une constante. Si en effet la fonc- 
tion U admet deux determinations IF et on aura 


et, par consequent, 
( 6 ) 


dU' _ dU" _ PEQ 
dy dy y F 5-1 ? 


d(W-U') 

dy 


= o. 


Si Ton se reporte main tenant a l’analyse de M. Picard dans son 
Ouvrage sur les fonctions algebriques de deux variables (t. I, 
p. 88 a 90), on verra que cela ne saurait avoir lieu sans que la 
periode cyclique soit nulle si le polynome F est le plus general 
de son degre et dans beaucoup d’autres cas encore. 

2" L’integrale ne presente pas de periode polaire. Elle pourrait 

en presenter si la fonction sous le signe J devenait infinie, mais 

cette fonction ne peut devenir infinie que pour x — o, on y — o, 
ou F = 0. 


Pour F = o, la foncLion est infiniment grande d’ordre s — 1 a 
moins que F ne soit divisible par un carre parfait, ce qui n’arrivera 
pas si F est le polynome le plus general de son degr^. La fonction 
n^tant pas infinie d’ordre»entier, il ne peut r^sulter de Ik une 
periode polaire. 

Soit maintenant x = 0. D^veloppons la fonction sous le signe J 

suivant les puissances enti&res positives et negatives de x de telle 
fagon que Pon ait 
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les C m etant des fonctions de y; alors nous aurons une periode 
polaire qui sera 

a i-C-!. 

En vertu de Tequation (6), cette periode devra etre une con- 
stante independante de y. Je clis qae cette constante est nulle, 

Soxty =y 0 une des racines simples de l’equation 

F(o ,y) = o. 

Bien entendu, avant de faire x = o, il faut multiplier F par une 
puissance convenable de x de fagon que F reste fini pour x = o. 

Faisons variery de facon que cette "variable revienne a sa "valeur 
initiale apres avoir tourne autour de y 0 ; alors y/F se changera 
en — y/F, en — ^ ? et, par consequent, C_ f en — C_ { ; mais, 
comme est une constante, on a 

C— i = — G_!, 

d’ou 

G_! = O. c. Q* F. D. 

II n 7 y a done pas de periode polaire pour x — o et l’on verrait 
de meme qu’il n’y en a pas pour y = o. 

Le raisonnement precedent ne s’appliquerait pas toutefois si F 
se reduisait a un carre parfait pour x = o. 

3° Mais S pourrait encore ne pas etre uniforme pour une autre 
cause. Supposons que Ton prenne Tintegrale U le long d’un con- 
tour enveloppant un nombre impair de valeurs de x qui annu- 
lent F. Alors y/F se change en — y/F quand on parcourt ce con- 
tour. 

Done Use change enU '=h — U, h etantune constante qui nepeut 
dependre de y. Ici liquation (6) ne s’applique pas et la difference 
U — U' n’est pas independante de y; car, y/F ayant change de 
signe, ^ et “^j7 ne sont P as % aux ? ma -is egaux et de signe contraire, 
de sorte que 

dlJ__dU dh_ 

dy ~ dy 3 dy ~~ 

Ainsi h est une constante absolue. De plus, si nous considerons 
un autre contour analogue qui change U en U v = h! — U ; je dis 
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que les deuxeonstantes h et ti sent ideniiques ; sans quoi FintegraL 
admettrait la periode h — h. 

La constante h etant la meme pour tons les contours, nous pon 
vons choisir la constante arbitraire dintegration de facon qm 
h — o. Alors notre integrate ne sera susceptible que de deu: 
valeurs L et — L . D'ailleurs le signe de U changera en menu 


temps que eelui de \ F, de telle sorte que 

est une fonction uniform e. 

4° S est un polynome entier en 


— — et par consequent 


i i 


En efiet, Fintegrale et S ne peuvent de\enir infmis que si le 
derivees 

dl J £U 
dx ’ dr 


de\iennent infinies, e’est-a-dire si 


x = =c. y = ac, x = o, y — o. F = o. 


Pour F = o, les derivees deviennent infinies d’ordre s — i etpa 
consequent L devient inlini d'ordre s — 2 et S reste fini. 

Voyons comment se comporte S pour x tres grand. Les terme 
preponderants de F, Q, Q se reduisent respectivement a x b , x> 
et multiplies par une fonction p et co etant les degr^s d< 

F ? Q, Q). Done ^-reduit a ses termes preponderants se reduira ; 


et L T a 
et enfin S a 


9 (y)xf^- g )" r i } ~ i E'P/^‘ W 


6( y E’FP^ 


73 {y)xi J -f~^, 


0(y) et ?i(y) £tant des fonctions de y faciles a former. 

On trouverait le meme r^sultat en cherchant comment S s< 
comporte pour x — o; il suffirait de changer /, p et co en — f, — / 
et — to. On raisonnerait encore de meme pour y tres grand ou tr& 
petit. 

Done S se comporte comme un polynome de degr6 
p — f — to = h — 'K f — 2 w . 
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Cest done un polynome de degre h — 2 f — 2 to, dependant de 

<2 h — 4 /' — 4 to — 1 p 

coefficients arbitraires. 

Done panni ]es 

21 2 h — if — 2 to — ii 2 

relations ( 4 ) il y en a seulement 

2 1 2 h — 2 / — 2 oj — 1 1 2 — 2/2 — 4 /* — 4 to - 4 - if 2 . 

qui sont distinctes. Done nous n’avons que 

( 2 h -+- 1 )“ -+- (2/2 — 4 f — 4 f » -e ij " 2 — 2 ( 2/2 — 2/ — 2 to ~ 1 » 2 

expressions II qui soient lineairement independantes. Or ce 
nombre est egal a 

8 (f^r- CO f 2 . 

II est done independant de h et s. Done, quelque grand que soit le 
degre de h du polynome H, nous n’aurons toujours, au plus, que 
^(/+ tL> )‘ expressions II distinctes. 

Remarquons que les expressions ou s = s 0 rentrent comme cas 
particuliers dans celles oiix = 5 0 H- 1 , puisque nous pouvons changer 
s en s-\- 1 et H en HF sans changer P expression II. Done toutes 
les expressions II, quel que soit le nombre s, le degre du polynome 
H et ce polynome lui-meme, peuvent s’exprimer lineairement a 
l’aide de 8(/+ co)- d’entre elles. 

Si nous prenons 

8 (/-+- to) 2 1 

expressions II quelconques, il y aura entre elles une relation 
lin^aire, dont les coefficients seront des fonctions rationnelles 
des coefficients des poiynomes F, 0 qui sont les mimes pour 
toutes ces expressions II ainsi que des divers poiynomes H qui 
ne sont pas les mimes pout * les differentes expressions II. Cela 
r^sulte de la fagon dont nos relations (4) ont 6t6 formees. 

291. Le nombre des expressions II distinctes peut s ? abaisser si 
les poiynomes pr<§sententcertaines especes de sym^trie. Supposons 
par exemple qu$ F, & et H ne changent pas quand on change x 
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en — x et y en — y. Si alors nous ecrivons ces polynomes sous la 
forme 

^ A 


nous voyons que les deux entiers a et b doi\ent etre de meme 
parite. 

Un polynome de degre m presentant cette symetrie ne contient 
plus que 

2 m- -h 2 m i 


coefficients arbitrages. Nous n’aurons done que 

2 A 2 -t- 2 h 1 

polynomes H. 

Si les polynomes P et Q presentent cette meme symetrie il en 
sera de meme du second membre de*(4). D’ailleurs nous obtien- 
drons de cette facon tous les polynomes H symetriques qui sont 
de la forme (4). Soit en effet 

*(P, Q) 


le second membre de (4). Considerons deux polynomes P et Q 
non symetriques de telle fa^on que P(#, y) ne soit pas egal a 


et supposons que 


p (— ^ —y) 


soit £gal a un polynome symetrique H de telle sorte que 


nous aurons a la fois 


H (*,y) = #[ p Q (&,y)] 

et 

H(27,jk) = H(—^, — y) = 4>[P(— a-, -y), Q( — — jr)], 
ou, en additionnant et divisant par 2 , 

qui est une relation de la forme (4) ou les polynomes P et Q sont 
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remplaces par les polynomes symetriques 

P ( x, v ) — P < — x , — v\ Q i x . y 1 — Q 1 — ,r, — y ) 

, : 1 

‘A > 

II suffira done de considerer les relations (4 ) formees avec des 
polynomes symetriques. Nous aurons done, non plus 

1 {'lk — 2 / — 20) -r- 1)2 

relations (4j, mais 

2[a( A — f — O) ) 2 -4- 2f A — /— 0) ) -4- t]. 

De meme, si nous nous reportons aux relations ( 5), nous verrons 
que, si P et Q sont symetriques, il en sera de meme de S ; nous 
aurons done 

2 ( A — 2/— 2 to y 1 -+• 2 f h — 2/ — 2 0))-i- 1 

polynomes S. 

En resume le nombre des expressions II distinctes ne sera plus 
(2A -j- I )" ( 2 A — 2/— 20)4- l >2 — 2 ( 2 A — / — CO -+- 1)2 = 8{/-+- CO) 2 , 

mais seulement 

6(A) -i- ft ( A — 2/ — 20)) — 2 0 ( A — / — co), 

en posant 

0 ( h ) =2 A 2 -+- 2 h — (— 1 . 

Ce nombre sera done 

4(/-eo)) 2 . 

292. Appliquons ce qui precede aux expressions (1) et (2), mais 
.auparavant il est preferable de les transformer en posant 



Alors, A 2 , qui etait un polynome contenant des termes en 

1, a? 2 , ar~ 2 , 7 2 , y- 2 , xy , ary-*, ar-iy, x~iy~i, x, y } x~i, y-\ 

aleviendra un polynome entier du second ordre en 



CHAIMTRE XXI. 


Ce polynome ne change pas quand on change ; et r, en — q 
et — t 4 . Ce sera lui qui jouera le role de F ; nous avons done 

/=*• 

Le polynome Q devient de son cote un polynome symetrique da 
premier degre en 


’’ £ 7 T ‘ j ^ 5 


de sorte que 


Dans la formule (i), bien entendu, ou l’exponentielle E a ne 
figure pas, on prendra co — o. 

Le role du polynome H sera joue dans l’expression (i) par x / n yj n r > 
et dans Y expression (2) par "~;- f , qui nous donneront Tune et 
F autre un polynome symetrique en 

1 1 

i, r 1. v 

Nos coefficients li mm ' et A mm * seront represents a un facteur 
numerique pres par Fintegrale 

r r u 


ou Q est du degr^ o ou 1 et ou 


ou enfin A joue le role de F 2 . Ce sont done des expressions de la 
forme II. 

Les coefficients de ces polynomes F, Q et H sont d’ailleurs des 
fonctions rationnelles des grands axes a et a', des excentricits e, 
e l et de y/i — e 2 , \J 1 — e' 2 , ou, si Y on aime mieux, de e et de e', en 
posant 


Ce sont egalement des fonctions rationnelles des lignes trigono-~ 
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metriques dependant de l'inclinaison etdes perihelies, c'est-a-dire 
de tang ^ tang ^ ? tang — 

Ce sont en resume des fonctions rationnelles des sept quantites 

a, a\ s, tang tang 15 , tang— , 

2 ‘2 2 


les longitudes des perihelies etant comptees dans les plans des 
deux orbites, a partir de Intersection mutuelle de ces plans. Ce 
sont ces sept quantites que nous appellerons les elements. 

Envisageons maintenant la derivee de ou de B mm ' par rap- 

port a l’un de ces elements, que je designe par a; elle sera ^gale a 


ou 


H' = 




*•(§ 


■»£) 


Ih 


d 42 
di 


sera un poljnome de m4me forme que H. Nous retrouvons done 
une expression de la forme II avec cette difference que 1’exposant 

I . 3 

s qui etait egal a - est maintenant -• 

Si nous considerons une derivee partielle seconde, cet exposant 
deviendra rien d’ailleurs ne sera change, et il en serait de 
meme pour les d^rivdes d’ordre superieur. 

Ainsi nos coefficients B mm < ainsi que leurs derivees 

partielles d’ordre quelconque par rapport aux elements sent 
des expressions de la forme U. 

Les conclusions du n° 290 sont vraies quand le poljnome F est 
le plus general de son degr 6. Elies subsisteraient bien entendu 
dans les cas particulars, par passage a lalimite; mais ici elles sont 
directement applicables. 

En effet le theor^me de M. Picard s’applique directement toutes 
les fois que le poljnome F est indecomposable. 

C’est ce qui arrive pour le poljnome A 2 dans le cas general du 
probl&me des trois corps. Quand l’inclinaison est nulle, le poljnome 
P. - II. 


9 
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A- se decompose, mais nous restons dans les cas tres etendus ou le 
theoreme de M. Picard est applicable. 

II n'y a done pas de periode cyclique. 

Pour montrer ensuite qifil n ; y a pas de periode polaire, nous 
avons suppose que F(j?, y) ne devenaitpas un carre parfait quand 
apres Favoir multiplie par une puissance convenable de x on y 
fait x — o. 

Si nous appliquons cette regie au polynoine 

Tj j, 

cela reviendra a reduire A- aux tenues en x 2 , rr, y 2 , ou enx~ 2 1 
, J'~ 2 : ou en x 1 , xy~ { , y“ 2 , ou encore en x~ 2 , JK 2 - Ce 

polynome ne se reduit pas ainsi a un carre parfait. La condition 
est done remplie. 


293. Appliquons d’abord cela a l’expression (i) et aux coeffi- 
cients B mm >] nous avons 

f = 2, 0) =0. 

Le nombre des expressions distinctes est done 
4 < /^ u >) 2 = [ 6 . 

Nous n’avons done que seize coefficients distincts. Done : 

Si Von envisage le developpement de la fonction perturba- 
trice suivant les anomalies excentriques, il y aura entre les 
coefficients des relations lineaires de recurrence dont les coeffi- 
cients seront des fo notions ratio nnelles des elements . Ces rela- 
tions permettent d 3 ex primer to us ces coefficients en fo notions 
de seize d 3 entre eux. 


Ces relations de recurrence sont assimilables a celles qui relient 
les coefficients de Laplace. Ces relations peuvent s’obtenir de la 
fagon suivante; ^crivons les deux identity 


d &/ 1 


d I . 


■ — ( L 

dy \ A 


-4- 2 y -r- A 2 = o. 
J dy 
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Remplacons dans ces identites A 2 et ~ par leurs developpements 

suivant les puissances entieres positives et negatives de x et y, et 
egalons a zero dans Tune de ces identites le coefficient de xPy$. 

On a retrouve dans le J\achlass de Jacobi, sans demonstration, 
un enonce d’apres lequel le nombre des coefficients distincts 
serait non pas 16 , mais i5. Cela est possible, car 16 n'est ici qu’un 
maximum. II serait interessant de pousser la verification jusqu’au 
bout. 

294. Supposons maintenant que Ton envisage un coefficient et 
ses derivees partielles des divers ordres par rapport aux elements. 
Ce sont encore des expressions II, elles sont done liees a seize 
d’entre elles par des relations lineaires. 

Ainsi chacan de nos coefficients B, consider e comme fonction 
de V un quelconque des elements , satisfait a une equation dif- 
ferentiate Uneaire du seizieme ordre au plus , dont les coeffi- 
cients sont des fonctions rationnelles des elements . 

Sii’onconsidere deux ouplusieurs elements comme des variables 
independantes, on peut former aussi un tres grand nombre d’equa- 
tions lineaires aux deinvees partielles. 

Tons nos coefficients B, et toutes les derivees partielles peuvent 
s’exprimer lineairement en fonctions de seize de ces coefficients, 
ou bien encore en fonction de Pun d’eux et de quinze de ses d&ri- 
vees partielles. 

293. Passons maintenant au cas de l’expression ( 2 ) et des 
coefficients A mm ' , qui sont ceux du developpement suivant les 
anomalies moyennes. Cette fois 

/=2, CO = I, 

4(/-H <*>) 2 = 36. 

Seulement le polynome Q depend de ni et de m 1 . II n’est done 
pas le m£me pour deux coefficients A mm ' differents. 

Les conclusions du n° 293 ne se g^n^ralisent done pas. 

Au contraire, le polynome Q est le meme pour un coefficieni 
A mm f et pour toutes ses derivees partielles, de sorte que le th^o- 
r&me dji n° 294 se g^n^ralise immediatement. 
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Chacun des coefficients A, cons^dere conime fonction de L un 
qnelconque des elements, satis fait it une equation di fferentielle 
lineaire du trente-sixieme ordre au plus. II satisfait en outre a 
de nombreuses equations lineaires aux derivees partielles, si Ton 
regarde les divers elements conime des variables mdependantes. 


296. Passons au eas ou les excentricites sontnulles; il en resulte 
de grandes simplifications. D’abord. en eftet, il ny a pas de dis- 
tinction a faire entre les anomalies excentriques et moyennes, de 
sorte que Ton a 

0 = 0, A mrn ,=z ^mm' m 

De plus on a 

A 2 = a- -+- a'~ -h act p {^z -4- l - ^ -+■ act'* ( w -l- —J > 

en posant 

3 = j = 'O 2 , w = xy = 

D’ailleurs les coefficients A mm : sont exprimes a un coefficient 
nuinerique pres par Tintegrale 



dz dw 
z w A 


, . II 

ou H est un polynome entier en • 

On a done 

/=!, C0 = 0. 

Le polynome qui joue le role de F est ici A 2 , et Ton voit que r 
pour z = o ? F expression z A se reduit a une constante; le raison- 
nement par lequel on demontre qu’il ne peut y avoir de periodes 
polaires n’est done pas applicable immediatement [puisque nous 
avons du supposer plus haut que le premier membre de liquation 
F(#, 7 ) = o ne se reduisait pas a un carr 6 parfait pour x = o]. 

Nous pourrions nous tirer d’affaire de plusieurs manieres diffe- 
rent es : 

i° En nous servant des fonction elliptiques, comme je Fai fait 
dans le Bulletin astro nomique, t. XIV ; 

2 ° En revenant aux variables x et y. 

Mais je pr^f&re me servir de la symetrie du polynome A 2 . 



RELATIONS DE RECURRENCE ET EQUATIONS DIFFERENTIELLES. 1 33 

Ce polynome ne change pas en eflet quand on change z en 
ou bien w en ~ Un polynome satisfaisant a cette double condi- 
tion sera ce que yappellerai desormais symetrique . 

Dans le developpement de ~ les coefficients de 


z a , z~~ a . Z a W~^, Z~~ a W~~b 


sont egaux, et peuvent etre represents a un meme coefficient 
numerique pres — par Pune des integrales doubles 


ou 



v=b 


dz dw 
zw A 


= f 

J J zw 1 


H 


Nous pouvons done toujours suppose r que le polynome H est 
symetrique . 


Nous devons ensuite rechercher quels sont les polynomes H 
synttriques qui nous conduisent a une integrale double identi- 
tiquement nulle. Ce sont ceux qui sont de la forme (4)- Dans le 
second membre de Pexpression (4), il faut, bien entendu, faire 
Q = o et remplacer x ety par les \ariables nouvelles 5 et pp. 

Je remarque maintenant que si, dans cette expression (4) ainsi 


modifte, on remplace P(s, (*>), Q(#, pp) par — wj, 
Q Pexpression H(s, w) se change en H vvj. De meme, 

si l’on remplace P(a,w), Q(*, vo) par p(z,dj, —Q(z,±), 

1’ expression H (z, w) se change en H (z. ~'j ■ 

Si Pon a alors 

P ( ,, W )=-P (!,«-)= 


Q(^< v)= q(i, 


nous dirons que P et Q ptsentent la sy me trie croisee. 
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II est clair que, si P et Q presentent la symetrie croisee, l’expres- 
sion H presentera la symetrie directe. 

Si H presente la symetrie directe, je dis que I’on pourra tou- 
jours, sans restreindre la generalite, supposer que P et Q ont la 
symetrie croisee. En effet, si H est symetrique, on pourra, sans 
changer H, remplacer P et Q par 



- p (i, «>j 

’ w )> 

ou par 

P (*, — ), 

- Q (z, —), 


' w ) 


ou par 

-I), 




v= ». 


ou enfin par les deux polynomes 


P(-s, w ) — p(l, w ) + p(^,-L)-p(l,i.) 

4 

q (3 , w) + q (L, w ) - Q(«,£)-Q(^ i) 


qui presentent la symetrie croisee. 

Si maintenant P et Q presentent la symetrie croisee, Texpression 



P 

F- s ‘-i 




supposee int^grable, presentera une troisieme espece de symetrie 
que j’appellerai la symetrie inverse, c'est-a-dire que S et U chan- 

geront de signe quand on changera sen^ ou bien w en ~ 

Dans ces conditions, 

dU Q 

dz ~~ F s ^ 1 z 


change de signe quand on change w en ~ * Si nous posons 


dXJ 

dz 



z m , 


notre p^riode polaire, si elle existe, sera ; elle sera ind£- 

pendante de w (qui joue ici le role dey). D’autre part, elle devra 
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changer de signe quand on changera w en elle est done nulle. 

c. Q. F. D. 

Examinons main tenant le degre de ces diff'erents polynomes; 
mais nous ne definirons pas le degre de la meme maniere quedans 
les numeros precedents. 

Soit un polynome en z. cp, — , 

Y A z a w b . 


Nous dirons qu’il est de degre m si 

| a I -r- | b j % m. 


Soit alors h le degre de H, nous voyons que celui de P et de Q 
est h — i et que celui de S est h — 2. 

Un polynome de degre h contient y A 2 -h 2A + 1 coefficients ar- 
bitrages, mais, s’il est symetrique, il n’en contient plus que 

^ 1 2 * Si s 7 m & vie est directe et, en effet, les coefficients 

de z ±a iv^ b se deduisent d’un seul d’entre eux, de sorte que nous 
n’avons plus a envisager que les termes ou les deux exposants sont 
nuls ou positifs. Si la symetrie est croisee, il arrive (pour P par 
exemple) que les termes independants de z sont nuls, de sorte que 
nous ne devons plus envisager que les termes ou l’un des expo- 


sants est nul, et l’autre nul ou positif. Cela fait 


h(h-L-i) 


coefficients 


distincts. Si enfin la symetrie est inverse, les termes ind^pen- 


dants soit de pp, 

soit de 

z, sont nuls, et les deux exposants doivent 

£tre positifs, de 

sorte qu 

’il reste — 

2 

— coefficients. Nous sommes 

done conduits au Tableau suivant : 


Polynome. 

Degr6. 

Symetrie. 

Nombre des coefficients. 

H 

h 

directe 

(A + i)(A + 2) 

2 

P 

h — t 

croisee 

h(h — 1) 

2 

Q 

h — 1 

croisee 

h(h — 1) 

2 

s 

h — 2 

inverse 

(A — 2)(A—3) 
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Le nombre des expressions D distinctes est done 

( h -f- 1 ) ( h h- 2 ( k — 2 ) (h — 3 ) . 

L 1 — h( h — 1 ) = 4 . 

2 2 

Done, entre les coefficients dudeveloppement de la fonction 
perturbatrice, ily a des relations lineaires de recurrence dont 
les coefficients sont rationnels par rapport aux elements et qui 
permettent d 3 exprimer tons ces coefficients a l 3 aide de quatre 
d 3 entre eux. Chacun d 3 eux } considere comme fonction de Vun 
des elements , satisfait d une equation differentielle lineaire 
du quatrieme ordre. Considere comme fonction de tons les ele- 
ments ; il satisfait a un grand nombre d 3 equations aux derivees 
partielles. Tons les coefficients et ioutes leurs derivees par - 
tielles peuvent s 3 exprimer lineairement a l 3 aide de quatre des 
coefficients , ou bien encore a Uaide de Cun d 3 eux et de trois 
de ses derivees partielles. 

297. On arriverait au meme resultat en conservant les variables x 
ety; on aurait alors 

A 2 = a 2 -+- a ' 5 -+- aa’ u ( - H- — ) -+- aa’v ( xy H — ) • 

\y * 1 \ xyJ 

Ce polynome serait de degr£ 1 au sens du n° 290 et non plus du 
n° 296. II serait d’ailleurs symetrique au sens du n° 291 et non 
plus du n° 296. Enfin l’expression x A 2 se reduirait, pour x = o, a 

aa'f+f), 

qui n’est pas un carre parfait; nous pourrions done appliquer les 
conclusions du n" 291, et, puisque 

/= i: w = o, 

le nombre des expressions II distinctes serait 

4(/n- <*0 2 = 4. c. Q. F. D. 

298. Les equations aux derives partielles aux quell es conduit 
l’analyse pr£c£dente ont d£j& et 6 rencontr^es. Ce sont celles que 
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nous avons etudiees an n 1 ’ 261. Nous avons vu que les onze deri- 
vees partielles de Z, designees dans ce numero par ( 8 ). ( 9 ), ( 10 ), 
( 11 ), ( 12 ). (i3'h s'expriment lineairement en fonctions des 
neuf expressions t j 4 ) qui sont elles-memes liees par les trois rela- 
tions (i5). II y a done, entre ces onze derivees, cinq relations. 
Soil 

Z = ^ U z h w k . 

Jmd 

Les coefficients de z h w k dans les onze derivees partielles seront 

M ^ w ,, TI , nl *u dn) # u 

dx 7 dp 7 dv 7 U ’ A ' ’ dp dv 7 dA 7 

d* U d*L d* U 

d%- 7 dec dp ? dz dv 

II y aura done cinq relations entre U et ses neuf derivees des 
deux premiers ordres. D’ailleurs, outre les relations (i5), les neuf 

expressions (i4) sont liees entre elles et a Z = ^ par la 

relation 

(H-a 2 )Z 7 — ‘ixpTd cos \ — ^avZ'cosir} = — sZ . 

Nous aurons done six relations entre U et ses neuf derivees, 
ctest-a-dire que CJ satisfait a six equations aux derivees partielles 
du deuxieme ordre. De U et de ses derivees, en tout dix expres- 
sions, il en reste done quatre qui sont independantes. 

Parmi ces six Equations, nous distinguerons celles qui ont 6t6 
etudtees a la fin du n° 261 et dont nous avons vu la parente avec 
les Equations de M. Appell. Nous avons vu que ces Equations ad- 
mettent quatre solutions distxnctes, de sorte que U consider^ 
comme fonction d’un seul element satisfait a une equation diffe- 
rentielle du quatrieme ordre. 

299. Les fonctions que nous avons etudiees dans ce Chapitre 
sont d^finies par des integrates doubles et ces integrates doubles 
sont prises le long des contours fermes. En d’autres termes, ce 
sont des piriodes de l’integrale double indefinie 

is 
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Supposons que nous fassions varier Fun des elements, que 
j'appellerai a, en lui donnant bien entendu des valeurs imagi- 
naires, et que cet element x decrive dans son plan un contour 
ferine. Notre coefficient, represente par notre integrale double, est 
une fonction analytique de Felement a; quand cet element aura 
decrit son contour ferme, nous retomberons sur une autre deter- 
mination de cette fonction analytique, et cette determination ne 
pourra etre qu'une autre periode de notre integrale definie double. 

Supposons que cette integrale definie ait k periodes fondamen- 
tales 

Pi, P 2 , ..., Pa; 

toutes les autres periodes et, par consequent, toutes les determi- 
nations de notre fonction, seront des combinaisons iin^aires a 
coefficients entiers de ces periodes fondamentales. 

Quand done la variable decrira un contour ferme, les diflerentes 
determinations de la fonction subiront une transformation lineaire 
a coefficients entiers. Considerons maintenant plusieurs fonctions II 
qui pourront differer par les polynomes H et Q, mais pas par le 
polynome F, toutes ces fonctions subiront la meme transforma- 
tion lineaire, puisque celle-ci ne depend que de la deformation 
des contours dfintegration. 

Soit done A le determinant forme par k determinations corres- 
pondantes de k fonctions II; ce determinant sera multiple sim- 
plement par un facteur numerique quand la variable decrira utx 
contour ferme; le rapport de deux de ces determinants demeurera 
done constant, ce sera done une fonction uniforme. 

Done entre k + i fonctions II, il y aura une relation lineaire 
dont les coefficients seront des fonctions uniformes des elements. 
Doric le nombre des periodes fondamentales est 6gal a celui 
des fonctions II en fonction desquelles toutes les autres peuvent 
s 3 exprimer par des relations lineaires dont Les coefficients sont 
des fonctions uniformes de s elements . 

Si les fonctions II ne pr^sentent que des singularites alg£- 
briques, ces fonctions uniformes se comporteront comme des fonc- 
tions rationnelles. Les coefficients de nos relations seront des 
fonctions non seulement uniformes, mais rationnelles. C’est ce 
qui arrive quand Q = o. 

Nous pouvons prdvoir par la que le nombre des periodes fonda- 
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mentales est de 16 dans le cas general et de 4 si ies excentricites 
sont nulles. (Test ce que montre le resultat obtenu au snjet du de- 
veloppement suivant les anomalies excentriques. 

Mais alors nous pouvons en tirer une conclusion relative au de- 
veloppement suivant les anomalies moyennes. II y aura entre les 
coefficients de ce de veloppement des relations lineaires de recur- 
rence dont les coefficients seront des fonctions non rationnelles 
mais uniformes des elements, et ces relations permettront de les 
exprimer tous a l 7 aide de seize d’entre eux. De plus chacun de ces 
coefficients satisfera a une equation differentielle lineaire, qui sera 
du seizieme ordre (et non plus du trente-sixieme), mais dont les 
coefficients seront non plus rationnels, mais uniformes. 

Ces relations de recurrence et ces equations differentielles sont 
analogues a celles que nous avons rencontrees dans l’etude des 
coefficients de Laplace. Nul doute qu’elles ne puissent rendre les 
memes services dans le cas ou les excentricites sont nulles ; mais il 
n’en est pas de meme dans le cas general; leur ordre semble trop 
eleve ; heureusement il est permis d'esperer que ces equations ne 
sont pas irreductibles et que l’etude des periodes de l’integrale 
double permettra d’en abaisser 1’ordre. 

J’ajouterai que M. Feraud, dans le Tome VIII des Annates de 
V Observatoire de Bordeaux 9 a montre que, par suite decertaines 
symetries, cet ordre pouvait s’abaisser de lui-meme dans certains 
cas particuliers. 
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300. Dans les Chapitres precedents, nous a\ons cherche a 
obtenir le developpement analytique des coefficients a l’aide de 
series procedant suivant Jes puissances des excentricites et des 
inclinaisons. L’emploi de ces developpements ne presenterait 
aucune difficulty si le nombre des termes n’etait pas si grand. 
Mais ce nombre a effraye beaucoup d’astronomes qui se sont 
efforces de calculer la valeur numerique des coefficients sans 
passer par Tintermediaire de ces developpements. 

Nous avons vu que les coefficients pouvaient se mettre sous la 
forme d’integrales definies. Telles sont les integrales (i) et( 2 ) des 


n os 280 et 289 


(I) 

4^ a f 

( 2 ) 

— 4^ S A -mm'= f f 


dx dy 

Xinyjn' ^ y) 

QE^ dx dy 
x myni' y ) 


On pourrait done calculer les valeurs de ces integrales doubles 
par quadratures m^caniques; nous pouvons ecrire en effet l’equa- 
tion ( 2 ) sous la forme 

( 3 ) 4'X-A. mm r=: f f E-Hmu-hm'u') da dll', 


ou bien encore sous la forme 


(4) 4^A mm - = j'J l-E-Umn-m'n dldl\ 

dTofi Ton tirerait les formules approchees 

n*A mm r= ^ 


(3 bis) 
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et 

< 4 bis ) 


A»> 




l,7TlI-i-m ' /') 


Dans la formule (3 il faut donner, sous le signe a u de 
meme qu 5 a ul , les n valeurs equidistantes 


2( l)lT 


Cela fera done, sous le signe des termes dont le nombre 

sera /z 2 , puisqiril faut combiner les n valeurs de u avec les n va~ 
leurs de u . 

Pour Ja formule (4 bis), on operera de la meme maniere avec 
cette difference que ce n’est plus a « et a u ! , mais bien a l et a l r 
qffil faut donner les n valeurs equidistantes 


2 7Z 471 

n 1 n ’ 


2 (n — l) 7C 


Le Verrier voulant calculer la grande inegalite Pallas-Jupiter, 
e’est-a-dire le coefficient A* 8 _ 7 , a employe la formule (4 bis). II 
aurait pu se servir avec plus d’avantage de la formule (3 bis). 
Mais il fallait la puissance de calcul de Le Verrier pour avoir le 
courage d’aborder une tache aussi dcrasante. Heureusement il y a 
des moyens d’dviter ces quadratures mecaniques, ou tout au moins 
de reserver les quadratures mecaniques pour le calcul, non plus 
d’une integrate double, mais d 7 une integrale simple. Les plus 
importantes de ces methodes sont celles de Hansen, de Cauchy et 
de Jacobi. 


301. Methode de Hansen. — Il s’agit d’obtenir les coefficients 
$mm r du developpement procedant suivant les anomalies excen- 
triques. Il sera aisd ensuite de passer aux coefficients A. mm ' du 
developpement suivant les anomalies moyennes en employant la 
formule du n° 240, ou bien encore de passer au developpement 
special de Hansen en employant la formule du n° 246. 

Prenons alors 1 J expression de A 2 ; e’est un polynome du second 
degre en E ±ZM , de sorte que nous pouvons ecrire 


(5) 


^2= A -h BE*"' -h B'E-*“'-h CE 2 *"'-+- G'E- 2 *"', 
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les coefficients A, B et C dependant seulement de u \ nous yerrons 
ensuite que Ton a simplement 


C = G' = 


de sorte que Cet C' sont independants de u et, de plus, du second 
degre par rapport aux excentricites. Nouspouvons en consequence 
negliger ces termes en premiere approximation et, appelant A* la 
somme des trois premiers termes du second membre de (5) et R la 
somme des deux derniers, ecrire 


(6) 


1 _ _ i R. 3 R 2 
A A 0 2 Aj ‘ 8 A§ 


La serie (6) converge tres rapidement a cause de la petitesse 
de R, de sorte que le developpement de ^ est ramene a celui de 

i i 

a7 I*' ““ 

Or nous pouvons poser 

Ajj = A -b BE*“'-j- B' E~ iu ' = H 2 [ r -h a 2 — 2 a cos (u' — (3 )], 

H, a et J3 etant des fonctions de u. On en deduit 


(7) = 

les b {k] etant les coefficients de Laplace formes en fonctions de a 
paries procedes du Chapitre XVII. 

A, B, B r peuvent etre regardees comme des fonctions connues 
de u\ il en est done de meme de H, a et (3 et, par consequent, 
des b {k) ; on aura alors 




si C hk est le coefficient de E ihu dans le developpement de 




e’est-a-dire si 

7 

(8) a*C** = y ^Lb^F-ii^hu). 
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On calculera Fintegrale (8 ) par quadratures mecaniques a Faide 
de la formule 


(9) 


J hk y , T T ■> C 5 




ou n est un grand nombre ; la fonction qui figure sous le signe ^ 

est une fonction de u ; on y donnera successivement a u les n va- 
leurs equidistantes 


o, 


2^ 


2 ( n — i 


Telle est la methode de Hansen; je me bornerai a renvoyer pour 
plus de details au Chapitre XXI du Tome IV de la Mecanique 
celeste de Tisserand et en particular aux pages 34 1 a 344- 


302. Methode de Jacobi. — Jacobi cherche aussi a former le 
developpement qui proeede suivant les anomalies excentriques ; il 
decompose aussi A- en deux parties A 2 et R, dont la seconde est 
du second degre par rapport aux excentricites et aux inclinaisons ; 
mais il n’a pas recours aux quadratures mecaniques. Sa maniere 
de mettre A 2 sous la forme 

Ag H- R 


n'est d’ailleurs pas la meme que celle de Hansen. Pour bien la 
faire comprendre, je suppose d’abord que Finclinaison soit nulle. 
On trouve alors, en appelant i\ et £ 7 , r{ les coordonn^es des deux 
plane tes, 

A 2 = (£ H- ir\ — — « V) (£ — ir \ — ?-*- 


D’autre part, si Fon pose 

2 e , 2 t' 

e = i e2 ’ e j g'2 1 

il vient 

a 

5 -4- ir\ = i ~ — ^ (E*»— 26 -+- e*E“ f »). 

On en d^duirait Fexpression de £ — nrj en changeant z et — /, 
et Fon aurait ensuite celles de 5'+ — *V en changeant a, s, 

as, zz en a', e 7 , m 7 , w 7 . 
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Si, dans P expression precedents nous negligeons le terme en 
e 3 E“'“, nous commettrons une erreur du second degre et nous 
trouverons 

= ( 26V — 26 s -7- b E*“ — b' E ^ z "' ) ( *ib\ £ f — 2 b x s — b\ E~* w — b\ E-*“') 
en posant, pour abreger, 

b— b'= 7-, b t = -» b\ = — - — — • 


Si alors nous designons par A' cette valeur approchee de A- et 
par R r erreur commise, nous aurons 

A*= Aj-h R, 

et R sera du second degre par rapport aux excentricites. 

Si Pinclinaison n ? est pas nulle, nous pourrons developper sui- 
vant les puissances de Pinclinaison et dans le developpement ne 
figureront que des puissances paires de Pinclinaison. 

Si done nous negligeons Pinclinaison, Perreur commise sera du 
second degre. Nous pourrons done conserver la meme expression 
pour Aj et Perreur R = A- — Aj sera encore du second degre par 
rapport aux excentricites et aux inclinaisons. 

Nous pouvons ecrire ensuite 

— H 2 (l — Y 0) — Y J£— i(u’— (O') j (j Y JZ+Uu'—Ui') ^ 


en posant 


H = 


v E —/t0 = — - R/(ar— ©') — ~ , 

T a' H -£ 2 ~ b ' 9 


y'E 1 '®' = 2 £ r — 2£ T7 • 


II est clair d’ailleurs que nous pourrions modifier les coeffi- 
cients H, y, to, y ; , to', qui figurent dans A^, pourvu que les modi- 
fications soient seulement du second degr^ par rapport aux excen- 
tricit^s et aux inclinaisons. IPerreur R resterait du second degr£. 

Jacobi profite de cette latitude pour faire disparaitre certains 
termes de R, a savoir les termes tout connus, les termes en 


cos 


{u — u f ). 


cos 
. u. 
sin 
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Je ren\errai pour les details de l'analvse a la Mecartique celeste 
de Tisserand ( t. IV, Chap. XVIII, p. 3oi a 3 06 >. 

Quoi qu'il en soit, nous avons 

1 __ 1 1 R 3 R 2 
A Ay 2 A'f 8 Ay 


de sorte que i R se reduisant a un nombre fini de termes) le deve- 
loppement de ^ est ramene a celui de jh - Si nous posons, pour 
abreger, 


Qi u—u'— co) — ^ £-/ u — ti)') = T -J j 


il vient 


— = — y; — y'-GJ ) —J> ( I — ~ Y'w -1 ) _ 


et il s’agit d'effectuer le developpement suivant les puissances 
entieres, positives ou negatives de z et de to. 

Il est aise de verifier que le coefficient de z A m* sera, au facteur 
pres y h y' k el a un facteur constant pres, une serie hypergeome- 
trique de deux variables de M. Appell par rapport ay- et y’ 1 . 

Ces series sont tout a fait analogues a celles que nous avons 
etudiees au Chapitre XVIII, seulement elles ne se reduisent pas a 
des poljnomes. 

Jacobi ne se sert pas de ces series, qui n’etaient pas encore con- 
nues de son temps: son analyse est un peu differente; on la trou- 
vera dans la Mecartique celeste de Tisserand (t. IV, p. 3o6 
a 3i 1). 

Les coefficients du developpement de ^3 > ^ > * ' * sont li^s 

par des relations de recurrence, et ces relations permettent de les 
exprimer tous a l’aide de quatre d’entre eux (de sorte qu’elles 
deviennent pratiquement utilisables). Il suffit, pour s J en con- 
vaincre, soit de se reporter k ce que nous avons dit au n° 264, soil 
d’appliquer les principes du Chapitre XXI, en remarquant qu’ici 
co = o, f = 1 et que le polynome pr^sente une symetrie parti- 

culiere, puisqu’il ne change pas quand on change z en ^ et to 


303. Method© de Cauchy. — Cauchy, corame Jacobi et Hansen,. 


P. — II. 


TO 
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cherche d'abord le developpement suivant les anomalies excen- 
triques pour en deduire, par le moyen des fonctions de Bessel, le 
developpement suivant les anomalies moyennes. Je nai pas a 
revenir sur le passage d’un developpement a 1'autre; je m’occu- 
perai done simplement de la maniere d'obtenir le developpement 
suivant les anomalies excentriques. 

La methode de Cauchy se rapproche egalement de celle de 

Hansen par un autre point. Cauchy commence par developper ^ 

suivant les anomalies excentriques de la seconde planete sous la 
forme 

=2 B n 'E ln ' u '. 

Les B /2 ' sont des fonctions de u et, si nous posons 
B „-=2 

il vient finalement 

Cauchy calcule les coefficients B /2 ' par des procedes analyliques 
et il en deduit ensuite les coefficients B /2/ / par quadrature meca- 
nique a Taide de Fint^grale 

(10) f B^E ~ lnu du. 

La difference pro vient surtout de la maniere d’obtenir les coef- 
ficients B/ 2 ',* reprenons la formule (5) : 

( 5 ) A 2 = A + BE m -h B'E-*“'h- CE 2 *“'~h C'E- 2 *'“'. 

Si nous £galons A 2 a o, nous obtiendrons une Equation du qua- 
trieme degr£ eny = E* v qui peut s’^crire 

B' 

{1 1) A By H— H— Gy 2 *4- C ; y 2 = o. 

Discutons cette equation ; nous observerons d’abord que ses coef- 
ficients ne sont pas r^els, ou du moins A, C, G sont r^els, mais B 
et B^ sont imaginaires conjugu^s; liquation ne change pas quand 
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on change y en — et i en — i. Les raeines peuvent done se repartir 
en deux paires : 

a =1 E zw , $= l 

ou a et {3 sont reels et plus petits que i, de telle sorte qu’on per- 
mute les deux raeines d’une meme paire en changeant j/- en — et i 
en — i. De plus, 

sy'2 p'-2 

C = C'= 

4 

de sorte que le produit des raeines est egal a-fi, ce qui donne 


11 en resulte que nous pouvons mettre A 2 sous la forme d’un pro- 
duit de deux facteurs : 

A 2 — H 2 [i — 2 a cos (u ' — to) -4- a 2 ] [i — 2 (3 cos( w'-f- to) -f- j 3 2 ], 

H etant facile a calculer quand on connait a, [3 et co. 

Nous remarquerons ensuite que, si Ton neglige e' 2 , le degr<§ de 
liquation s ? abaisse, car C et (7 s’annulent et liquation ( 11 ) de- 
vient simplement 

A -s- By -4 = 0. 

J y 

L’equation du quatrieme degr£ se r£duit done au deuxieme, une 
des raeines etant devenue nulle et l’autre infinie. Liquation se 
r^sout done aisdment pour C^C'^o. Nous pouvons ensuite 
d6v£lopper les raeines suivant les puissances de C, en consid^rant 
A, B, B' et C comme des variables independantes et faisant Q= C. 
Comme C est du deuxieme degr£ par rapport aux excentricit^s, la 
convergence sera tres rapide. Nous vojons en m£me temps que (3 
est une quantity du deuxieme degre par rapport aux excentricit^s. 

II reste a effectuer le d^veloppement de et pour cela on ddve- 
loppera les deux facteurs 

[i — aacos(^~ w) + a 2 ] 2 — ^ j c p Wp u ' , 

[l — 2p COs(tt'-h U>) H- [3 2 ] 2 =s 2 dqW**'. 
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On aura ensuite le coefficient B fi en faisant le produit des deux 
series; on trouve ainsi 


( 

Les quantifies 


c P c n’—p‘ 

s/H ~*p r 

c p E 1 pu, 


ne sont pas autre chose que les coefficients de Laplace 



b[?> 


calcules respectivement en prenanl a = a et a = j 3 . On les deter- 
mines par les procedes du Chapitre XVII. Cauchy prefere em- 
ployer, pour cette determination, les series procedant suivant les 

puissances de ( ou de particulierement commodes 

dans le cas ou p (ou q) sont grands, ce qui est le cas ou il est 
place. 

Quant a la serie 

Op C n ’— py 


elle converge tres rapidement a cause de la petitesse de [ 3 ; en 
elfet, le deveioppement de c ! n ,_ p suivant les puissances de [3 com- 
mence par un terme en done c' n ._ p est du degre 2 \n f — p[ 

par rapport aux excentricites et est par consequent tres petit. 

Ayant calcule ainsi B n ' par la formule (12), Cauchy aurait pu 
oblenir ensuite B n/i ' par la formule (10) et en deduire ensuite A /2/ ^ 
par la formule (12) du Chapitre XVI. Mais ce n’est pas tout k fait 
comme cela qu’il opere, il passe par l’intermediaire d’un d^velop- 
pement mixte procedant suivant Fanomalie excentrique u! et 
i’anomalie moyenne l : 

de sorte que 

J /*2TC ~27U 

I B/i/E — in/ dl= I B n 'E~ £rL/ ( t i — = e pos u ) da . 

0 . " 
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Cauchy calcule O nn ’ a I aide de cette derniere formule et pour 
cela il applique les quadratures mecaniques, c'est-a-dire qu'il 
prend (K etant un entier suffisamment grand ) 

^B n >E- e ” / (r — e cos u), 

B„ et E _z/2 ^ sont des fonctions de u] on donne a u sous le signe 
les K valeurs equidistantes 

it 2 ( K — r ) 7Z 

°, 

On calcule enfin A nr y a l’aide de la formule 

nn ' — Gp,p r n & )i 

analogue a la formule (12) du Chapitre XVI. Cauchy a applique 
cette methode au calcul de la grande inegalite Pallas-Ju piter, en 
faisant n = — 18, n f = y, il n’avait done a calculer que les fonc- 
tions de Bessel de F argument unique 7 eh 

Je renverrai pour plus de details a la Mecanique celeste de 
Tisserand, t. IV, Chap. XVII. 


304 . On pourrait ^videmment fonder d’autres methodes ana- 
logues sur les propri^tes des fonctions elliptiques et sur Femploi 
de formules, telles que celles dont nous avons fait usage au n° 236 . 
Je me bornerai a cet egard a quelques indications sommaires. 
Supposons d’abord les excentricit^s nulles et que nous ayons a 
calculer, par exemple, l’integrale 


<i3) 


Sf 


dz dw z~ h w~ l 


A B f z 4- 


j) 


C ( w - 


ofi A = a 2 -4- a!" 1 , B = aa'v, C = aa'u. Int^grons d’abord par rap- 
port k w; nous avons une int^grale elliptique que nous pouvons 
prendre par le proc£d£ de la moyenne arithm^lico-g£om£trique. 
Ce proc£d£ est fond£ sur F£galit6 


/ 


dw w~~ l 


[/* 


■ b l w 


y = V- 


dw w~ l 


■ b ! ( w h 

w 
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qui a lieu quand 

v/<z 4- 2 6 — a, \J~a — 2 b — 3 ; \J a! 26' = z* — 

\J a! — ‘ib’ — \JH§ . 

Si nous prenons 1 ’integrale qui donne le coefficient de Laplace 

dw w~ l 


liizb 


}“/ 


^/n_ a *_ a («,+ 1) 


les quantites qui jouent le role de y/a + 26 sont 

i+a, 1 — a; 

en prenant les moyennes arithmetique et geometrique, on trouve 


I, yfi — a 2 ; 

en faisant une seconde fois la meme operation, 

s/j — a 2 4 j t > 

, V 1 a 2 ; 

une troisieme fois, 


[ 4- y^i — a 2 \ 2 


> * • * 9 


d’oi 


ou 


"1 =/ 




y/i — a> 


x/*I = — 

V 2 14-/1 — a 2 

ce qui est la formule du Chapitre XVII; nous avons vu quelle 

approximation donne cette formule quand a est plus petit que ^4 * 

Appliquons la m£me m^thode a l’integrale (i 3 ).Nous pouvons 
la mettre sous la forme 


// 


z~ h w~ l dz dw 


a 2 -+- a ' 2 4- ad ^w 4- L 
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en posant 


en posant 


— B — -1 J ±iC = ( a rz cl f- — a 





I-2C 

/ A+bG--! 

SC 

V *) 



Notre integrate, par Papplication de la regie precedente, se 
reduira done a Pintegrale simple 


r 8 i'jzz ~ h dz 

J 


ou a! et oc sont des fonctions de z definies par les formates prece- 
dentes et que Pon pourra calculer par quadratures mecaniques, ou 
mieux de la maniere suivante : 

Comrae B est tres petit de Pordre du carre de Pinciinaison, on 
pourra developper la fonction sous le signe J suivant les puis- 
sances de B [z -f- ce qui nous donnera en meme temps te d£ve- 

loppement de cette fonction suivant les puissances entieres posi- 
tives et negatives de z. 


305. On pourrait faire quelque chose d’analogue dans le cas 
general; il s’agit toujours de calculer Pintegrale 



dx dy 

x mym f ^/R(x, y) 


On commencera par exemple par integrer par rapport hy; Pinte- 
grale & calculer estalorsune integrate elliptique que je puis ecrire 


/ 


A y~ m ' dy 

— a )(r — f*)(.r -r)(r — °) 


et il faut calculer Pune des p^riodes de cette integrate elliptique. 
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Seulement A, a. 3, y. 3 sont des fonctions de x . Quand e est nul, 
urte des quatre raeines a, 3, ■*', o est nulle et une autre infinie. On 
est done ramene au calcul de Tintegrale 


( O) 


ym dy 


\ f r'y — 7 ) 0 ' — 


Le calcul des raeines y et 3 est alors immediat et Ton peut appli- 
quer directement les procedes du Chapitre XVII et en particulier 
ceux du n° 2ofL Toutes les integrales peuvent d’ailleurs se deduire 
de deux d’entre elles par les relations de recurrence du Cha- 
pitre XVII. 

Si e ] n’est pas nul, le cas est plus complique ; il faut d'abord 
determiner les raeines a, {3, y, 3; nous avons explique au n° 303 
comment cela pouvait se faire, grace a la petitesse de e r . D’autre 
part, nous n'avons plus coin me au Chapitre XVII a calculer Tinte- 
grale 



— e 3 ) m du, 


mais Tintegrale plus compliquee 


(ioj 



du , 


ou a et b sont des constantes et m un entier positif ou ndgatif. Et 
en eflfet, pour ramener Tintegrale (i4) a la forme canonique, il faut 
poser 


p(u) etant la fonction de Weierstrass. 

La periode de Tintegrale (i5) depend de celles des quatre int6- 
grales suivantes : 


fdu, f 


^ {u ± u 0 ) du, 


j [£(m 4- w 0 ) — ■+• m 0 )] du. 



( U H- Ui ) — £ ( 24 — W! ) ] du , 


Uq et Ui ^tant d^finis par les Equations 


p(u o ) = >( «i ) = b. 
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On s’en assurerait en decomposant la fonction doublement perio- 
dique ( I en elements simples. 

Or ees quatre integrales ont pour periodes 

27 tU 4 4^1 Ml- 

On a vu au n° 2 o 6 comment on pouvait calculer co< et r i{ ; on 
trouvera dans F excellent ouvrage de M. Schwarz ( Formules et 
propositions pour Vemploi des fo nations elliptiques, d'apres les 
lemons de Weierstrass, Paris, Gauthier-Villars, 1894) des for- 
mules tout aussi convergentes pour le calcul de u 0 et de u K 
(pages 67 a « 3 ). 

Onvoitqu’ici les formules de recurrence permettent d’exprimer 
toutes les integrales en fonction non plus de 2, mais de 4 d’entre 
elles. 

D’autre part 2co 4 , 2 r H , 4 r n Uo, 4^i^i ne sont plus ici des 
constantes, mais des fonctions de x , et il faut maintenant les deve- 
lopper suivant les puissances entieres, positives et negatives de x. 
Ce developpement peut se faire, soit par quadrature m^canique, 
soil par des procedes analjtiques ainsi que je l’ai expliqud a la fin 
du numero precedent. 

Nous remarquerons qu’il est preferable de faire la premiere inte- 
gration en prenant pour variable non pas r comme nous Pavons 

fait plus haut, mais il arrive alors que nos quantites u><, 7p, etc. 

varient peu quand on fait varier x , les variations etant de l’ordre 
des excentricit^s. 

306 . Seulement ce n’est pas la fonction perturbatrice elle-m£me 
qui figure dans les equations differenti elles du mouvement; ce sont 
les derivees partielles de cette fonction si Pon emploie la methode 
de la variation des constantes ; ce sont les composantes de la force 
perturbatrice avec d’autres m^thodes. 

Si nous avions les coefficients de developpement de la fonction 
perturbatrice sous la forme analytique, il serait aise d’en deduire 
les coefficients correspondants dans Je developpement des d^riv^es 
partielles, ou des composantes de la force., Mais il n’en est plus de 
mdme si nous possMons seufement la valeur num^rique de ces 
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coefficients, et c’est la tout ce que nous donnent les methodes 
exposees dans le present Chapitre. Cela nous oblige a examiner la 
question a ce point de vue nouveau. 

Pas de difficulte en ce qui concerne la derivee partielle suivant 
Tune des anomalies mojennes l ou l'. Si l’on a 

on trouve immediatement 

^ \ ~ S im 

Pour les autres derivees partielles, il s’agit (a etant Pun quelconque 
des elements) de calculer 

d_ £ _ JP 

doi A A 3 
OU 

P== _ i. ff*!. 

2 dec 

Si Ton veut les composantes de la force suivant trois axes rectan- 
gulaires, on aura a developper 

6-r ZLZL 1 7 kzl 

A 3 9 A 3 ’ A 3 5 

en designant par ?1 7 , V, 'Q les coordonnees rectangulaires 

des deux planetes. 

Dans lamethode de Hansen ( voir Tisserand, t. IV, Chap. XXI, 
p. 340 011 considere les composantes de la force suivant trois axes 
particuliers et l’on est amend a envisager les combinaisons 

A 3 * A 3 

Dans tous les cas, il s’agit de developper une expression de la 
P i 

forme -r? et ~ est encore de la meme forme en faisant P = A 2 . Ce 
A 3 A 

qui nous int^resse c’est que, quand on d^veloppe P suivant les ano- 
malies excentriques, e’est-a-dire suivant les puissances entires, 
positives et negatives dex etdejp, ce developpement ne contiendra 
qu’un nombre fini de termes. 

Cela est vrai des coordonnees rj, t /, ZJ et par consequent 
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de leurs combinaisons 



Cela est vrai de A 2 et par consequent de 
II faudra done operer de la facon suivante : 

i° On developpera suivant les anomalies excentriques. Les 

methodes des numeros precedents qui visaient particulierement i 

s’appliquent sans changement a et en particulier a -i • 

2° On multipliera le developpement de ~ par celui de P. 

Ce dernier developpement se reduisant a un nombre fini de 
termes, cette multiplication ne presente aucune difficulty. 

p 

3 ° On passera du developpement de — suivant les anomalies 

excentriques a celui de cette meme quantite suivant les anomalies 
moyennes a Paide de la formule (12) du Gbapitre XYI. 

Une derniere observation toutefois. Les derivees dont 

dx A da 

nous venons de parler et que Pon rencontre dans la mdthode de la 
variation des constantes, sont les derivees prises par rapport a un 
systeme de variables parmi lesquelles figurent les anomalies 
moyennes. 

Designons au contraire par avec des d ronds, les derivees 

prises par rapport a un systeme de variables parmi lesquelles 
figurent les anomalies excentriques. Nous avons immediatement A 2 

et il est aise d’en deduire et r 7; mais ce qu’il nous faut e’est 

dx da A * 

d 1 

^ -* Tout d’abord, si Pdlement a n’est pas Pune des deux excen- 

tricitds e ou e\ on a simplement 

d 1 __ d 1 
da A da A 

Si ot = e, l = u — e sin w, il vient 

d 1 __ d 1 d 1 idl _ d 

de A de A dl A de de A dl A 
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Nous connaissons le de\ eloppement de et de ~ —• J1 reste done 

a tromer celui de sin;/ — 

dl A 

Soit done 


I =2B,p'E rpu-p-u = Va„„„. E‘ ml+m-n 

avee 

A. turn' ~ C A p pt 5 C = ^—^7 3 , n — p( me ) 3 m ' — p r ( m e ), 
il viendra 


d_ i 
oe 


I =- V —.IUL E Hpu+p-u-) = V 

A Jmmi. Oe 


D'a litre part 


a^ ec 


II reste done 


avee 


et 


i-w ~2d C ~de 


d_ r 
de 

dA m 

m' 


L — y dA mm ' Tpn mt +_ m ' n 
A de 


, X7t> dC 
~de~~ ““ 2 d "^T" 2j pp ~dZ' 


sinu Ti i =2 G " iw ' E,,,/n/w/,> 




rfC 


<^C pp r , 

= -j^r !/«-/> (^) J, (/n e ). 


La formule pr^eddente est analogue a la formule (ia) du Cha- 
pitre NVL Nous connaissons et ; les formules pr(§c4- 

dentes nous permettent d’en d^duire A mm, &mm', Gmm' et par con- 
sequent le d^veloppement de j su ^ vant les anomalies mojennes. 
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307. On peut etre amene a rechercher une valeur approchee du 
coefficient A mm > quand m et m ! sont des entiers tres grands en 
valeur absolue et cela pour deux raisons : 

i° On peut se rendre compte ainsi de la rapidite de la conver- 
gence des series; 

2 ° t n terme d'ordre eleve A mm peut donner lieu a une pertur- 
bation notable si, le rapport des moyens mouvements n et rt ! etant 
sensiblement commensurable, le diviseur mn -f- m! n' devient tres 
petit. Le coefficient de la perturbation de la longitude est alors de 
I’ordre de 

/// m ' 

( mn -+- m' n' )- 7 

et il convient alors de calculer le coefficient A mm * sans avoir besoin 
des coefficients precedents. Le plus souvent, une fois ce calcul 
termine, on reconnaitra qu : il etait inutile, parce que la valeur 
trouv^e pour A mm ' est trop petite. II y a done inter^t a se faire 
d 'avarice une idee de Fordre de grandeur de ces coefficients, et 
m^me a pouvoir en trouver rapidement une valeur approchee. 
G’est ce but qu’on peut esp^rer atteindre par Pemploi de la m 6- 
thode de M. Darboux pour le calcul des fonctions de tres grands 
nombres. Cette mdthode repose sur les principes suivants : 

i° Si une s^rie 

?(■*>= 

est convergente dans un cercle de rayon p\ et que sur la circonffi- 
rence de ce cercle elle poss&de un point singulier tel que dans 
le voisinage de ce point la fonction ©(^)'soit developpable suivant 
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les puissances entieres, T'eelles croisscintes et d'ailleurs fraction- 
naires ou meme incommensurables, positives ou negatives, de 

i — —3 de telle facon que le premier terme du developpement 
dont Fexposant ne soit pas entier soit 

A 

on aura approximativement pour n tres grand 

A 71 s - 1 

an= rTT? 

2° Si le point singulier est logarithmique, c’est-a-dire si 'f(z) 
est de la forme 

?{-> = log ( l_ 5^) 

<p 1 (5) et cp 2 (s) etant developpables suivant les puissances reelles et 
croissantes de 1 — — ? il faudra opererde la fagon suivante. Soient 

Z Q> 

s { et $ 2 les exposants du premier terme du developpement de cp< 
et de dont Fexposant n’est pas entier. Soit s 0 l’exposant du 
premier terme de les exposants entiers n’etant pas exclus. Si 
Si<s Q , on pourra appliquer la formule precedente en changeant s 
en Si . Si Si^s 0 , il faudra envisager le premier terme de cp 2 



Si s 0 n’est pas entier negatif, on a approximativement 

A 0 n s 0— 1 log/i 
if ~ Vis,) ? 

et, si s 0 est entier negatif ou nul, 

a n — i)*o+i r(i — s 0 )n*o-i ; 

z o 

3 ° Si la fonction cp(s) presente plusieurs points singuliers sur la 
circonf^rence du cercle de convergence, la valeur approximative 
de a n sera la somme des valeurs approximatives partielles que l’on 
obtiendrait en consid^rant isolement les divers points singuliers; 
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4° Supposons que s(V) soit de la forme 

et que cette serie double converge dans une couronne comprise 
entre les deux cercles | s |= p < , [ J3 1 = p 0 , p t > p 0 - Alors la serie 

cL n z a convergera pour |s[<pi et ne presentera aucun point 
singulier sur le cercle |s|=p 0 ni a l’interieur. Au contraire la 
serie ' S ^J) n zr n convergera pour |^|>p 0 et ne presentera aucun 

point singulier sur le cercle |js|=p 4 ni a l’exterieur. La valeur 
asymptotique de a n s’obtiendra d’apres les regies pr^cedentes en 
envisageant les points singuliers de 0 ( 5 ) sur le cercle |js|=p 4 ; 
celle de b n s’obtiendra d’apres les memes regies en envisageant les 

points singuliers de © (z ) , ou pi utot de © ^ j sur le cercle | z j = p 0 . 


308. Comment ces principes peuvent-ils s’appliquer au pro- 
blem e qui nous occupe? II semble d’abord qu’ils sont faits unique- 
ment pour les fonctions d’une seule variable. Aussi M. Flamme 
a-t-il commence par decomposer le terme a evaluer en une somme 
dont chaque element etait le produit de deux facteurs dependant 
chacun d’une seule variable, qui etait l’anomalie moyenne de la 
premiere planete pour le premier facteur et celle de la seconde 
pour le second facteur. 

Mais on peut operer d’une autre maniere. Soit 
( 1 ) m = an -+- 6, ml — cn •+• d. 


a, 6, c, d sont des entiers finis et donnes une fois pour toutes; 
n est un entier tres grand; a et c sont premiers entre eux. II faut 
calculer 


Nous avons 


— 4tc 2 A 


A mm' — k-an-hb, cn-i-d • 

E& Qd&dy 


■/; 


x m ytn' y / R(a ? 1 y) 

n me / 1 \ ml e' / 1 \ 


avec 



i6o 

Soit 
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ae / i , ce' ! i \ 

be / i \ de' / i • 


Notre integrale deviendra 


U) 


r r / EQq \ n E StQdxdy 

J J \x a y c ) x b y d y/R (x,~y) 


Considerons la fonction auxiliaire 


*(*; = - 


■ 4 


y^-A, 


Sous le signe ^ nous ne donnons pas a m et m! toutes les 

valeurs, mais seulement celles qui sont de la forme (i); mais nous 
pouvons encore operer de deux manieres : 


i° Nous pouvons donner a n toutes les valeurs entieres posi- 
tives , outre la valeur zero; nous obtenons ainsi V integrate de 
Feraud 


( 3 ) 


*(*) = 


si 


E^iQ dx dy 



2 ° Nous pouvons donner a n toutes les valeurs entieres positives, 
negatives ou nulles; la fonction 4>(x?) peut alors se mettre sous la 
forme d’une integrale simple. Nous pouvons trouver deux entiers. 
a et y tels que 

a a -+- c v = i, 


puisque a et c sont premiers entre eux. Posons alors 
E i{=z<*tr c , E ir z=z ztt u . 

Nous avons le ddveloppement 
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Considerons alors Pintegrale simple 

( 4 ) j' i fbc-ad - 1 z ~{zb+-ydj c h 

prise le long de la circonference |*j=i. Cette integrate peut 
s’ecrire 


W dt-> 

X = i[m — b ) -+- p(m' — d) 1 

jjl = a( m! — d) — c ( m — b ) — i . 


L’integrale f tv dt est nulle a moins que p. = — i et alors elle 

est egale a 2 irz; pour que [jl = — i, il faut que m et m ! soient de 
la forme ( 1 ); mais alors onal = /i; de sorte que Pintegrale (4) se 
reduit a 


Ainsi, pour resoudrele probleme qui nous occupe, nous n’avons 
qu’a rechercher La position et la nature des points singuliers de la 
fonction $( 5 ), et pour eviter toute confusion nous distinguerons 
la fonction <h(s) de M. Feraud definie par Pintegrale double (3) 
et la fonction <£> 0 ( 3 ) definie par Pintegrale simple (4)* 


309. Appliquons ces principes au calcul des coefficients de 
Laplace qui sont donnas par la formule 


avec 


*ziTzb ( s k) = J F ~ s z k ~ l dz , 

F =(1 — a«)(i- 2). 


Considdi'ons d’abord s comme fixe et faisons croitre k , il s’agit de 
trouver le coefficient de z k dans le d^veloppement de F 


F-* = (i — a z)-* 



Cette fonction pr^sente deux points singuliers, Pun sur le cercle 
ext^rieut k la couronne de convergence, l’autre sur le 


P. — r* II. 


II 
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cercle interieur, ^ = a; si nous supposons k positif et tres grand, 
c’est le premier qu’il faut considerer. 

Dans le voisinage de ce point, on a sensiblement 

F-j? = (1 — a -g ( 1 — a- , 

ce qui nous donne asymptotiquement 

b$ k) — ~ — a ; qi — a 2 )-* k~ s . 
s T(s) v J 

Supposons maintenant k fixe et s tres grand ; soit 



Alors 2 ir. b\ h '' sera le coefficient de $ n dans le developpement de 
Pintegrale 



Les points singuliers de la fonction sous le signe j nous sont 
donnes par 

F = 0 , F = (3. 

II faut ou que ces deux equations aient Line racine commune, ce 
qui donnerait (3 = 0, solution a rejeter, ou que l’une d’elles ait une 
racine double. Nous pouvons exclure la premiere qui ne depend 
pas de (3; il reste la seconde, qui a une racine double si 

z = =h 1 , (3 = (1 dr a) 2 . 

La racine qui convient c’est (1 — cl)' 2 . 

Quelle est la nature de ce point singulier? Pour [3 tres voisin 
de (1 — a) 2 , les parties les plus importantes de Pintegrale seront 
celles qui correspondent aux valeurs de z voisines de 1 . On a 

alors sensiblement 5=1, a, de sorte que Pintegrale 

s’ecrira sensiblement 



La foncLion sous le si 



est une fonction rationnelle dont il 
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faut ohtenir le residu; or ce residu est 


‘fi—; 


^ elant donne par l'equation 

Mais cette equation donne sensiblement 




— 1 ) i — r p 

? 

a 


Le residu est done 


(r — a)v/a 


W(i — a ) 2 — p 


et 1’int^grale est egale a ce residu multiplie par 2 /tt. Done b { s k) est 
le coefficient de J3" dans le developpement de 


ai/a [’ 0-«)*J 


_ J_ 

* 


e’est-a-dire 


-if* 


2/a(i — a) 25 — 1 p 


formule qui, on le remarquera, est independante de k. 


310. La determination des points singuliers de la f’onction <3>(s) 
ne presente pas de difficuitb; on n’a qu’a appliquer aux integrates 
(3) et (4) les procedds du n° 282. Pas de difficulte non plus en ce 
qui concerne la nature de ces points singuliers, je me bornerai & 
renvoyer aux Methodes de la Mecanique celeste } t. 1, p. 3^2. 

La difficulte provient de ce fait que tous les points singuliers ne 
conviennent pas et n’appartiennent pas a la branche considerbe de 
la fonction <£(,s). Nous avons vu en effet aux n os 282 et 283 quelle 
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est la condition pour qu’une singularity de la fonction <£(s) con- 
vienne. Cette singularity se presente quand deux points singuliers 

de la fonction sous le signe j se confondent et, pourque la singu- 
larity comienne, il faut que ces deux points singuliers soient, 
avant de se confondre, depart et d'autre du contour ddntegration. 

Les points singuliers qui conviennent sont dits admissibles et 
il nous faut choisir, si nous adoptons Fintegrale ( 3 ), celui des 
points singuliers admissibles dont le module esl le plus petit, et, si 
nous adoptons Fintegrale (4), celui des points singuliers admissibles 
dont le module est le plus voisin de i. 

La discussion pour reconnaitre F admissibility des points sin™ 
guliers est assez delicate et a ete jusqu’ici le principal obstacle a 
Femploi de cette methode. J’ai indique dans les Methodes nou- 
velles de la Mecanique celeste les principes generaux qui per™ 
mettent de faire cette discussion. Mais je n’ai applique ces prin- 
cipes qu’au cas ou Finclinaison est nulle, Tune des excentricitys 
nuile et Fautre tres petite. M. Hamy, dans le Journal de Liouville 
(1894 et 1896), a traite le m£me probleme sans s’astreindre a la 
troisieme condition. M. Cocuiesco, dans sa these, 1895, s’est occupy 
du cas ou Finclinaison est nulle, les deux excentricites petites et 
differentes de zero et ou la longitude du pyrihelie est la meme. 
Enfin M. Feraud, dans sa these, 1897, a traite le cas ou Finclinaison 
est finie et les deux excentricitys nulles; il a d’ailleurs retrouve les 
rysultats de M, Hamy. 

Il semble d’abord que la discussion doive etre plus facile avec 
Fintegrale simple ( 4 ) qu’avec Fintegrale double (3); il n’en est 
rien, a moins que l’une des deux excentricitys ne soit nulle, parce 
que la fonction sous le signe j n’est pas une fonction uniforme 

de Jf, mais possede une infinite de determinations, de sorte que la 
discussion ne pourrait se faire que par la considyration d’une 
surface de Riemann a une infinity de feuillets. Aussi, M. Fyraud, 
en introduisant Fintygrale ( 3 ), a-t-il realisy un syrieux progres. 
Mais il faudrait faciliter la discussion des intygrales doubles et 
pour cela ytudier les propriytys de leurs pyriodes. Nous avons dyja 
vu aux Chapitres XX et XXI Fimportance que pourrait avoir cette 
ytude. 

Si 7 au lieu d’envisager le dyveloppement suivant les anomalies 
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moyennes, on envisageait le developpement suivant les anomalies 
excentriques, le problerae serait considerablemenl simplifie. On 
reconnaitrait, par exemple, qae la fonction <E>(s) satisfait a une 
equation difTerentielle lineaire a second membre dont les coeffi- 
cients et le second membre sont des fonclions rationnelles de 

Je n'insisterai pas davantage sur cette question, me bornant a 
renvojer aux Alemoires cites et en particulier au Cbapitre VI 
des Methodes no uve lies de la Mecanique celeste . 
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